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FOREWORD 


This document was prepared by Rocketdyne Division, Rockwell International Cor- 
poration, in accordance with Article I and Line Item Nos. 2, 3, and 4 of the 
Data Requirements Lise of Contract NAS9-14315 with the National Aeronautics and 
Space Administration. The contract period of performance was 6 August 1974 to 
1 September 1975. The NASA/ JSC Technical Monitor was Mr. F. D. Freeburn. The 
Rocketdyne Program Manager was Mr. R. H. Helsel for the first three months; he 
was replaced by Mr. R. D. Paster for the remainder of the program. Mr. J. A. 
Nestlerode served as the Principal Engineer, assisted by Dr. D. R. Kahn. 

Several technical people at Rocketdyne performed work or served as consultants 
regarding specific areas of the various program tasks: Mssrs. J. K. Hunting, 

R. L. Nelson, and L. E. Sack with respect to the feed system hydrodynamics, Mr. 
F. R. Linow with respect to combustion dynamics, Mr. M. D. Schuman with respect 
to combustion dynamics, chamber dynamics, engineering model formulation, and 
computer programming, and Mr. K. W. Fertig with respect to numerical analysis, 
computer programming, and checkout. 


ABSTRACT 

This report is an operating manual for the Feed System Coupled Stability Model. 

It is submitted as partial fulfillment of an 11-month program designed to de- 
velop, verify, and document a digital computer model that can be used to anal- 
yze and predict engine/feed system coupled instabilities in pressure-fed stor- 
able propellant propulsion systems over a frequency range of 10 to 1000 Hz. 

The first section describes the analytical approach to modeling the feed system 
hydrodynamics, combustion dynamics, chamber dynamics, and overall engineering 
model structure, and presents the governing equations in each of the technical 
areas. This is followed by the Program User's Guide, which is a complete de- 
scription of the structure and operation of the computerized model. Last, ap- 
pendixes provide an alphabetized F0RTRAN symbol table, detailed program logic 
diagrams, .computer code listings, and sample case input and output data listings. 
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INTRODUCTION 


Historically, during the development of pressure-fed propulsion systems, feed 
system/engine coupled instabilities have been frequently encountered. Resolution 
of these problems usually included increasing injector pressure drop to decouple 
the feed system from the combustor, the result being substantial system weight 
penalties. A dynamic computer model would be a useful tool in obviating coupled 
stability problems during the development of the Space Shuttle Orbit Maneuvering 
System (SS/OMS). A model could be used both as a system design tool to optimize 
component location and pressure profile (minimize system weight) and a system de- 
velopment tool to define test programs for assessing stability margins of the OMS. 

This document is an operating manual for the Feed System Coupled Stability Model 
(FSCSM) and is submitted as partial fulfillment of an 11-month program conducted 
by Rocketdyne to develop and verify an engineering digital computer model for the 
NASA/ JSC which can be used to analyze feed system/engine coupled instabilities in 
pressure-fed, storable propellant, propulsion systems over a frequency range of 
10 to 1000 Hz (frequencies lower than the chamber transverse frequencies). The 
model is sufficiently general so that it may be readily applicable to present and 
future engine and propulsion programs. For scaling purposes, the baseline config- 
uration chosen is the OMS engine. The model has been written for use on the NASA/ 
JSC Univac 1110, EXEC-8 computer system, and provides NASA a tool which can be 
used to: 

1. Conduct preliminary design tradeoff for feasibility studies prior to 
propulsion concept selection. 

2. Guide the design of propulsion systems to ensure stability at all oper- 
ating ranges and with minimum penalties. 

3. Guide testing programs by predicting the least stable operating regimes 
thereby reducing the number of stability tests required. 

4. Provide stability verification in the event system changes are made and 
hot-fire verification is impractical. 

5. Diagnose problems on existing systems and evaluate potential solutions. 
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The work performed in completing the requirements of the program's technical ef- 
fort is described in a separate companion document, entitled OME Feed System- 
Coupled Stability Model, Final Report (Ref. 1), This includes the mathematical 
formulation of the model, development of the model into an overall engineering 
structure, and verification of the model's operation and capabilities by compar- 
ing the model's theoretical predictions with experimental data from an QMS engine 
and test rig with known feed system/engine chugging history. 

The present document contains a detailed description of the structure and opera- 
tion of the FSCSM. In the first section, the mathematical formulation of the 
model is reviewed. The analytical approach to modeling the feed system hydrody- 
namics, combustion dynamics, chamber dynamics, and overall engineering structure 
is described and the equations utilized by the model in each of the technical 
areas are presented. The reader may consult Ref. 1 for more details pertaining 
to the derivation of the equations. 

The Program User's Guide section contains the instructions necessary to operate 
the computer model and interpret the results. First, the structure and logic of 
the main program and all subroutines are described, followed by a description of 
the input data required to run FSCSM. The input is divided into four major sec- 
tions: (1) main control, (2) nozzle admittance control, (3) hydrodynamics con- 

trol, and (4) combustion dynamics control. The format, content, and description 
of each input data card is clearly specified for each control section. The output 
of the FSCSM computer program is then discussed in terms of each tabular page of 
printout. Finally, additional details on program operation are presented, includ- 
ing program size, overlay structure, computer time, and program input/output data 
set file information. Appendixes provide an alphabetized FORTRAN symbol table, 
detailed program logic diagrams, computer code listings, and sample case input 
and output data listings. 
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MATHEMATICAL FORMULATION OF MODEL 


INTRODUCTION 

During certain periods of a rocket engine's operation, conditions 
within the combustion chamber and feed system are time variant, i.e., 
the operation is not steady with respect to time. Prime interest of 
this computer model is focused on abnormal transient operation during 
unstable combustion, i.e., pressure oscillations in a combustion device 
which are driven by the feed system and sustained by the combustion 
process. Start and stop transients are not considered. 

The deviations from steady-state combustion which occur durina unstable 
burning depend upon the kind of instability experienced. Liquid rocket 
instabilities are classified according to their dominant time-varying 
processes. They may be divided initially into two categories, depending 
upon whether the instability oscillation wave length is long or short 
compared with the chamber dimensions. 

If the instability wave length is considerably longer than the chamber 
length and diameter, pressure disturbances propagate rapidly through the 
combustion space compared with rates of change due to the instability. 

As a result, wave motion in the chamber may be neglected and chamber pres- 
sure can be considered to vary only with time but not to vary spatially 
(i.e., is a lumped parameter). These instabilities depend upon a fluid 
mechanical coupling between the propellant feed system(s) dynamics (fluc- 
tuating injection rates), the propellant combustion rates (delay times), and 
the combustion gas exhaust rates (pressure relaxation). Such instabilities 
can be further subdivided into various categories depending on the extent 
of wave motion in the feed system. 

The breakpoint at which chamber wave motion becomes important is not abrupt. 
In reality, chamber wave motion is always present and, in effect, lumped 
chamber- instabil i ties are really "zero order mode" limtts of more general 
wave motion instabilities. In practice, it is found that the chamber gases 
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can be considered to act as a lump until the frequency of oscillation 
exceeds roughly one-fourth of the frequency of the lowest chamber 
acoustic resonance mode. At and above such frequencies wave motion 
becomes important and cannot be neglected in analysis. Chamber wave 
motion instabilities are characterized by the wave-length of the oscil- 
latory motion being comparable to the chamber dimensions. As with 
lumped chamber instabilities, the driving energy comes from oscil- 
latory spray combustion. With wave motion instabilities, however, in 
addition to the effects of injection rate fluctuations, there is the 
combustion response of burning propellant sprays as they are disturbed 
by passage of a pressure wave through them. Wave motion may increase 
local burning rates by any of several mechanisms: (!) a pressure effect 
on the drop vapor gas phase burning rates; (2) enhanced mixino between 
gases and between sprays and gases; and (3) increased spray gasification 
rates. Increased spray gasification may be due to transient increases 
in convective flow velocities, to increased temperature or concentration 
gradients, and/or spray droplet shattering. The instability amplitude 
depends upon the magnitude of the response, and vice versa; typically, 
the interacting processes are driven to a limit represented by abrupt, 
essentially complete consumption of the propellant sprays. This direct 
response can be so great that injection rate fluctuations may be of 
secondary importance. As a result this class of instability can also 
be further subdivided as to the importance of feed system coupling. In 
the absence of feed system coupling, the instability is referred to as 
classical acoustic instability." Only longitudinal chamber modes with 
feed system coupled instabilities are considered in this program. 
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FEED SYSTEM DYNAMICS 


Development of the Waterhammer Equations 

Consider the differential control volume of a fluid element in a duct 
shown in Fig. 1. 



Fig. 1. Differential Pressures Developed Across the 
Incremental Length of a Fluid Element 


Fluid compressibility and Newton's second law leads to the following 
pair of differential equations: 


There are several ways in which to solve these equations. The solution 
method presented here follows that of Ezekiel (Ref. 2). The general 
form of**the solution that satisfies either of equations^!) and ( 2) i 

P = F, (t + f) + Fj (t - f) 
where F-j and F^ are arbitrary functions. 


_e _ ^2 3v 

3t 3X 3X 


= 

3X 


3V 


3t 


B 3V 
_2 3t 


where 


p = fluid pressure 
V = fluid velocity 
B = fluid bulk modulus 
p = fluid density 
c = acoustic velocity = 
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Taking the partial derivative of p with respect to x and t separately 
and substituting the results in equations { 1) and ( .2) gives: 


where 


3V _ 

1 

i£ = 

1 

3X 

6 

3t 

e 

3V _ 

_ i 


_ 1_ 

3t ~ 

P 

3x 

pc 

\ 

F' 

' (?) 

'^1 

II 

(0 

3? 

ixpression 

for 

V is ( 


The expression for v is obtained from either equation ( A) or ( 5); 
zv = - Fi (t + |-) + (t - 


where 


i- 


Z = ( p 6 )'* . 


( 4) 
( 5) 


( 6 ) 
( 7) 


Letting the subscript o denote x-0, the upstream position, and the subscri 
L denote x=L, the downstream position, and defining t = L/c as the signal 
propagation time between the two positions, equations ( .3) and i .6) become 

Po = + F^Ct) 

Pl = Fi(t + t) + F2(t- t) 

2Vq = - F^(t) + F2(t) 

zVl = - F^(t + T) + F2(t-t) 

Combining Eqs. { 8) and ( 10), and Eqs. (9) and ( 11) separately, yields 
four additional relations: 

Po-:^ = 2 F2(t) 

Po - ^''o = 2 F^(t) 


Pt 


( 8 ) 

( 9) 

( 10) 
( 11 ) 


( 12 ) 
( 13) 
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= 2 F2(t-T) 


- zVj^ “ ^ • 


( 14) 
( 15) 


Eliminating the functions and F2 gives the final result as: 

[pL - = Po - -0 


Consider now Fig. 2, which depicts a generalized line segment forming a 
portion of a feed system with many such segments. 


''n- "n- ”n 


w 




1 
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Flow 


w 


n+1 

n+1 


Figure 2. Generalized Line Segment 


The equations which describe the pressure and flows as functions of time 
and of each other for the generalized line segment are obtained from 
Eqs. (16) and (17): 


( 18) 


( 19) 
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The expression, t - indicates values at t seconds before, 
and " 

P’n+1 = Pn+1 ^ Pn'^K 


( 20 ) 


"'n = Pn ^ ''n • 


( 21 ) 


Equations ( 18) and ( 19) are solutions of the wave equation, and equation 
( 20) is the flow through a nonlinear fluid resistance. Letting 


°'n ~ A g 

n ^ 


(22) 


these equations can be combined to give: 

Pn “n '^n “ [Pn+1 ^ *^n ^'^n+l ^ ^'^n+1 " “nO, . 

(t-Tn) 

P. * «n-1 * “n-l «„ ' [p„-1 * Vl Vl]„ 


(23) 


(24) 


Eliminating and rearranging into quadratic form results in 
Vl * (°„-l * - [P„_, + Vl “n-l] 


* [Pr.+1 * ''n!"n*ll'"n+l ' “ni] 


'‘■Vl' 




= 0 


which can be solved for the appropriate solution using the quadratic 
formula. The tank end parameters are obtained using a solution of 
Ep- ( 23) only. The injector end solution is obtained using the 
quadratic formula for equation (25). 


(25) 


The linear model incorporated in the Hydrodynamics subprogram utilizes 

the same basic equations, (23) and (24 ), but in the following linearized 
form: 
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R = 


2R^ 



(31) 


At the tank end, the term 6 is zero for constant tank pressure. At the 
injector end, 6 p^^-| is the independent variable. 
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Injector Dynamics 

The injector dynamics are included by treating the injector as a lumped 
compressible volume as shown in the figure below. 



Figure 3. Schematic of the Injector 
as a Lumped Compressible 
Volume 


The pressure in the injector manifold, pj, is related to the entering 
flow, from the upstream pipe segment and the injector flow, Wj, 
as follows: 

<lPl _ '^I^ . . 

dt“ " K ■ Wi) (32) 


where Vj is the injector volume and Cj is the fluid sonic velocity. 


The injector flow is controlled by the differential pressure across the injector as 
well as by the resistance and inertia of the injector orifices. Thus, 


.2 

Pi - Pc = Rj Wj + 


Jl a ■ 

Ag dt '^I ’ 


(33) 


where p^ is the thrust chamber pressure, Rj is the injector hydraulic 
resistance and £/Ag is the equivalent inertance of all the injector 
orifices combined, i.e., 

n 

^/Ag " ^ • 


(34) 
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In the preceding equation, and are the length and area, respectively, 
of an individual injector orifice. 


An additional factor which can have a significant effect on the response 
of the feed system to chamber pressure oscillations is injector face 
flexibility. This effect can be expressed as a change in injector volume 
proportional to a change in injector pressure drop: 


d V, /d p, d 

dt \ dt dt 


Also, 

^ _ c^ ^ 

dt ~ g dt 



which can be rewritten as 


2 2 

^ _ C CP 

V? -~W 


V . 


(35) 


(36) 


(37) 


Combining Eqs. ( 35 ) and ( 37 ) gives 


d pj 

dt 


^ I 

^'^n - ) - ~Y^ 


V 


dPi dp^' 

3t“ " W~ 


(38) 


which can be rewritten as 


1 + 


Cj p j K 


dp 


I _ 


2 

Cj pj K dp^ 


ar = ^«n - ^ “ Vj g ar 


(39) 


This expression reduces to Eq. (32 ) if no injector flexibility exists 
(K = 0). 
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Two-Phase Flow Acoustic Velocity 

In the waterhanmer equations the acoustic velocity, c, of the fluid 
appears in two places; (1) directly in the constant relating flow to 
pressure, and (2) indirectly in the time delay value, t, which equals 
Vc seconds, where ^ is the pipe segment length. The acoustic velo- 
city of a fluid is a property of that fluid. However, its effective 
value can be reduced by the elastic walls of the flow passage or by 
the entrainment of gas and vapor in the liquid (two phase flow). 

Gas in the liquid can appear from two sources. One source is direct 
entrainment from mixing of gas and liquid in* the propellant tank, 
while the other can result from the evolution of dissolved gas as 
the pressure drops along the feed system. 


Given the steady- state pressure at each point in the feed system and 
data on the solubility of the pressurant gas in the propellant as a 
function of pressure and temperature, the amount of gas in the fluid 
can be determined for each feed system segment. Then, knowing the 
amount of gas in the liquid, the effective acoustic velocity of the 
mixture may be calculated. 


Assuming isentropic compression of the gas, the change in volume of the 
gas is 

and for the liquid 



V 

Defining a constant, a = 

the following relation is obtained: 




g 

a 


Kp 
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The bracketed term is the compressibility of the mixture. The density of 
the mixture can be shown to be 


“Pq + Pjl 
‘^m" (1 + a)' 


The acoustic velocity of a liquid in an elastic pipe is 


c = 


1 


e eE 


(43) 


(44) 


Using the above expressions for density and compressibility, the acoustic 
velocity, can be written as 


c = 


1 


Pm 

/a 1 1+a 

Dcf\ 

T+^ 

\ Pg V ^ 





(45) 


This expression can be used to define the acoustic velocity of a feed 
system segment with two phase flow. For an all liquid system, a = 0 
and the same equation can be used. 


In th^^Hydrodynamics subprogram the effect of the wall compressibility 
term, on the fluid acoustic velocity is handled automatically 

t)Cf , 

(assuming input value of are provided for each feed system segment). 

However, the program does not compute the effects of two-phase flow. If 

such flow occurs in the feed system being modeled, an effective fluid 

acoustic velocity must be pre-calculated for each affected segment. 

Dcf 

Equation ( 45 ) above, with the term set equal to zero can be used 
for this calculation. 


Simulation of Branch Lines 

In the Hydrodynamics subprogram, branched lines are handled by assuming that 
each branch has zero internal volume and that its flows are incompressible. 
Thus, the pressures at the end of all segments which meet at a branch are 
set eqijal . The continuity of flow is then used to pro^de the additional 
equations in combination with the waterhammer equations to solve for the 
overall feed system dynamic response. 
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Generalized Feed Systen Model 


A schematic of the generalized feed system which is modeled by the hydro- 
dynamics subprogram is shown in Fig. 4. The configuration chosen is 
based on design and operating mode data for the OHS, PBK, and RCS feed 
systems obtained from McDonnell Douglas/St. Louis. The system is comprised 
of 30 individual line segments, each denoted in Fig. 4, as the lines be- 
tween the black dots. A continuous parameter representation of each line 
segment is obtained through the use of separate sets of waterhammer equa- 
tions. Each line segment can have a different line length, area, wall 
compliance, fluid acoustic velocity and resistance, and hence can model 
a wide variety of feed system components by merely choosing the appro- 
priate values from these parameters. Also included in the generalized 
model are lumped parameter descriptions of two injectors (designated 
"0" and "F" on Fig. 4). Parameters for the injectors are volume, re- 
sistance, inertance, fluid acoustic velocity and face flexibility. 

The system of 57 equations describing the generalized Fig. 4 feed system 
is listed in Table 1. The equations are shown in the linearized, LaPlace 
transformed format required by the frequency response subroutine. 
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Figure 4. Generalized OME Feed System Schematic 
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TABLE 1. HYDRODYNAMIC EQUATIONS 


N , -2T,s. 

+ (R^ + a^) + [P^ + (R^ - a^) vv^J e = 0 1-1 

^1 ■ ®3'^3 ■ ^^2 ^ ^*^3 ■ ^3^ ® 1-2 

-T-s 

?2 + (R 3 + a^) - [P^ + a^w^] 8 = 0 1_3 

“2T 5 

P] - (R2 ■*■ ^2) (w-| - w^) “ [P] - (P'2 ■ 32) (w^ - w^)] e ^ = 0 1-4 

-TqS 

^2 ■ ^5^6 ■ *^'^3 ® 0 1-5 

-T,s 

P2 - ^4 ^'^4 ■ '^6^ ■ ^*^4 ■ ®4^ ”15^ e ^ = 0 1-6 

Ps (^4 + w^5 - [P2 + a^ (w^ - Wg)] e ^ = 0 1-7 

P3 + (Rg + dg) W7 - [P^ + 3g Wg] e = 0 1-8 

-■■■fiS 

P 3 - ag W 7 - [P 4 + (Rg - ag) wg] e =0 1-9 

P4 - (R7 + 37) (wg - w^q) - [F^ - (R^ - 37) (wg - w^q)] = 0 1-10 

^4 ■ ^8 '^10 ■ ^^8 “ ^8^ ”11^ 8^=0 1-11 

*^5 ^*^8 ^8^ '^11 ■ '*' ^8 ”10^ 8 ® = 0 1-12 

P5 ~ ^9 '^ 14 -“ [-Pg ■*■ (Rg “ 3 g) W^2l 8^=0 ^ 1-13 


R-9808/16 



TABLE 1. (Continued) 


-TgS 

^6 ^ ^ ^9^ *^12 ” ^^5 ^ ^9 ''ll^ ^ ~ 0 


”^10^ 

^6 ■ ^10 ”l 2 ■ ^^7 ^^10 " ^ 10 ^ '^ 13 ^ ® ® 


“^IqS 

^7 ^*^10 '*' ^ 10 ^ '^13 ■ ^*^6 ^10 '"^ 12 ^ ® 


^7 “ V '^14 ■ “^0 '^14 ^0 


^^7 - *13 " l^gvy„ *14 " ''o ='’7 = [gM % ^’’0 

“^1 

^8 ^ ^14 '^17 “ ^^12 ■ ^*^14 ■ ® 14 ^ '^ 24 ^ ^ ° 


-T14S 

^12 ■ ^*^14 ^14^ '"^24 ■ ^^8 ■ ^14 ‘‘'17^ ® ~ ° 


T^^s 


^8 ■ ®11 ^'^15 ■*■ '^17^ ■ ^^9 ^*^11 ■ ^11^ '^18^ ® ° 


-T-j 1 s 

Pg + (Rn + a^^) w^3 - [Pq + (w^5 + e = 0 

•^■^13" 

Pg - (R ]3 + 3 ^ 3 ) (w^g - W^g) - [Pg - (R ^3 - 3 ^ 3 ) (w^g - W-jg)] 6 =( 


"T'lzS 

^9 " ^12 ''^19 ” ^ ^12 ’ ^ 12 ^ '^ 20 ^ ^ ~ ^ 


■2T^5S 

^n * ^*^15 ^15^ '''22 '*’ ''’ ^^15 ■ ^15^ ”22^ ® ~ ° 


-TigS 

^11 ■ II6 '^22 “ ^^12 ^^16 ■ ^ 16 ^ ^”24 * ^ 25 ^^ ® 


1-14 

1-15 

1-16 

1-17 

1-18 

1-19 

1-20 

1-21 

1-22 

1-23 

1-24 

1-25 

1-26 
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TABLE 1. (Continued) 


^12 ^ ^ ("^24 ^ ''25^ ‘ ^Pn + a^g e = 0 1-27 

P^3 + (R22 + 822) W34 + [P-,3 + (R22 - 822) VV34] e = 0 1-28 

P-,3 - 324 W3g - [P^7 + (R24 - 324) W37] e = 0 1-29 

2 T 

^13 ■ ^*^23 ^23^ ^^^34 ' ”36^ ' ^^^13 " ^^23 ‘ ®23^ ^”34 " '^36^^ ® 

P12 - (R^7 + 3^7) W25 - [P^4 - 3^7 W2g] e = 0 1-31 

Pl4 + 8 i 7 «2g - [P^2 " (^17 ■ 3^7) W2g] e =0. 1-32 

'"ly ■ ^21 ^”37 " '^38^ ■ ('^21 ■ ^21^ ”27^ e = 0 1-33 

Pl4 + (R21 ^21^ '''27 ■ I^’^17 ^21 ^'^37 ' '^38^^ e = 0 1-34 

Pl4 - a^g (w2g + W27) - [P-|g + (R^g - a^g) W29] e = 0 1-35 

^15 ^'^IS ®18^ ^29 ■ f^^l4 ^ ^18 ^^26 *^27^^ e = 0 1-36 

Pl5 - (R^9 + a^g)(w2g - Wg^) - [P^^ - (R^g - a^g) ( w 29-W3T )] e =0 1-37 

^15 ~ ^20 *^31 ■ (^20 ■ ^20^ '^'32^ e = 0 1-38 

*^16 '*’ ^*^20 '*’ ^20^ ”32 " ^^15 ''' ^20 '^31^ ® ® 
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TABLE 1 . (Continued) 


■^25^ 


^17 ■ ^25 '^38 ■ ^ ^'^ZS " ^25^ '^39^ ® ° 


^18 ^^25 ^ ^25^ '^39“ ^^^17 ^ ®25 '^38^ e - 0 


■^26^ 


'" l8 ■ ^26 '^39 " ^^19 ^*^26 “ ^ 26 ^ ” 40 ^ ® ^ 


P ]9 + (R^e ^26^ '''40 ■ *^*^18 ^26 '^39^ e - 0 


- 2 T 27 S 


^19 ^^27 ^ ^27^ ^'^40*^42^ ” ^^^19 ” ^^27”^27^^'^40”'^42^^ ^ 


■^28^ 


Pl9 - ^28 w^2 - ^^20 ' ^28^ '"43^ ® ° 


^20 ^*^28 ^28^ '^43 " ^^^19 ^28 '^42^ ® ° 


-T 29 S 


P20 - 929 '''43 ■ ^^21 ^*^29 ■ ^29^ '^44^ ® ° 


^21 ^ ^*^29 ^29^ '"44 ■ ^^20 ®29 '"43^ ® 


-T 29 S 


-■^30" 


^21 ■ ^30 '■'44 ■ l^^22 ^^30 ■ ^30^ ''^45^ ® ^ 


■^3n^ 

^22 ''’ ^*^30 '*’ ^ 30 ^ '^45 ‘ 1-^21 '*' ®30 *^ 44 ^ ® ^ 


P 22 ZpS w^g - Rp w^g Pp 


_ /c2\ . . A2\ 

^ ^22 ; (gV/r ^'■'45 ■ ^ l^/p "^F ^^22 


Igvjp ^''f 


1-40 

1-41 

1-42 

1-43 

1-44 

1-45 

1-46 

1-47 

1 —48 

1-49 

1-50 

1-51 

1-52 
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TABLE 1. (Concluded) 


^ ^ ^ ^6^ '-^8 - tPs " ^6 = 0 

PlO + (R]2 ^ ^12^ '^20 “ ^^9 ^ ^12 ”19^ e = 0 

Pl 7 + (R24 + 324) W37 - [P ^3 + 324 ‘^ 36 ^ e = 0 


1-53 

1-54 

1-55 

1-56 

1-57 


^ 


NOTE: 3jj 


A, ^9 
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It should be noted that the hydrodynamics subprogram solves the complete 
system of 57 equations (describing the complete Fig. 4 feed system) 
each time it is called. Thus the frequency response of the entire system 
is calculated each time. It has been shown, however, that simpler feed 
systems, representing only a portion of the Fig. 4 schematic, can be 
modeled by merely assigning values to the parameters of the unneeded 
line segments which will exclude them from any effect on the system fre- 
quency response. This is accomplished automatically by the hydrodynamics 
subprogram via the assignment of very large resistances and very short 
lengths to all line segments for which no data is entered. 
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COMBUSTION DYNAMICS 
Analytical Approach 

In the past, the combustion response has been modeled with a simple time 
delay(s) (Ref. 3 through 9). This time delay represents the tine required 
for the propellants to travel at their injected velocity from the point 
where they are injected to another point where they burn, and implies the 
burning is concentrated at a fixed nlane some arbitrary distance from the 
injector face. The procedure outlined above is obviously an oversimolifi- 
cation of the burning process which is distributed in some fashion through- 
out the combustion chamber. 


Steady-state combustion models (Ref. 10 and 11 for example) provide insight 
to determine the droplet burning distribution as well as additional infor- 
mation required to relate the distribution to a combustion response as a 
function of frequency. Combustion models are designed to march incremen- 
tally down the combustion chamber from a set of specified initial conditions. 
In so doing, the model calculates the rate at which the propellants are 

consumed as a function of the axial position in the combustion chamber (burn- 
ing rate profile). 


The analytical technique selected to describe the corrfcustion dynamics is 
based on employing the mathematical expressions used in the steady-state 
combustion models (in particular the JANNAF DER program. Ref. 11). These 
mathematical expressions are expanded into time averaae and oscillatory 
components and are described in the following sections. 


Atomization Process 

A very essential part of the combustion field initialization is the assign- 
ment of propellant spray droplet sizes and flowrates. Analytical descrio- 
tions of the atomization process are not available but empirical correla- 
tions that relate droplet diameter to injector geometry and flow conditions 

are available (Refs. 12. 13, and 14). For like-doublets, one empirical 
relationship is (Ref. 12). 




4.85 X 10^ 


Vj-°-^' (Pc/Pj)'°-'' d. 


0.57 


(46) 
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where Vj is the liquid jet velocity and dj is the liquid jet diameter 
at the atomization plane. (For steady-state analysis, the velocity is 
the injection velocity and the diameter is the orifice diameter.) 


For purposes of the current analysis, the atomization process is des- 
cribed by: 


D. =K(d.)® (V 

“ ^ „ J x=x. 

imp imp 


where x^.^p is the location of the atomization plane or the impingement 
point. Expanding Eq. 47 into time-averaged and oscillatory parts, yields 
the oscillatory droplet diameter 



In order to evaluate the oscillatory droplet diameter, the oscillatory 
liquid jet diameter and velocity (and therefore the jet flowrate) are 
required at the atomization plane. Therefore, the dynamics of the fluid 
from the injector to the atomization plane is required and outlined in 
the following section. 

Klystron Effect 

The dynamics of the liquid propellant jet from the injector face to any 
location in the chamber are described by the continuity and momentum 
equations: 


Assuming 

p. = constant 
J 


({)=(}) + $ 


(({> any variable), 
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where 


(time average value) 
(oscillatory value) 


(54) 


<t> = f(x) 
$ = 

<P'= g(x). 


(55) 

(56) 


the preceding egeations can be expanded into time average and oscillatory 

parts and integrated between the injector face and any location in the 
chamber to yield: 



where « is the angular freguency and the oscillatory injection rate. (sT.) 
is determined by the feed system dynamics. Equation 59 is the oscilla-^ 
tory jet flowrate at x and is usually referred to as the Klystron effect 
(Ref. 15 ). The Klystron time delay, is therefore given by 

Xu 


inj 


T 



(60) 
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Considerable amplification of the injector face flow oscillations are 
possible when the Klystron effect is present and could explain the periodic 
burst of acoustic resonances called resurging and the steep-fronted 
waves seen in low and intermediate frequency instabilities. 

Droplet Vaporization 

Theories of droplet combustion (Refs. 10, 16, 17) are available which 
may be used to evaluate the extent of coupling between droolet burning 
rate and local pressure and velocity fluctuations. In general, droplet 
burning is enhanced by increased turbulence levels or by periodic direc- 
tional variations in velocity, because droplets are relatively heavy and 
resist following gas streamlines. 

Calculation of the spray heating and vaporization is usually accomplished 
through specification of the corresponding individual droplet processes 
and suiranation over all the droplets that constitute the spray (s) being 
analyzed. The calculation of single droplet evaporation is usually based 
on a spherically symmetric model of simultaneous heat transfer and mass 
transfer across the gas side boundary, or film, separating the liquid 
droplet from the surrounding hot combustion gas. Forced convection and 
resultant nonspherical transfer processes are accounted for through 
empirical Nusselt number correlations for both heat and mass transfer. 

For the fuel or oxidizer spray, the droplet continuity equation is 


d7 \ ^ap^ 

and the vaporization rate is (Ref. 


(61) 



ttD|^ kf^ Nu,^^ Z,^ 


( 62 ) 


where pj,- is the spray density (mass of spray per unit chamber volume), Nj^ 
is the number of droplets per unit chamber volume, and 
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c_ Nu„ p MW 3f 
\ ^ ''k \ 

■ k- NUm R T, 

h \ h 


P-P. 


(63) 


Noting that 


ttD. 

"k = ^"V = \ -T 


(64) 


the droplet number flowrate can be written as 

A v.P, 


"k ' »k '' "k 


'k k 


m, 


(65) 


Therefore, Eg. 61 can be written as 


ST '"k^k* ' - "vap|, 


( 66 ) 


For steady-state combustion models, the preceding equation (along with 
Eg. 62) is numerically integrated allowing the droplet diameter, D|^, 
to vary along the length of the combustion and maintaining constant 
droplet number flowrate (N^). Combs (Ref. 18) has shown that changing 
from a variable droplet diameter to a variable droplet number flowrate 
yields approximately the same results for steady-state vaporization. 
Therefore, in order to simplify the integration for stability analysis, 
the droplet diameter was held constant and the dronlet number flowrate 
was assumed to vary. 


Summing Eg. 61 over all fuel or oxidizer droplet size groups yields 


I 






Pk (6)Z, 




(67) 
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which can be written as 




Ap, 


A tn. 


vap. 


( 68 ) 


wnere 




Pc- = 

k k 




S (^k^) 


(69) 

(70) 


_ P|.(6)Z k. Nu„ 

X - xy_L_l__L_k 

^ P. ^k S 

^k Py, 


(71) 


Letting 7^^, , Cp be independent of k and assuming 


k 


Tj = f(t) 


yields; 


r\ Pr- 

, - /s \ ^ 

^s (6)Zg k^ Nu 


s «s 


where 


•"s = (^^•nj 



(72) 


(73) 


(74) 


From Eq. 68 

J(ApsVs) 


dx 

^s"s 


(75) 
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Integrating Eq. 75 between x^, the start of vaDorization olane, and any 
location x yields 


0 


I ^ 

^s'^s 


(76) 


substituting Ep. 76 into Eq. 68 yields the fuel or oxidizer spray 
vaporization rate: 


'"vapj AXjVg 


X 

r dx 


0 s'^s 


(77) 


Using perturbation techniques, the time average vaporization rate can be 
written as 


(m. ) 

s'x=x. 


m 


vapj ATjV^ 


exp 


(x-x^) 

^s-^s 


(78) 


and the oscillatory vaporization rate can be written 


as 




s'x=x 


m = m. 


vap, "Vap^ Tirr: 




S'X=X. 


Xo 


V \ ^ 

^1 Vs j 


(79) 


Assuming 


(Ve)v= 


S'X=X, 


(80) 


yields 


= (=^ ) 

Vs 


x=x, 


(81) 
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Letting p. , c , and 
fore be expressed as 


be constant, the oscillatory time delay can there- 



Nu, 



^"^s X MR 
3MR ' 


The oscillatory spray droplet diameter (D^) is given by Eq. 48 and the 
oscillatory flowrate is given by Eq. 59. The above formulation results 
in a linear oscillatory vaporization model similar to, but more realistic 
than Crocco's n-x model (Ref. 4). The formulation includes the effects 
of: (1) distributed energy release, (2) oscillations in the injection 

rate, (3) oscillations in droplet diameter, (4) oscillations in droplet 
temperature, (5) gas pressure and velocity oscillations, and (6) oscil- 
lations in the local mixture ratio. 


Nusselt Number. It may be observed that one of the dominant terms in both 
the expressions for the average and oscillatory time delay is the Nusselt 
number. The Nusselt number, for longitudinal modes, is (Ref. 19). 

= 2.0 + 0.6 

In ortfer to evaluate the oscillatory Nusselt number, the oscillatory droplet 
spray velocity is required. The droplet spray velocity can be obtained from 
the drag equation. 


pD. 

— |v - V I 

p ' S' 






's dt 


^ pD^ |v-vj (v-v^) Cr 
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Letting 


"s 

r - S 

drag^ " (18} v 


(85) 


where 


o = 


24 


'D. P 


V - V. 


(86) 


the oscillatory droplet spray velocity can be written 


as 


'^s 


1 + i 


0) T 


drag^ 


1 + (w T. )' 
drag^' 


V = P V 
''s 


(87) 


Definina 






v-v. 


-iV2 


caM. 




( 88 ) 


where 


A^s = 


|y-vj 


J Steady state , 


(89) 


the Nusselt number can be written as 


Nu„ = 2.0 + 0.6 Pr 

“s s 


1/3 


-, 1/2 


pD 

— CAM 
P s 


1/2 




(90) 
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Expanding the preceding equation into time average and oscillatory 
parts yields 



F ^ 

if = Rp , (94) 

Fp p p 

the oscillatory Nusselt number is 
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For small perturbations in the pressure, the linear response factors and 
the time average values for and F,, are 

(96) 


(97) 



Calculations have been made which indicate that, for laroe droplet diameters, 
the average and oscillatory Nusselt numbers are Quite sensitive to pressure 
and velocity oscillations. Therefore, the Nusselt number can have a signi- 
ficant effect on engine stability. 


Dr o £ let Heat Transfer Blockage Term . The oscillatory combustion time delay 
given by Eq. 82 requires the evaluation of the heat transfer blockage term 
(Zg) which is related to the combustion gas and liquid vapor properties by 
Eq. 63. Because the vapor pressure (Py ) at the droplet surface is related 
to the droplet temperatore, the blockage^era also depends on the oscillatory 
droplet surface temperature inside the droplet which is given by: 




3 




3 r / 


(98) 


Therefore, the oscillatory heat transfer rate to the droplet can be related 
to the oscillatory droplet surface temperature by 



The droplet heating rate can also be written as (Ref. 10) 

(T - T J 


^ *^f 

s s 


VdD, 


(e 5 - 1) 


s J 


(^D^) 


Assuming that 




0 (droplet at "wet bulb" temperature) 
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(99) 


(100) 


( 101 ) 



and 



etc. for the other variables, the response factor for the heat transfer 
blockage term can be related to droplet and gas properties and flow conditions. 

Examination of the response factor for the heat transfer blockage term indicates 
that this term is not important at low frequencies (also see Ref. 20 and 21). 
Therefore, the oscillatory heat blockage term has not been included in the 
computer program. Detailed eouations for this term are presented in Appendix B 
of Ref. 1. 


Generalized Vaporization Rate Expression 


In order to maintain generality in representing the combustor dynamics, the 
spray vaporization rates (fuel and oxidizer) were written as: 



! .C5 (MR)^^o .c (MR) 

s\ ■■ / s s 

X 






/ 


•Xq u n P /x=0 


s \ P 




-^"13 


x=0 


I P 




dx 


''s 


(103) 


Combining the expressions of the preceding sections with this generalized 
vaporization rate expression yielded the combustion coefficients: 


c 







(104) 


c 




(105) 
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where the subscript s denotes the fuel or oxidizer and 



0 06) 

(107) 

(108) 

(109) 

( 110 ) 

( 111 ) 

( 112 ) 

(113) 

(114) 

(115) 


In the above expressions is the oscillatory injection rate divided 

by the oscillatory pressure al the injector face and is calculated by the 
hydrodynamics suborogram. 

The main function of the combustion dynamics subprogram is the calculation 
of the combustion coefficients. The general soray vaporization rate ex- 
pressions-3re used in the chamber dynamic subprogram whichjs discussed in 
the following section. 
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CHAMBER DYNAMICS 
Analytical Approach 

Two methods of approach were considered for solving the chamber dynamics. 

The first method used a linear lump chamber coefficient. This method is 
valid only at low frequencies (less than 500 Hz) and results in a set of 
nonlinear algebraic equations to be solved. 

The second method employed a first-order perturbation model to define the 
chamber frequency and growth coefficient along with the oscillatory pressure 
distribution in the chamber. This method is valid for all frequencies of 
interest in the present program (10 to 1000 Hz). For the oscillatory 
variables, solutions of the form $ = 4*'s where w is the complex fre- 
quency, were assumed. These forms yielded a set of nonlinear differential 
equations which were numerically integrated between the injector face and 
the nozzle inlet plane. Using iteration techniques and the requisite boundary 
conditions at the injector and nozzle inlet plane, the chamber frequency and 
growth coefficient are obtained. 

Consideration of the degree of complexity in solving the governing equations 
by each of the above methods as well as the range of validity of each approach 
resulted in choosing the first-order perturbation models as the best method 
for describing the chamber dynamics. In the following paragraphs, the derivation 
and solution to the first-order perturbation model stability equations are 
presented. 

First-Order Perturbation Model 

In this section, chamber model equations are stated without showing their detail 
derivations. Complete derivation of the basic equations is presented in Ref. 

1. Assumptions used in the derivation of the basic equations are: (1) ideal 

gas flow is a valid state equation; (2) dilute sprays occupy a negligible 
fraction of chamber volume; (3) the spray can be represented by a finite number 
of dropsize groups; each dropsize group contains a large number of locally 
identical drops; and, each size group constitutes a separate liquid phase and 
exchange terms between liquid phases are not included; (4) drag contributes only 
kinetic energy to the spray energy equation; (5) secondary "shear" breakup of 
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drops is not included; (6) negligible coupling between diffusion and thermal 
gradients; and (7) no body forces. 

The following equations can be formulated for the gas phase: 

Gas Continuity 

V(P(I) = J ^ 
n j 

Gas Momentum 

(PU) + V- (pu; u) 

- 1 s [>; 

n j 

Equation of State 


<"j "vap ) 


= -Vp + V* T 


.n _n n 

- “4 '-I 


(116) 


(117) 


p = pRT 


Shear Stress 



( 118 ) 


(119) 


Gas Energy 

It |^p{h + \[] + V . j^pu (h + ^)^ 


= - 7.q + V. (uTj + || 



(120) 
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Gas Mixture Ratio 


(pMR) + v-(pif MR) 



Heat Transfer Rate 


i 


(p 


‘’eff) 




Drag Force 



( 121 ) 


( 122 ) 


(123) 


Assuming 


(1) Diffuser, thermal and viscous gradients are negligible, 

(2) Droplet drag forces and heat transfer to the droplets are negligible, 

(3) Droplet velocities are approximately equal to the gas velocity, 
ana letting 



(124) 

(125) 


R 


3R 




(MR - MR^) , 


(126) 
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( 127 ) 


m = V V .,n n 
^ox j vap. 


n j 
fu 


m 


vap 


fu 


^ 21 ".-lap. 


n J 


(128) 


where the subscriot $ denotes that the properties are evaluated based on the 
overall injection mixture ratio during steady-state operation, the preceding 
equations can be simplified for longitudinal modes to 


Gas Continuity 

^ It h v) = A (m„,„ + m„_ ) 


Gas Momentum 


p 4+ pv 

^ at ^ ax ax 


''"Pox ''"Pfu 


= 0 


Equation of State 


Gas Energy 


* M ^ '"'If + ' 

[‘"ox - <2 

*”vap,„L‘"fo* (|W)^ 


'"vap 


Gas Mixture Ratio 


„ aMR . 3MR 

^ JT ^ 


= (MR + 1) 1^ 


•"vap - (MR) m, 
^ox 


n 


''"Pfu J 


(129) 


(130) 


(131) 


(132) 


(133) 
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Because of the complexity in solving nonlinear partial differential equations, 
perturbation techniques were used to simplify the governing dynamic equations. 
Assuming 


<() = (|> + (ji {(\) any variable). 


(134) 


where 

♦ = f (x) 


(135) 


and 

= g (x, t), 


(136) 


each perturbation quantity was taken to he of order (e), where (e) is a small 
ordering parameter that is a measure of the wave amplitude. The perturbation 
expressions for each of the independent variables were substituted back into 
the nonlinear partial differential equations, where all terms of the order 
(e^) or higher were neglected. The resulting time-averaged equations were 
solved for the time-averaged variables and the oscillatory equations were solved 
by assuming solutions of the form 

♦ = (137) 

where <!>' = f(x) and o) is the complex frequency. The resulting equations 
form a system of ordinary differential equations in terms of the variables 
4>' and can be numerically integrated by employing boundary conditions 
and iteration techniques. 


Following this approach the perturbation equations were expressed as: 

p H p [l + p' e'''“^3 (138) 
V = V + c^ V’ e"''“^ (139) 
T = T [l + T’ e‘''“^ il-M) 
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|j + p' 


p = p U + p' e 

MR = MR + MR’ 


-io)t 


"'’Pox ' >Pox * "’«Pox " 
"vaPfu ' "vaPf, * "'vap,^ p 


(141) 

(142) 

(143) 

(144) 


The time-averaged equations were determined to be: 


(a) Gas Continuity 

(A p 7^ = ^ ) 


dx 


(b) Gas Momentum 




PPPpX ™Pfu 


(145) 


(146) 


(c) Equation of State 


p = p ^ (wp - «p^)] 


(d) Gas Energy 
dx 


^ * V p ^ <P''> ' 


'V ’> M ™vap, 


OX 


‘V- <5W)x ’ 




'>] 


+ m. 


vap 


fu 




(147) 


(148) 


(e) Gas Mixture Ratio 


p 7 d MR _ ,n,r 


dx 


= (MRr,) 


(149) 
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and the oscillatory equations were determined to be: 


(a) Gas Continuity 


t /-iw\ dv' , v' d(Ap) 


P Ap 

iK.,, ) 

ox 


vaP„v vap^^j 


p c 


0 


(b) Gas Momentum 


^ ^ dv' _ p ' dp 

c , dx cv dx " —'2 dx 

00 0 P 

. p .'_. ^ ^ ..!_ ^ = 0 

-^2dx -^2dx ^ 

P P 


(c) Equation of State 


P’ = P' T’ + ^ MR' 


(d) Gas Energy 


f] 

*aV’ k = 

0 p 




■\ap 


OX 


Q'-OX 


(150) 


(151) 


(152) 
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- mK (2MR.V 


"vap,„ L""fu 


1 ■ ^ ""vap^, ^3 Mr\ 


2 ^ (^)^ MR' 


Gas Mixture Ratio 


MR' (^) Mf ) ^ * [(f) 


p' + v' 


(MR + 1) r^. 

= A_ L 

p c^ ^ Pox 

:r^ L^ao ■ 


- RFT m. 


d W 
dx 


(2 MR + Dm J (MR') 


In the computer model, the preceding set of ordinary differential equations 
are numerically integrated between the injector face and the nozzle inlet plane. 
The method of calculating the complex frequency for the perturbation model, 
based on nozzle admittances calculated from upstream and downstream variables, 
is discussed in the Engineering Model section, 

Steady-State Solution 

The boundary conditions for the steady-state differential equations are 

® X = X,, 


''xq ''x=0 


(if v„ . t 0) 


(Ap V (if f 0) 
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2 

Assuming small Mach numbers, i.e. , M « 1, the steady-state differential 
equations can be integrated between the start plane for vaporization (xq) 
and any location (x) to yield 


p = constant = p 

c 



(Av) 


V = 


P A 


- <C.v) Ah 


inj 


ox' ox 


" '' ^R)x=o ( ImT 


-MR (-|(^) 
^ 3MR ^ 




'i 


(flp v)^.„ + (1 - ♦„.) 


inj 


ox' 


fu 


A V 


(flp (, . 


inj 


ox 


' 0 - %,) ] 


( 160 ) 


(161) 


(162) 


(163) 


T = 


E 

P 

1 + -L / ^ 
lsHR/0 

(MR - MR^) 


(164) 


where 


(}, = e 


(x"Xq)/ Tj 


(165) 
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If the gaseous injection velocity is equal to zero (7 = o). the steady- 

state mixture ratio and density at x = are determined by 



( 166 ) 


(167) 


These equations were developed by taking the 1-imit as x ^ x^. from 
a downstream distance. ^ 


Oscillatory Solution 

The hoe, .dan, conditfons for the oscillatory differential eoeations are 

0 x = 0 
p' = Ap 


(168) 

(169) 

(170) 


From these boundary conditions and the oscillatory differential equations 
the oscillatory conditions at the start plane for vaporization (Xq) can be 
determined and are: 


'0 - 


Ap COS n 


W Xr 


W X- 




„2J20 

e S “ 


(171) 
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where 



If the gaseous injection velocity is equal to zero ~ ®)> oscil- 

latory mixture ratio at is determined by 


MR' 

^0 




^ nj /'It md 


A p c > 

Xq 0 


(1 + mr„ ) 


(MR^ + 1) 




inj 


1 



''fu^ 


I 



m 

m 

vapf^j 


(v' ) MR 
^0 0 


(MR^ + 1 ) 

(MR^ + 2) 
^0 



( 177 ) 


This equation was developed by taking the limit of the mixture ratio 
equation as x Xq from a downstream distance. 
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The ordinary differential equations describing the oscillatory solution are 
solved using a second order implicit finite difference method. This method 
has the advantage of being simple to implement and modify, as well as being 
unconditionally stable for systems of equations which do not have exponen- 
tially grov/ing solutions. The method as applied to the first order system 

Y' = AY + g 

where Y and g are nxl vectors and A is an nxn matrix is as* follows: 


'i+1 




M 

2 




(y,- + y,.+, ) + 


H+h 


(179) 


Here, the subscript i refers to the i'th mesh point in the finite difference 
scheme, e.g., x^. = Xq + iAx. The y. approximate the Y vector at x^-. That 
is y-j ' " Y(x^). The subscript i+Jj refers to evaluation at x. + ax/2; 

e.g.. A^^,_ = A(x. + -y). 

That the above method leads to a second order approximation (error is prooortional 
to Ax^) can be shown as follows: 


Solving for y.^^ yields 

.Vi+l = (I - ^ A.^jJ-'' (I + f y. + AX (I - g.^j, (IRO) 

Without loss of generality, assume i = 0. 


From the two expansions 


’'l ' ''•i * T 


(181) 


'‘o = \ - f K 


Ax' 


■0 '2 '), • 'ii 


+ Y ■ + o(Ax'’) 

2 ^2 


(lot) 


the following are obtained 


Yi = Yq +-AX Y^ + o(ax-^) 


(183) 


and 


Y,. = (Yq + Y^)/2 + o(Ax^) 
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(184) 



3 

Let Yq = Yq; it is necessary to show y-| “ Y-j + o(Ax ) in order to demonstrate 
second-order accuracy. 

From the differential equation 


9.. = V., - A,, Yj, 

2 '2 2 ^ 

Substituting(185)into(lb0)and noting Yq = Yq, 


( 185 ) 


y, = (I - f A^)-' (I + f \) V„ + AX (I - f - V \) 


AX 


ax . A-1 


= (I -f v' 


Yq + AX Y^^ Aj, AX (H Yq - 


(186) 


Using (183) and (134). gives the result 
Yl = Y^ + (I - ^ 


= Y-j + o(ax'^) 


(187) 


Consider now the stability of the finite difference formula (179) for 
systems which do not have exponentially increasing solutions; that is, 
the real part of each of the eigenvalues of A is negative. To prove 
that they are stable for this situation, define the error e- = Y. - y^. 
and consider the two equations given by (179) and 


Y = (T - — A 1“^ (T + — A ) Y + 

^•+l 2 ^ 2 i 


AX (I - f + o(Ax^) ^ 


the latter resulting from 


E 


i+1 




(187). Subtracting (179) from (188) 
(I + ^ 


(188) 


(189) 
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Let B = . The method is stable if and only if the matrix 

(I+B) has a spectral radius less than one, for this would produce (Ref. 22) 

lim [(I-B)"^(I+B)]" = 0 

Since the eigenvalues of (I-B)'^ (I+B) are just equal to (1 +b)/(1-b) 
where g is an eigenvalue of B, the spectral radius of (I-B) (I+B) 

is just 


(190) 



max |(l+3)/(l-3) | 
3 

(191) 

For 

this to be less than one. 



4 - 

A 

1 

(192) 

for 

all 3. This implies 


or 

1+3 + F + 3F < 1-3-J+ 3 3 

(19J) 


3 + 3 < - (3+ D 

(194) 


Real ( 3 ) < 0 

(195) 


Since 3= ^ , where a is an eigenvalue of A, the method will be stable if 
all the eigenvalues of A have real parts less than zero, that is, the 
solutions to(178) are not exponentially increasing. 
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Nozzle Admittance 

The nozzle admittance based on downstream conditions is calculated based on 
the following analysis. 


The gas flowrate of the nozzle inlet plane is 


m = 


P g 


= Ap V 


( 196 ) 


where the characteristic velocity is 


ygv RT 


p 2 -]{y+1)/2(y-1) 

LwrJ 


(197) 


For short nozzles, the oscillatory mass flowrate can be written as 


iV = T' / 8c* 

m "p 7 ^ ^ T \ 3MR 


MR' 


(198) 


Assumi ng 


£l = i £!. II 
7 " ^ p ’ T 


Y-1 


p" 


(199) 


the nozzle admittance for a short nozzle can be written as 


= P c V 


r (^-1) . 

he* \ 

MR’ “1 

L - 

2y p 

\3MR / 

^ P' 


(200) 
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Assuming 



/ 

I 



MR = constant 


( 201 ) 


the nozzle admittance based on downstream conditions becomes 




1- 2Y 

\3 MR / c* p' V ■' 


\r = 


constant 


( 202 ) 




where A^. is calculated using the admittance program developed 

MR = constant 

by-Bell (Ref. 23). 


The nozzle admittance based on upstream conditions is 



(203) 


where (v'/p')^^ _ ^ is calculated from the oscillatory solution. For solutions 
to the chamber dynamic equations, the nozzle admittance based on upstream and 
downstream conditions must be equal. The method of calculating the complex 
frequency which satisfies this condition is discussed in the Engineering Model 
section which follows. 
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ENGINEERING MODEL 


Analytical Approach 

The overall model structure had the greatest variety of factors influencing 
its nature. Some of these factors were related to the overall confidence in 
the success of the effort. Factors relating to cost included the solution 
time and numerical stability, which bears on the number of runs which will be 
required for a solution. Still other factors were related to the JSC Uni vac 
1110 capabilities. The remaining factors concerning the overall model structure 
reflect on its conversion cost applicability to the resolution of propulsion 
system problems. Its accuracy has direct bearing on the design margins 
involved. The type of impact and the obtainability of characterization 
parameters could not limit the accuracy and useability. The type and useability 
of the output was also given due consideration as well as the degree of 
generalization such that the model can be applied to a range of systems. 

The structure of the Engineering Program was based on a trade-off of setup 
time, storage capabilities, and solution time. General input data to the 
program includes geometric factors, engine operating conditions and propellant 
properties. An equilibrium gas properties program similar to NASA ODE 
computer program, and the DER combustion model program are executed external 
to the stability program. The control program then executes the nozzle 
admittance and hydrodynamics programs to calculate the admittance and 
oscillatory injector flowrate as a function of frequency and stores the 
results on tapes. Steady-state distributed combustion parameters calculated 
from the DER Model are inputs to the combustion dynamics subprogram which 
are iterated with the chamber dynamic subprogram until the nozzle and 
injection admittance conditions are satisfied. The solution method for 
obtaining solutions for the complex frequency is outlined in the following 
section. 
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Determination of Complex Frequency 


The complex frequency, uj , is determined such that the boundary condition in 
the nozzle is satisfied. Specifically, the admittance is required to be 
continuous across the interface between the combustion zone and the zone 
immediately downstream of the combustion zone. In the downstream zone, 
the nozzle admittance, , is computed from a nozzle admittance program. 

In the upstream combustion zone, the nozzle admittance, , is computed 
from the oscillatory flow parameters determined by the ^ chamber dynamics. 

The complex frequency must be such that 

\ " \ (204) 


Let w = X + ly and F = A^, - A., 

N|, N 

X and y such that ^ 


u + iv. The numerical problem is to find 


u(x.y) = 0 (205) 

v(x,y) = 0 ( 206 ) 

Several methods were considered for solving this system of equations. 

Because F is not an analytic function of co, the complex form of the Newton- 
Raphson method may not always work. On the other hand, one could use the 
two-dimensional form of Newton-Raphson (Ref. 24), but since the derivatives 
of u and v with respect to x and y must be computed numerically, the two- 
dimensional Newton-Raphson method will require three functional evaluations 
of F at each ( 0 , i.e. , (x,y),(x + Ax,y),and (x,y +Ay). Alternatively, a 
far more efficient method is to use the two-dimensional form of the secant 
method (Ref. 24) since this does not require the evaluation of any deriva- 
tives. Specifically, this method approximates the u and v surfaces with 
linear functions Uj^ and V|^ (planes) based on three previous guesses for w , 
(xi,y-|), (x 2 ,y 2 )» (^ 3 »y 3 ) • The next guess forw, (x^,y^), is determined 
from the equations Uj^ (x^.y^) = Vj^(x^,y^) = 0. The new value of tu then re- 
places one of the previous three values, normally the one with the largest 
error as measjj.red by the absolute value of F(x.,y.), and the iteration 

xJ J 
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is continued until convergence is reached. The actual equations 

for the above process take the following form. Let = Xj + iy^., Uj = 

u(x.,y.), and v. ■ v(x.,yj for j = 1,2,3. 

j j j j j 


1 . 


2 . 

3. 

4. 


5. 


6 . 

7. 


Determine ir., j = 1,2,3, such that 

J 

"'l + ’'2 + ’'3 " ^ 


TTi + TT2U2 + ir3Ui = 0 


T^l V 3 + TT 2 V 2 + ^ 3 V^ = 0 


Compute = 11^0)3 + ''^ 2^2 ^^3^^! " 

Compute u^ and v^ based on x^ and y^. 


(207) 

(208) 

(209) 

(210) 


Test for convergence, i.e., require that min L. - co, l/u.l <e, 

' J 4 ' 4 1 

and IF4I <£2- 

If the process has not converged, continue with steps 5, 6, and 7. 


2 2 

Determine the j between 1 and 3 such that u. + v. is maximum. 

V J 

Replace u^, and v^ by u^, and v^, respectively. 


Go to 1 , 


Operationally, steps 4 and 5 may be altered to replace the w's 
cyclically, i.e., w^_-j,u^. v^. In fact the 

computer program as written alternates between these two procedures 
every three iterations in order to avoid any possible cycling that 
may occur. 

The above algorithm has been found to be very efficienr"when the first 
three guesses are relatively near an actual solution. The difficult 
problem was to develop a searching algorithm which determines the regions 
in the u plane where solutions exist. 
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One possible procedure would be to utilize the fact that the surface u^ + 
has an absolute minimum at each solution. Using any reasonable value of oj 
as a first guess, one might be tempted to employ a gradient, or modified 
gradient, method to march along the surface until one comes near a relative 
minimum. Unfortunately, this procedure fails because the surface u^ + v^ 
has many relative minima which are not actual solutions. The reason for 
the large number of relative minima (and maxima) for this surface is undoubt- 
edly due to the coupling between the combustion processes and the feed system 
oscillation in conjunction with the very rapid change of the feed system 
response as a function of frequency. The searching algorithm must be able to 
discriminate between those relative minima that are not solutions and those 
that are. Such a procedure was developed for this program. It takes 
advantage of the fact that a large portion of the computations required 
are only a function of the real part of w, i.e., they use x as an independent 
variable and do not depend upon y. Thus, y may be changed without having to 
redo many of the calculations within the program. 

Intuitively, the idea is to increment x through a range of values, while de- 
termining y at each x according to the criterion mentioned below, until it is 
determined that a solution has been crossed. This determination employs the 
use of a test function which changes sign when a root is crossed in the same 
manner that a single equation in one unknown changes sign as it goes through 
a zero. The task of developing a defining criterion for y and a test function 
for X would be easy if, for example, v were a strong function of y and u were 
a strong function of x. Then, each x, y could be chosen such that v(x,y) = 0 
and, as x is incremented, a solution would be crossed when u[x,y(x)] changes 
sign. Unfortunately, neither u nor v behaves this way. 

To develop functions that do behave this way, the following procedure was 
developed. First, for each x, choose y such that the absolute value of F is 
minimized. This can be done in several ways. The program uses a method that 
always guarantees finding a value if one exists. Essentially, the absolute 
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value of F squared and its gradients are computed. The value of y is altered 
in the direction indicated by the gradient until either the gradient changes 
sign or is so close to zero that convergence has been reached. Once the gradient 
changes sign, Muller's method (Ref. 25) is used to converge on the root. This 
is essentially a bisection method followed by inverse parabolic interpolation. 
For this searching process, it is not necessary to make the convergence criteria 
very tight, since only rough estimates are eventually needed in order to start 
the two-dimensional secant method described earlier. 


Now that a criterion for y has been established, it is only necessary to find a 
test function that will change sign when a solution is crossed while incrementing 
X. Such a function is given by 


This function acts as a very good test function because it represents the 
coordinate direction in the u,v plane along which the vector (u,v) changes 
most with X. When this coordinate changes sign as one goes from, say, x^ 
to X 2 with y^ and chosen so that the length of the vector (u,v) is minimized, 
then it is very likely that a solution has been crossed. Exceptions to this 
rule occur when one is near relative minima of the surface u^ + v that are 
not zero. To see this, consider the actual equations that are being solved. 

In order that the vector (u,v) is minimimi for each y, it is necessary that 
3(u^ + v^)/3y = 0. That is uu + vv^ = 0. Combining this with the above 
equation, we see we are finding an x and y such that the matrix equation 



( 212 ) 
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is satisfied. The matrix is Just the transpose of Jacobian of u and v with 
respect to x and y. 

This equation can be satisfied if either u and v are zero, or the Jacobian 
is singular. The Jacobian is necessarily singular at all relative minima of 
the surface u^ + except those at u = v = 0. In order to differentiate 
between those solutions to (212) that are due to singularities of the Jacobian 
and those that are due to u and v vanishing, we employ two different tests. 

First of all, when a singularity point is crossed, the determinant of the 
Jacobian should change sign. If this occurs, then the program rejects this 
point as a possible solution. Sometimes, however, the determinant does not 
change sign because either the convergence criterion used in the searching 
algorithm is too loose or because the singularity has a double root. In 
either case, the procedure is to test the condition number* of the transpose 
of the Jacobian matrix in the region near the suspected solution. If the 
condition number does not exceed a given input limit (e.g., around 80), 
then the point in question is usually a solution. 

Once it is determined that a potential solution has been crossed between x.| 
and x^, for example, the procedure is to (a) determine x^ based on the method of 
false position using the test function given in (211), (b) determine y^ to 
minimize |F|, and (c) use (x.j,y.j), and ss the required first 

three guesses for the two-dimensional secant method. 


The above procedure has been found to be most satisfactory for the conditions 
tested in this program. The search algorithm described above has several 
salient features. First, as mentioned earlier, the search method takes 
advantage of the fact that many of the computations are not a function of the 
imaginary part of namely y; This allows the minimization of |F| with respect 
to y to proceed with high efficiency. Secondly, and more importantly, 

*The condition number of a matrix. A, is a measure of how sensitive a solution 
to the system Ab=c is to perturbations in c. It is equal to the square root 
ratio of the absolute value of the largest eigenvalue of A'A to the smallest 
eigenvalue of A'A. For singular matrices, the condition number is infinite. 

For matrices that are nearly singular, the condition number will be quite large. 
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the procedure has been automated to the extent that the user only has to specify 
a frequency range and a maximum number of roots desired in that range. 

The algorithm will start at the lower end of the frequency range and will 
increment through it until either the maximum number of roots are found or the 
upper end of the frequency range is reached. This is a very powerful property 
since it does not require the user to have a clear knowledge of the location 
of any of the roots in the w plane. 
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PROGRAM USER'S GUIDE 


FSCSM MAIN PROGRAM 

The main program for the FSCSM computer model performs much of the 
input/output activity of the model and controls the sequencing through 
the various major subprogram blocks of the model. After writing the 
header page (Fig. 5), the main program reads in and writes out the data 
described in Table 2 of the input section of this report under the heading 
Main Control Section Input. The program then computes the area profiles 
through a call to subroutine AREA. It then begins its main do loop 
controlling the number of solutions to the nozzle admittance boundary 
equation that are desired. During the iterations for each solution to 
be found, the main program proceeds with successive calls to subroutines 
N0ZADM, HYDRDY, C0MBDY, and S0LVW in order to compute the downstream 
nozzle admittance factor, the feed system response parameters, the 
combustion dynamics coefficients, and performance calculations neces- 
sary to solve the nozzle admittance boundary equations respectively. 
Subroutine S0LVW also causes a call to subroutine CHAMDY which computes 
the oscillatory profiles. Further, during the first iteration in the 
do loop for the main program for the first solution, the main program 
also calls subroutine STEADY in order to obtain the steady-state profiles. 

The variable ISCNT.is the F0RTRAN variable set by the main program and 
altered in S0LVW which controls the type of iteration being performed. 

When ISCNT equals one or four, the search algorithm described on pages 54-56 
is called out. This is the initial condition at the beginning of each 
set of iterations to solve the nozzle admittance boundary equation. When 
ISCNT equals five, the two-dimensional secant method is being performed 
in S0LVW. 




I’AGB blank 
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Figure 5. Logic Diagram for FSCSM Main Program 
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The variable KWHERE, which is set by subroutine S0LVW, controls the logical 
flow in the main program subsequent to a call to S0LVW. When ISCNT equals 
one or four and KWHERE equals one, the control is returned to that portion 
of the main program which starts the calculation for the next value of u 
being tried by the search algorithm. When KWHERE equals two and ISCNT 
equals one or four, control is passed to that portion which will perform 
calculations for a perturbed w in order to compute the derivative of 
^Nu ■ ^Nd respect to When ISCNT equals 5 and convergence has 
not been reached, the normal exit from S0LVW also sets KWHERE equal to 
two. For this case however, no derivatives are calculated. The program 
will Just proceed with successive calls to N0ZADM, HYDRDY, and C0MBDY 
in order to compute the downstream nozzle admittance factor, the feed 
system response, and the combustion dynamic coefficients, respectively. 

It then proceeds to S0LVW in order to obtain an updated estimate of w 
using the two-dimensional secant method. 

When KWHERE equals three or four, convergence on a solution to the nozzle 
admittance boundary equation has been attained. For this case, control 
is passed to that portion of the main program which prints the final 
results for that solution. 

When KWHERE equals five, control is passed to the end of the main do loop 
in the main program. The output portion is bypassed. This only occurs 
when an error was detected by subroutine S0LVW. 


R-9808/61 


S0LVW 

This subroutine performs many of the calculations and controls most of the 
logical flow required to match the downstream boundary condition on the noz- 
zle admittance. The FORTRAN variable ISCNT controls the logical flow with- 
in subroutine S0LVW (Fig. 6). If ISCNT equals one, then that portion of 
the subroutine used for searching the w-plane for possible solutions to 
the nozzle admittance boundary equation is used. Two calls to S0LVW are 
used for this purpose. During the first call, the imaginary part of w 
is determined so that the absolute value of the upstream nozzle admittance 
minus the downstream nozzle admittance is minimized. The second call is 
made in order to complete the computation of the Jacobian of this dif- 
ference with respect to the real portion of w. (The derivatives of the 
difference with respect to the imaginary part of w are comouted during 
the first call.) When ISCNT equals one, tests are also made to determine 
if a solution is nearby. The actual test function and the logic employed 
is described on pages 54-56. 

Once it is determined that a possible solution is bracketed by two suc- 
cessive frequencies, the variable ISCNT is set equal to four. Subroutine 
S0LVW performs the same calculations for this value of the variable ISCNT 
as it does when ISCNT equals one. The only difference occurs at the end 
of the second call to S0LVW. At that point, checks are made to ensure 
that the potential solution is in fact an actual solution and not due to a 
singularity in the Jacobian. If the error passes certain criteria and at 
least three iterations have been performed with ISCNT equal to four, then 
ISCNT is set equal to five for subsequent calls to S0LVW. Between each 
iteration for ISCNT equal to four, as well as for the first iteration for 
ISCNT equal to four, the real part of the frequency is modified using the 
method of false position or the bisection method, depending upon the value 
of the iteration counter, KSCNT4. 

When ISCNT equals five, no derivatives are computed. For this situation, 
subroutine S0LVW only computes the difference between the upstream and 
downstream nozzle admittances based on the current value of w. It then 
checks to determine whether convergence has been obtained. If not, the 
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Figure 6. Logic Diagram for Subroutine S0LVW 
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Figure 6. (Concluded) 
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value of 0 ) is updated according to the two-dimensional secant method 
described on page 53. 

The variable which controls the flow in the main program subsequent to a 
call to S0LVW is KWHERE and is described in the section of this report 
dealing with the main program. 

CHAMDY 

This subprogram is called by XIMA6F in order to compute the oscillatory 
pressure, temperature, velocity, mixture ratio, and density profiles. 

From these quantities (Fig. 7), it determines the upstream nozzle admit- 
tance. It solves for the oscillatory profiles by solving the linearized 
set of differential equations presented on pages 41 and 42. This is done 
using a second order implicit finite difference scheme. Those integrals 
appearing in the vaporization expression which cannot be integrated analyti- 
cally are numerically integrated using the trapezoidal rule. 

Because the differential equations represent a linear initial value problem, ' 
the finite difference equations are also linear and one can "march off" the 
solution from the initial plane. Since the four differential equations 
are coupled, replacing them at each axial position by their finite difference 
approximation results in a four by four system of complex linear equations. 
Because of the nature of the differential equations, the resulting matrix 
equations are essentially diagonally dominant and can therefore be solved 
very quickly using Gaussian elimination with the diagonal element used for 
pivoting. 

XIMAGF 

This subroutine is called by S0LVW and ZER0. Its main function is to compute 
the difference between the upstream and downstream nozzle admittances (Fig. 

8). When the program is still performing the search algorithm, this routine 
computes the derivative of this difference with respect to the imaginary 
part of u) as well as the derivative of the absolute value of this difference 
squarednvi th respect to the imaginary part of m. — 
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Figure 7. Logic Diagram for Subroutine CHAHDY 
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Figure 8. Logic Diagram for Subroutine XIMAGF 
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CHMC0N 

This routine is called by CHAMDY and calculates certain parameters 
necessary for determining the coefficients used in CHAMDY. 

C0MMAT 

This routine solves the four by four complex system of linear equations 
discussed in the section describing subroutine CHAMDY. It uses Gaussian 
elimination with pivoting on the diagonal. 

ZER0 

This routine is called by subroutine CHAMDY. Its function is to find the 
zero of a given functional when that zero is bracketed both above and below. 

The functional in this case is the derivative with respect to the imaginary 
part of 01 of the absolute value squared of the difference between the 
upstream and downstream nozzle admittances. Finding the zero of this func- 
tional is done in order to minimize the error in the difference between the 
nozzle admittances with respect to the imaginary part of u . The method used 
by subroutine ZERO is due to Muller (Ref. 25). It essentially involves a 
bisection step followed by inverse parabolic interpolation to determine the 
next guess. 

STEADY 

This routine (Fig. 9) is called by the main program to determine the time inde- 
pendent solution to the set of differential equations given on page 40. These 
equations have been analytically integrated on pages 42-44. This subroutine uses 
these latter equations to determine the steady state profiles. Also, several 
parameters which are a function of these steady state variables are computed 
and saved for subsequent use by the chamber dynamics subprogram, CHAMDY. If 
the FORTRAN variable IPRSTE is greater than zero, a printout of the steady 
state profiles will be given. 

C0MBDY 

This subprogram (Fig. 10) calculates the fuel and oxidizer o^mbustion coupling 
coefficients required for the determination of the time oscillatory vaporization 
rates needed to solve the chamber dynamics. The equations for these parameters 
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Figure 9. 


Logic Diagram for Subroutine STEADY 
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Figure 10. Logic Diagram for Subprogram C0MBDY 
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are given on pages 33-34. During the first pass into this program, the com- 
bustion dynamic input variables are read in from input device 5 and written 
out onto device 6. A discussion of these variables is given in the Model 
Input Section. The logical flow in C0MBUY is controlled by the F0RTRAN 
variable INPC0M. If this variable is less than or equal to one, the combus- 
tion coefficients are computed for the current frequency only. If INPC0M 
is greater than one, these coefficients are computed for the entire fre- 
quency table, FREQT (e.g., from 10 to 1000 Hz), and saved on tape/disk 
ITAPE for subsequent use by the main program. 

AREA 

This subroutine is called by the main program. It computes the area pro- 
files and axial distance profiles necessary for solution of the steady-state 
and transient profiles. 

L0CFAC 

This routine is used to determine the subscript, I, within an ordered array, 
TX, such that the input argument, X, is in the interval TX(I), TX(I+1). 

This routine also returns the interpolating factor FX = (X-TX( I)/(TX( I+l )- 
TX(I)) for subsequent use in linear interpolation. 

HEAD 

This subroutine is called by the main program to print the heading page which 
gives the title of the program, by whom and where it was developed, and the 
program sponsor. 

HYDRDY 

Subroutine HYDRDY (Fig. 11) is called by the main program to calculate the 
frequency domain characteristics of the feed system. Functions performed 
by HYDRDY are (1) reading of input data describing the physical attributes 
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of the feed system components, (2) generation of a matrix of linear dif- 
ferential equations representing the complete feed system, (3) solution 
of the feed system equations to yield the amplitude and phase response of 
all feed system pressures and flowrates as a function of chamber pressure 
oscillation amplitude and frequency, and (4) generation of tabulated output 
of injector flowrate frequency response for use by the main program. 

A basic assumption of subroutine HYDRDY is that the feed system being 
modeled can be represented by the generalized schematic of Fig. 4 
(or by some portion of this schematic). This assumption is necessary 
because HYDRDY sets up and solves the complete set of simultaneous 
equations representing the Fig. 4 schematic. By assigning very high 
resistance and very short length attributes to any of the 30 numbered 
line segments of the generalized schematic, whose segment can effectively 
be excluded from having any effect on the frequency response character- 
istics of the rest of the system. With this approach a wide variety of 
feed systems can be modeled with no changes to the program other than 
the input data. 

Figure 11 shows the functional block diagram for subroutine HYDRDY. When 
called, the subroutine initially zeroes the values of all of the elements 
of the coefficient and time delay matrices C and TD in labeled coninon 
block F. The values of various fixed input arguments required by the fre- 
quency response subroutine (FRESP) are then initialized. 

Input argument IR directs the reading of subroutine HYDRDY input data. If 
IR is zero or less, the program assumes that all required data has pre- 
viously been read and the data read function is bypassed. If IR = 1, the 
program assumes that no hydrodynamic data has been read and proceeds to ini- 
tialize all hydrodynamic input variables to values which will exclude all 30 
line segments and both injectors of the generalized Fig. 4 feed system sche- 
matic. Control is then passed to statement 45 for reading of input data. If 


R- 9808/73 



IR = 2, the assumption is made that most input data is already set up 
(such as from a previous case during the same program run). Control 
is passed directly to statement 45 for reading any changes to the 
input data. 

Input data reading for HYDRDY is in the NAMELIST format (NAMELIST 
name HYD) and is normally in the form of card input on logical unit 
number 5. However, if the program is run in a timesharing environment, 
an option is provided for reading data from a timesharing terminal. 

This option is controlled by variables IRFLAG and ITERM. Both of these 
variables are stored in labeled common block/F/ and can be changed by 
input to NAMELIST HYD. Variable IRFLAG is tested in statement 45 and 
if non-zero specifies reading NAMELIST HYD data from unit 5. If time- 
sharing terminal input is desired, variable ITERM is set (by input 
data read or block data initialization of labeled common block/F/) to 
the logical unit number to the timesharing terminal. If variable ITERM 
is non-zero, statement 46 sets IRFLAG to 1. Thus, once terminal data 
input has been specified, all subseouent data reads will default to the 
terminal. Card input can be respecified (for a subsequent data case) by 
entering IRFLAG = 0 in the terminal data input. 

Input variables for subroutine HYDRDY are described in detail in the Hydro- 
dynamic Input Section. The variables include the length (L), cross-sec- 
tional area (A), propellant sonic velocity (V), propellant density (RHOL), 
hydraulic resistance (R), and wall compliance (CW) for each of the 30 
numbered waterhammer segments in the generalized Fig. 4 feed system 
schematic. Input variables for the left ("tf") injector of Fig. 4 are 
resistance (R0) ,inertance (Z^), volume (V0L(», propellant sonic velocity 
(V0), and injector deflection constant (K0). The corresponding input 
variables for the right ("F") injector are RF, ZF, V0LF, VF and KF. 

The designation of the two injectors as "0" and "F" is a notational 
convenience for cases in which the feed system being modeled has only 
one injector and sufficiently simple flow paths so that both oxidizer 
and fuel ^sterns can simultaneously be laid out on the Fig^4 schematic. 
Such cases have the advantage of reduced computer time because the 
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frequency response of both fuel and oxidizer feed systems is obtained 
with a single call to subroutine HYDRDY. Of the fuel and oxidizer 
feed systems overlap when laid out on the Fig. 4 schematic, subroutine 
HYDRDY must be called twice - once for each feed system. 

When data input is complete, a value of 2 (or greater) for variable 
IWRITE specifies a printout of all input data on logical unit 6. 

IWRFTE = 0 is the default and specifies no printout of input data. 

Next, control is passed to the D^ loop at statement 100 in which the 
input values for propellant density, propellant acoustic velocity 
and segment wall compliance for each of the 30 waterhammer segments 
of the Fig. 4 schematic are combined to yield an effective acoustic 
velocity for each segment. The subsequent statements, up to statement 
400, combine the input variables as required to yield the constant 
coefficients of the 57 linear waterhammer and injector equations des- 
cribing the complete Fig. 4 feed system. Simultaneous solution of these 
57 equations, at each specified input frequency, yields the oscillatory 
amplitude and phase response of all pressures and flowrates in the feed 
system to inputs via chamber pressure oscillations at that freouency. 

At statement 500 a call to subroutine FRESP is made to obtain the frequency 
response solution of the feed system equations. Initially, however, at 
statement 400 the value of input argument INPHYD is tested to determine 
the desired output from FRESP. If INPHYD is greater than 1, HYDRDY will 
call FRESP to calculate the feed system frequency response for each of 
the frequencies in array FREQT. The total number of frequencies is given 
by variable NFREQT and may range from 1 to 100. Both the variable, 

NFREQT, and the array, FREQT are in labeled common block/COMTAP/. If the 
value of INPHYD is less than or equal to 1 HYDRDY will call FRESP to 
calculate the feed system frequency response for the single frequency 
given by input argument, FRE. 

a 

Output^ata from subroutine HYDRDY consists of a pair oT complex numbers 
for each specified input frequency. If INPHYD was specified as £ 1 then 
the output numbers GINJ^iX and GINJFU are returned in labeled common 
block/F/ and also in the HYDRDY argument list as GIN() and GINF. The real 
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and imaginary parts of complex number GINJ0X (GIN^) represent the 
amplitude and phase angle respectively of frequency 

FREQ (FRE). Similarly, the complex number GINJFU (GINF) re- 
presents the amplitude and phase angle of /p^ at frequency 

FREQ (FRE), W0IN, WFIN and PCIN are the input normal values for 
the oxidizer injector flowrate, fuel injector flowrate and chamber 
pressure, respectively, from the HYDRDY argument list. 

If INPHYD was specified as >1, then rather than a single pair of complex 
numbers representing oscillatory injection flowrates, HYDRDY returns two 
arrays, 6INJ0T and GINJFT, containing the oscillatory injection flowrate 
amplitude and phase data for each of the NFREQT frequencies in array 
FREQT. The output arrays GINJOT and GINJFT are stored in labeled coimon 
block /C0MTAP/ and are also written out on the output device whose logical 
unit number is designated by variable ITAPH in labeled common block /C0MTAP/. 

The order of storage on the output device is GINJ(>T(I), GINJFT(I), for 
I values from 1 through NFREQT. After writing the GINJ(>T and GINJFT 
arrays on the output device HYDRDY sets the value of variable INPHYD to 
3. Also, before returning control to the main program, HYDRDY tests the 
value of variable IWRITE. If IWRITE is non-zero, each specified frequency 
and the corresponding value of GINJ0T and GINJFT are written out on logical 
unit 6. If only one frequency was specified (INPHYD 1), then only the 
single point values of FREQ, GINJ0X and GINJFU are written out. 

It should be noted that although output from a single call to HYDRDY con- 
tains values for both "oxidizer" and "fuel" oscillatory injection flowrates 
(at one or more frequencies), the output values actually refer to the "0" 
and "F" injectors of the Fig. 4 schematic. Thus, unless both oxidizer and 
fuel feed systems can simultaneously be modeled with the Fig. 4 layout, 
it is necessary to call HYDRDY twice - once for the oxidizer feed system 
and once for the fuel feed system. 

FRESP 

Subroutine, FRESP (Fig. 12) is used to obtain the frequenc y, d omain character- 
istics of the feed system indirectly from input data that describes the physical 
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Figure 12. Logic Diagram for Subroutine FRESP 
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characteristics of the feed system. (The actual input to FRESP is 
generated by the subroutine HYDRDY which orders the physical character- 
istics of the system into specific matrices of coefficients that FRESP 
can use as input data.) 

FRESP merely solves for the variables X. in the following relation- 
ship: 

C {X} = a Y 


where Y is a single input variable that represents a unit value of 
the injector end combustion chamber pressures. The matrix a then 
relates the specific pressure input to each applicable equation that 
contains combustion chamber pressure ( a may contain both static and 
dynamic terms.) The matrix C is simply the coefficients of the linear 
differential equations that represent the physical system. The values 
of the coefficients for the a and C matrices are computed by the 
subroutine HYDRDY. 


The FRESP matrices can be expressed as: 


C 


ijk 



(X.} 



• Y 


with the differential operator defined as S = Ju, where J = yTT~ and 
oj is the frequency. The matrices may be broken down to provide real 
matrices and imaginary matrices. 


C 


ijl 



+ C, 


ij5 


4 

0 ) 


+ J 







+ J 2 03 


^•4 



Y (213) 


Since the time delay coefficients used in the differential equations are 
of the form e"'^^-X, which is equivalent to e""^^“-X, and since e"^^ = 
cos(y) + j sin(y), these terms may be added to the previous Informed 
real and imaginary matrices to give: 
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■ ^1j3 ^ +...•+ cos to) 

f 3 

+ IJ ^ " ^ij4 ^ ^ •••• SI 


n(,,j «.)] [Xj| 


3«^ ^i3 ^ ^ •••• cos 


+ J I ^ ■*" 


and solved for 




[*<)= 


Ijl - S-j3 " •••• I 


J I ^*|j 2 ^ " ^'^j 4 ^ •••♦ sin w) 


dii - 3i3 (I) + + COS (t^. Hi) I + 


J I w ~ a^.^ 0 ) +....+ sin (to)) 


( 214 ) 


(215) 


The ipatrices are multiplied and then solved for j j in the subroutine 
C0GAEL which employs the standard Gaussian elimination procedure for 
solving linear equations. The | j solution is still separated into 
real and imaginary components, and are simply combined to form a vector 
for each variable. The procedure is repeated for each frequency being 
considered. 
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C06AEL 

This subroutine is called by the hydrodynamic frequency response sub- 
routine, FRESP, to triangularize the complex matrix of feed system 
equations. Back substitution into the triangular system of equations 
is subsequently performed by subroutine FRESP to yield the real and 
imaginary portions of each feed system variable. 

The conventional method of Gaussian elimination is employed by C0GAEL 
to triangularize the system of equations. The reduction process pro- 
ceeds in column order from left to right. First the complex element 
with the largest absolute value in the current ("oivot") column at 
or below the diagonal is located. Then the rows are interchanged if 
required to move this maximum element (oivot element) to the diagonal. 

The row interchange serves to minimize the round-off errors from the 
subsequent reduction process. The pivot row (row containing the oivot 
element) is then divided by the pivot element yielding 1.0 from the 
pivot element. Finally, the elements in the pivot column below the diagon- 
al are eliminated by subtracting the appropriate multiple of the pivot 
row from each row below it. The subtraction is not actually performed 
on the elements below the diagonal since these elements do not enter 
into the subsequent back substitution process performed by subroutine 
FRESP. 

It should be noted that the above discussion refers to the complex 
matrix as if the elements were single numbers. The actual elements 
are stored as two numbers in each row, the real portion to the left 
of the constant term, and the imaginary portion on the right. This 
distinction does not alter the elimination process except that two 
separate numbers must be operated on at each step. 

TDPL0T 

If the value of I CRT is greater than 0 this subroutine is called by the 
frequengy response routine FRESP to generate CRT plots Qi_the gain and 
phase of the output variables as a function of frequency. 
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The input arguments to TDPL0T are, W, an array of up to 101 freouencies; 
Y, an array of gain or phase values; NFP, the number of data points in 
arrays W and Y; TL, the lowest desired frequency grid line; XR, the 
highest desired frequency grid line: and LL, a flag indicating a gain 
plot if 1 or a phase plot if 2. 

Initially TDPL0T scans the Y array for the maximum and minimum values 
and generates values for the Y axis grid scaling. The first value of 
Y is not included in this scan. This allows an initial very low value 
of frequency to be used to approximate the system DC frequency response 
without upsetting the plot frequency scaling. TDPL0T uses the standard 
graphics package routines for the SC-4020 to generate the plot orids 
and plot the data points. If the value of LL is 1 the CRT frame is 
advanced and a plot of Y(I) versus W(I) is made on the bottom half of 
the page. If LL is 2, the frame is not advanced and the plot is made 
on the top half of the page. 

In addition to the plots, TDPL0T prints the numeric value of the first Y 
array element immediately above the corresponding plot. This element 
typically corresponds to a very low frequency value (default value of 
.001 cps in subroutine FRESP) which is well below the frequency range 
desired for the plot and approximates the DC value of the output 
variable. 

TDPL0T does not generate any titles or identifi eating information on the 
plots. 
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N^ZADM, (RKTDIF, RKTZ, RKZDIF, TADAMS, ZADAMS) 

This routine and its support routines, is called by the main program to deter- 
mine the nozzle admittance based on downstream conditions. The programs were 
developed and programmed by Georgia Institute of Technology and the user is 
referred to Ref. 23 for a complete description of these routines. The main 
nozzle admittance program was modified by Rocketdyne so that input data could 

be read if required and also the nozzle admittance saved on a tape unit ITAPN 
(Fig. 13). 


R-9808/82 



SUBPROGRAM N0ZADM 
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Figure 13. Logic Diagram for Subroutine N0ZADM 










FSCSM COMPUTER MODEL INPUT 


This section, and the Hydrodynamic Input section, describes the input necessary 
to run the FSCSM computer model. The input is broken up into four major control 
sections. These are the main control section input, the nozzle admittance con- 
trol section input, the hydrodynamics control section input, and the combustion 
dynamics control section input. Table 2 lists all the variables that are to be 
input for each control section. This input is in the usual 80 character card form. 
Listed in Table 2 for each control section are each card's number and format, the 
variable names appearing on each card, and a brief description of each variable 
appearing in the list. 

The main control section requires either eight or nine cards depending upon the 
input value of INPN0Z. (If INPN0Z is less than or equal to three, the last card 
in this section should be input). The first two cards, cards 1 and 2, should 
contain title and case identification information. These are read in alpha num- 
eric format and printed at the top of almost every page of output to identify the 
case being considered. 

The third card contains control information for various run options, file numbers 
for the auxiliary storage of datasets used by the program, print codes, and the 
number of mesh points to be used in the analysis. The control flags are INPHYD, 
INPC0M, and INPN0Z. These flags allow the user to save the results from the cur- 
rent case or use the results from a previous case for the current case. It is 
recommended that these datasets be set up as permanent files in order to use them 
for subsequent job submittals. 

The first of these flags, INPHYD, controls the hydrodynamics feed system dataset. 

If it is input as one, no action will be taken with respect to saving or reading 
information on or from a dataset. The hydrodynamics coupling terms will be re- 
computed each time the frequency changes. If INPHYD is input as two, a table of 
hydrodynamic coupling terms will be generated for the frequency range specified 
by the input variables NFREQ, FREQMI, and FREQMX. This table will be saved on 
file ITAPH (also input on this card) and used to linearly interpolate in each 
time the frequency changes. If INPHYD is input as three, the program assumes a 
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TABLE 2. FSCSM INPUT DATA 


CONTROL 

SECTION 

Main Program 
Input 


CARD NO./ VARIABLE 

FORMAT NAME DESCRIPTION 


THESE CARDS MUST ALWAYS BE INPUT 


1,2 ( (TITLE(I ,J) , Title information is input on the 

(18A4) 1=1,18), J=l,2) first two cards in the input data- 

set. These should be used to 
identify the case being run. 

3 

(1216) INPHYD Code used for hydrodynamic calcula- 

ti ons . 

1; Hydrodynamic coupling terms (HCT) 
will be computed each time the 
frequency changes. 

2; A table of HCT's will be computed 
using frequency table, FREQT, 
and saved on file ITAPH. Each 
time the frequency changes in 
trying to satisfy the nozzle 
admittance boundary condition, 
the HCT will be interpolated for 
in that table. 

3: Table of HCT's already resides on 

file ITAPH from an earlier run. 

When the frequency changes , the 
HCT will be interpolated for in 
that table. 

INPC0M Code used for the combustion dynamics 

coefficients (CDC). It can take on 
values from one to three. It has the 
same meaning with respect to the CDC 
as INPHYD does with respect to the HCT. 
The tables are saved on file ITAPC. 

INPN0Z Code used for the downstream nozzle 

admittance term (DNAT). 

1: The DNAT is computed each time 

the frequency changes. 

2: A table of DNT's versus fre- 

quency using frequency table, 

FREQT, is computed and saved 
internally. Each time the fre- 
quency 'changes in trying to 
satisfy the nozzle admittance 
boundary condition, the DNT will 
be interpolated for in that table. 
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TABLE 2. (Continued) 


CONTROL CARD NO./ VARIABLE 

SECTION FORMAT NAME 


Main Program 
Input (Cont.) 


ITAPH 

ITAPC 

ITAPN 

IPRHYD 

IPRC0M 


IPRN0Z 


IPRCHM 


DESCRIPTION 


3: Same as 2, but the table is 

also saved on file ITAPN, 

4: A table of DNT's versus frequency 

already resides on file ITAPN 
(from a previous run). Inter- 
polation is the same as 2 and 3. 

File used to save hydrodynamic coup- 
ling term table. Must be specified 
if INPHYD^2. 

File used to save combustion coeffi- 
cients table. Must be specified if 
INPC0M^3. 

File used to save nozzle admittance 
terms table. Must be specified if 
INPN0Z>^3. 

Code used for hydrodynamics coupling 
term table printout. 

0: Table will not be printed. 

1: Table will be printed. 

Code used for combustion coupling 
coefficients printout 

0: Coefficients will not be 

printed. 

1: Coefficients will be printed 

at each freouency found to 
satisfy the nozzle admittance 
boundary condition. 

Code used for the downstream nozzle 
admittance terms table printout. 

0: Table will not be printed. 

1: Table will be printed. 

Code used for the oscillatory profiles 
printout. 

0: Profiles will not be printed. 

1: Profiles will be printed at 

each .freouency found to satisf, 
the nozzle admittance boundary 
condition. 
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TABLE 2 . (Continued) 


CONTROL CARD NO./ VARIABLE 

SECTION FORMAT NAME 

Main Program IPRSTE 

Input (Cont.) 

NXP 

4 

(6E12.8) 

XN0Z 

RINJ 


DESCRIPTION 

Code used for the steady-state profiles 
printout. 

0: Profiles will not be printed. 

1: Profiles will be printed. 

Number of axial positions to be used 
between the XO plane (start of vapori- 
zation) and the nozzle inlet plane. 

(Both the XO and nozzle inlet planes 
must be counted.) 

Axial coordinate of start of vaporiza- 
tion plane, inches. 

Nozzle inlet plane, inches 

Radius of combustion chamber at injector 
face, inches. 


5 

(616) 


GAM0 

C0 

deLp 

NR00T 

IWRT 


Ratio of specific heats of combustion 
gas (Cp/Cy) evaluated at overall mixture 
ratio, umtless. 

Sonic velocity evaluated at overall 
mixture ratio, ft/sec. 

Oscillatory non-dimensional pressure 
amplitude at injector face, unitless. 

Absolute value of this variable is the 
number of frequency solutions that will 
be searched for starting at the frequency 
specified by the real part of 0MEGA and 
ending at FROMAX. This variable may be 
input as either positive or negative. 

(See input of0MEGARfor explanation.) 

Intermediate output dump code used to 
write the oscillatory profiles solved 
for in CHAMDY for each iteration. 

0: Oscillatory profiles will not 

be printed out between iterations. 
1: Oscill^ory profiles will be 

printed for each iteration. 
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TABLE 2. (Continued) 


CONTROL CARD NO,/ VARIABLE 

SECTION FORMAT NAME 


Main Program IWSKP 

Input (Cont.) 


KNTMX 


KNTRMX 


KNTSMX 


6 0MEGAR 

(6E12.8) 

0MEGAI 

FRQMAX 

DELFRQ 


DESCRIPTION 

Intermediate output dump code used in 
subroutine S0LVW. 

0: No intermediate output will 

be printed from S0LVW during 
iterations. 

1: Limited intermediate output 

will be printed by S0LVW 
2: Extended intermediate output 

win be printed by S0LVW. 

Maximum number of iterations allowed to 
minimize the error in the nozzle admit- 
tance boundary eouation with respect to 
the imaginary part of u. 

Maximum number of times the frequency will 
be allowed to be changed by DELFRO during 
the searching algorithm between each 
solution. 

Maximum number of iterations allowed for 
the convergence of the two-dimensional 
secant method used in S0LVW. 

Starting value for the real part of complex 
frequency of NR00T>0. This should be in- 
put in units of Hertz times Zir. If NR00T<0, 
this should be input in Hertz. 

Starting value for the imaginary part of 
complex frequency. This should be input as 
the growth coefficient if NR00T>0. It 
should be input as the decrement if 
NR00T<0. 

Maximum freouency above which no solutions 
to the nozzle admittance boundary equation 
will be sought. Hertz. 

Increment used to adjust the frequency 
during the searching portion of the 
algorithm to solve the nozzle admit- 
tance boundary equation. Hertz. 
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TABLE 2. (Continued) 


CONTROL 

SECTION 

Main Program 
Input (Cont.) 


CARD NO./ VARIABLE 

FORMAT NAME 


DELMX 


CTEST 


7 EPSF 

(4E12.8) 


EPSX 


DESCRIPTION 


Maximum allowable change in the growth 
coefficient between two successive 
iterations in the portion of the pro- 
gram that minimizes the error in the 
nozzle boundary condition eouation with 
respect to the growth coefficient, sec"'. 

Upper bound on the condition number of 
the transpose of the Jacobian of the dif- 
ference between the upstream and down- 
stream nozzle admittances with respect 
to the complex frequency. If the condi- 
tion number of that matrix exceeds CTEST 
for a given frequency, then it is assumed 
the Jacobian is singular near that fre- 
quency and hence a solution will not be 
sought at that point. Unitless. 

Relative error criterion used during the 
search algorithm/or the portion of the 
program that minimizes the eror, HN, in 
the nozzle admittance boundary equation with 
respect to the growth coefficient, (dj. To 
obtain convergence, it is necessary that 



Uni tl ess. 

Relative error criterion used during the 
search algorithm for the portion of the 
program that minimizes the error, HN, in 
the nozzle admittance boundary equation with 
respect to the growth coefficient, wj. 

To obtain convergence, it is necessary 
that 

|(i)f - «o. I / |to, |<EPSX 

M 4 / ^2 

where the subscripts 1 and 2 refer to two 
successive iterations. Unitless. 
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TABLE 2. (Continued) 


CONTROL CARD NO./ VARIABLE 

SECTION FORMAT NAME 


Main Program EPSFS 

Input (Cont.) 


EPSXS 


8 

(3E12.8) 


PC 

MB0XI 


MBFUI 


DESCRIPTION 


Tightened relative error criterion used 
to determine if convergence has been 
obtained while iterating to solve the 
nozzle admittance boundary equation. 

To obtain convergence, it is necessary 
that 




Unitless. 


A. < EPSFS 


Tightened relative error criterion used 
to determine if convergence has been 
obtained while iterating to solve the 
nozzle admittance boundary equation. 

To obtain convergence, it is necessary 
that 



where the subscripts 1 and 2 refer to 
successive iterations. Unitless. 

Steady-state chamber pressure, psia. 

Oxidizer injection flowrate, Ibm/sec 
Fuel injection flowrate, Ibm/sec 


INPUT THIS CARD ONLY IFINPN0Z<3 


9 

(I12,2E12.8) 


NFREQT 

Number of points in frequency table. 

FREQMI 

Minimum frequency in frequency table. 


Hertz. 

FREQMX 

Maximum frequency in frequency table. 
Hertz. 
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TABLE 2. (Continued) 


CONTROL CARD NO./ VARIABLE 

SECTION FORMAT NAME DESCRIPTION 


Nozzle 

Admittance 

Program 


INPUT THE CARD IN THIS CONTROL SECTION 
ONLY IF INPN0Z<3 


1 RCCX 

(4E12.8) 

RCTX 

ANGLEX 

CRR 


Ratio of the radius of curvature at the 
nozzle inlet to the chamber radius at 
nozzle inlet, unitless (see Fig. 15 ). 

Ratio of the radius of curvature upstream 
of the throat to the chamber radius at 
nozzle inlet, unitless (see Fig. 15). 

Nozzle convergence half angle, degrees, 
(see Fig. 15 ). 

Contraction ratio, cross-sectional area 
of chamber/throat area, unitless 
(see Fig. 15 ) 


Hydrodynamics 

Program 


THIS CONTROL SECTION READS IN ITS INPUT DATA 
IN NAMELIST FORMAT. THE NAMELIST NAME IS 
/HYD/. INPUT THIS DATA ONLY IF INPHYD<2. 


&HYD* 

Input these characters starting in column 

two of the first card of the 

input. 

NAMELIST 
Variables 
in any 
order 

See Table 3 for a listing of the NAMELIST 
data input names. The accompanying text 
describes the meaning of these variables. 

&END* 

Character string denoting the end of the 
NAMELIST input block. 

*For Uni vac 1110 systems, use $HYD 

and $END 
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TABLE 2. (Continued) 


CONTROL CARD NO./ 

SECTION FORMAT 

Combustion 

Dynamics 

Program 

(6E12.8) 


VARIABLE 

■ name description 


INPUT THE CARDS IN THIS CONTROL SECTION 
ONLY IF INPC0M<2 


XK0X 


Klystron constant for oxidizer jet, 
inches . 


TAUB0X 

VB0X 

DELH0X 

TDRAG0 

ADV0X 

2 ADD0X 

(5E12.8) 


Steady-state oxidizer vaporization 
time delay, sec. 

Steady-state oxidizer injection 
velocity, ft/sec. 

Pseudo energy term for oxidizer, 
Btu/lbm. 

Steady-state oxidizer drag time 
delay, sec. 

Velocity exponent for the oxidizer 
atomization process, unitless. 

Oxidizer liquid jet diameter exponent, 
unitless. 


DELV0X Steady-state velocity difference between 

oxidizer droplets and gas stream normal- 
ized to the sonic velocity at the overall 
mixture ratio, unitless. 

Steady-state oxidizer Nusselt number used 
in vaporization expression, unitless. 

DT0XDM Partial derivative of oxidizer vaporization 

time delay with respect to mixture ratio, 
holding the vaporization blockage term, 
drop diameter, and Nusselt number constant, 
sec. 


XIMP0X 

3 XKFU 

(6E12.8) 

_ TAUBFU 

VBFU 


Oxidizer jet impingement point, inches. 

Klystron constant for fuel jet, inches. 

Steady-state vapori^ion time delay, sec 

Steady-state fuel injection velocity, 
ft/sec. 
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TABLE 2. (Concluded) 


CONTROL 

SECTION 

CARD NO./ 
FORMAT 

VARIABLE 

NAME 

DESCRIPTION 

Combustion 
Dynami cs 
Program 
(Cont. ) 


DELHFU 

Pseudo energy term for fuel, Btu/lbm. 


TDRAGF 

Steady-state fuel drag time delay, sec. 



ADVFU 

Velocity exponent for fuel atomization 
process, unitless. 


4 

(5E12.8) 

ADDFU 

Fuel liquid jet diameter exponent, 
unitless. 



DELVFU 

Steady-state velocity difference between 
fuel droplets and gas stream normalized 
to the sonic velocity at the overall 
mixture ratio, unitless. 



NUBFU 

Average steady-state fuel Nusselt 
number used in vaporization expression, 
unitless. 



DTFUDM 

Partial derivative of fuel vaporization 
time delay with respect to mixture ratio 
holding the vaporization blockage term, 
droplet diameter, and Nusselt number 
constant, sec. 



XIMPFU 

Fuel jet injection impingement point, 
inches. 


5 

(5E12.8) 

MWG 

Steady-state molecular weight of the 
gas at the overall mixture ratio, 
Ibm/lbm-mole. 



CS 

Characteristic velocity at the overall 
mixture ratio, ft/sec. 



DRGDMR 

Partial derivative of gas constant with 
respect to mixture ratio evaluated at 
the overall mixture ratio, ft-lb/lb/°R. 

— 


DCSDMR 

Partial derivative of characteristic 
velocity with respect to mixture ratio 
evaluated at th^ overall mixture ratio, 
ft/sec. 



DHDMR 

Partial derivative of gas reference 
enthalpy with respect to mixture ratio 
averaged over the mixture range during 
steady-state operation, Btu/lbm. 
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table of hydrodynamic coupling terms versus frequency already resides on file 
ITAPH in the format used during generation of such a table when INPHYD is input 
as two. The program will interpolate in this table in order to obtain the hydro- 
dynamic coupling terms each time the frequency changes. The use of interpolation, 

once a table has been generated, substantially reduces the computer run time for * 
each case run. 

The other two control flags input on the third card control the datasets for the 
combustion dynamics coefficients and the nozzle admittance factors. These flags 
are similar to INPHYD. Their description is given in Table 2 . 

Also input on the third card of the main control section input are the file num- 
bers of the datasets discussed above and print control flags for the various forms 
of output one can obtain. These are all self-explanatory. The user need only 
refer to Table 2 in order to determine the values that should be input for the case 
being considered. 

The final entry on the third card is for the variable NXP, the number of points 
to be used for the axial distance and area arrays. This controls the step size 
that will be taken during the integration of the chamber dynamics equations; i.e., 
step size = (XN0Z - X0)/(NXP-1) where XO is the axial location of the start of 
vaporization plane and XN0Z is the axial location of the nozzle inlet plane. The 

values of XO and XN0Z are both input on the very next card read in (card 4 in 
Table 2). 


The start of vaporization plane (XO) is calculated by plotting the percent un- 
burned of both fuel and oxidizer that is calculated by the DER program (or equiv 
alent steady-state combustion model) as a function of distance from the injector 

face (Fig. 14). These plots are then extrapolated back to 100% unburned and the 
axial location of this point is XO. 


Also input on card 4 are the radius of the combustion chamber at the injector 
face, RINJ, the ratio of specific heats (Cp/C^), GAM0, the sonic velocity, C0, 
and the oscillatory non-dimensional pressure amplitude desirecLat the injector 
face (AP/P), DELP. The variables GAM0 and C0 should be evaluated at the overall 
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UNBURNED PROPELLANT 


= 2.318 IN. (0.059 m) 

OX 

TEST #12 FUEL ox " 0-00^69 SEC 



DISTANCE FROM INJECTOR, METERS 

f4gure 14. Percent Unburned Propellant as a Function of Distance 
from Injector Face 
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mixture ratio. The variable DELP will scale the amplitude of the oscillatory 
waves solved for in subroutine CHAMDY. A value of 0.1 is recommended. 

The first variable input on the fifth card of the MAIN program control section 
input is NRf)0T. The absolute value of this variable controls the maximum 
number of frequency solutions the program will try to find. The program begins 
its search at the frequency impl ied by the input variable 0MEGAR. It will stop 
looking once it has found |NR00T| solutions or if the frequency is above FRQMAX. 
The variables IWRT and IWSKP are the next two variables input on this card. 

They control the amount of intermediate output one desires. Their exact 
function is described in the Program Operation Section of this manual. 

The last three variables input on the fifth card are KNTMX, KNTRMX, and KNTSMX. 
Their meaning is explained in Table 2 . Recommended values for these variables 
are 50, 100, and 20 respectively. 

The sixth card of the MAIN program control section input contains the variables 
0MEGAR, 0MEGAI, FRQMAX, DELFRQ, DELMX, and CTEST. The first two of these 
variables specify the starting guess in the u plane for solution. No solutions 
will be sought below the frequency implied by OMEGAR. Note that 0MEGAR and 
0MEGAI can be input as the frequency in Hertz times 2 ir and the growth co- 
efficient or the frequency in Hertz and the decrement depending upon whether 
NROOT was input as positive or negative. The variable FRQMAX, as mentioned 
earlier, is the maximum frequency allowed for the search algorithm to find 
solutions. The variable DELFRQ specifies the "stepsize" used by the search 
algorithm. Since there are sometimes many areas in the u plane which 
contain solutions, a fairly small stepsize is recommended, e.g., 5 Hz. The 
variable DELMX controls the maximum allowable change in the growth coefficient 
during successive iterations to minimize the error in the nozzle admittance 
boundary condition as a function of the imaginary part of w. A recommended 
value for this variable is 50 sec'^ . The last variable on this card, CTEST, 
is the upper bound on the condition number of the transposed Jacobian used to 
solve the nozzle admittance boundary condition. If the calculated condition 
number exceeds CTEST, then the search algorithm assumes that there is a sing- 
ularity near the current value of u and hence, does not proved further in 
that area to try and find a solution. A value of 50 to 80 is recommended. 
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The seventh card of the MAIN Program Control Section input contains error 
tolerances used in solving the nozzle admittance boundary condition. The first 
two, EPSF and EPSX, are used during the search algorithm and should be fairly 
large, e.g. 0.01 to 0.05 (1% to 5% error). The last two, EPSFS and EPSXS, 
control the final stages of iteration and should be fairly tight, e.g. 0.0005 
(0.05%). 

The eighth card of this control section contains the variables PC, MB0XI, and 
MBFUI. The first is the steady state chamber pressure, in PSIA, and the next 
two are the oxidizer and fuel injection flowrates, respectively (Ibm/sec). 

The last card, card number 9, in the MAIN Program Control Section input should 
be input only if the variable INPN0Z is less than or equal to three. If this 
is the case, the program needs to know the size and range of the frequency table 
it will use to generate tables for linear interpolation as discussed in the 
section describing the input variables INPHYD, INPC0M, and INPN0Z. The input 
variables on this ninth card are NFREQT, FREQMI, and FREQMX. Their meaning 
is described in Table 2. 

The next control section to read data after the MAIN program is the Nozzle ' 
Admittance Program. The data for this control section should be input only if 
INPN0Z<3. Otherwise, the information is not needed since the nozzle admittance 
information will be on tape ITAPN. Even when INPN0Z < 3, there is only one 
card input. This card contains information describing the nozzle geometry. 

Refer to Table 2 to determine the meaning of the variables on this card. 

Figure 15 shows exactly what portion of the nozzle each variable is appli- 
cable to. 

The next control section that requires data is the Hydrodynamics Program section. 
This control section uses namelist input. The data for this control section are 
only input if INPHYD ^ 2. Otherwise, the hydrodynamics information will be on 
tape ITAPH. The Hydrodynamic Input section (page 102) describes the meaning of 
the variables to be input for the control section. 

The last control section to require input data is the Combustion Dynamics 
Program. This input is contained on five cards. It should be omitted if 
INPC0M is greater than or equal to three, since then the combustion dynamics 
information will reside on tape ITAPC. 
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The first two cards for this section contain information specific to the oxidizer; 
the next two cards contain information specific to the fuel; the last card con- 
tains information specific to the combustion gas. Much of the information 
required in this section is obtainable from the output of the DER program 
(Ref. 11). This program solves for the steady state behavior of the chamber 
under consideration. The single stream tube option of the DER program is 
sufficient for this application. 

As mentioned above, the variables input on the first two cards of the Combustion 
Dynamics Program control section are applicable to the oxidizer. These variables 
are XK0X, TAUB0X, VB0X, DELH0X, TDRAG0, ADV0X, ADD0X, DELV0X, NUB0X, DT0XDM, and 
XIMP0X. The first variable XK0X, is the Klystron constant for the oxidizer jet. 
This variable controls the distance downstream from the injector face that the 
Klystron effect will be allowed to occur. The exact method of calculating XK0X 
has not been determined but it is recommended that a distance corresponding to 
approximately 45% of the oxidizer vaporized be used. 

The next variable, TAUB0X, is the steady state oxidizer vaporization time delay. 
Figure 14 gives the percent unburned for the sample input case of liquid 
oxidizer and fuel remaining in the chamber as a function of axial distance plotted 
on semi log graph paper. This output was derived from a DER computer run. The 
reciprocal of the average slope of the oxidizer curve during the initial burning 
phases should be taken as the oxidizer vaporization distance delay. For Figure 14, 
this is 2.318 inches. The time delay is obtained by dividing this by the 
average injector velocity, which, for the sample case presented in Figure 14 
is 28.859 ft/sec. The result is 0.00669 seconds. The next variable, VB0X, is 
the average oxidizer injection velocity. This is also given by the DER program 
from injector orifice and steady state AP considerations. 

The variable DELH0X is the pseudo energy of the oxidizer droplets and is dis- 
cussed at the end of this section. The variable TDRAG0 is the steady-state 
oxidizer drag time delay and a value of zero is currently recommended. If 
this variable is different from zero, the computed pressure and velocity 
profiles calculated in the chamber are unrealistic (see Conclusions and 
ReconmendaTions in Ref. 1). 
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ADV0X and ADD0X are the velocity and orifice diameter exponents for the oxidizer 
atomization process, i.e.. 



ADV0X 

V 


orif 


ADD0X 


(216) 


These variables are obtained from cold flow tests and for like-doublets -0.75 
and 0.57 are recommended. It should be noted that these variables are on 
different input cards. The next input variable is DELV0X which is the steady- 
state velocity difference between the oxidizer droplets and the combustion 
gas stream normalized to the sonic velocity. Since the vaporization rate is 
highest near the injector face, this variable should be evaluated near the 
injector face. Because the gas and the droplet velocities are approximately equal 
to each other at this location, a value of 0.01 is recommended based on turbulence 
levels in a rocket chamber (Ref. 26). 

The next variable, NUB0X, is the steady state oxidizer Nusselt number used in 
the vaporization expression. This variable should be computed from the relation 

0 2 

N., = 8/6 “^ox 

%x — ; (217) 

Sx ox 

where is the average oxidizer droplet diameter near the injector face, k' is 
the vaporization coefficient used in the k'-model evaluated at a mixture ratio 
near the injector face and x is the oxidizer time delay. Both D and k' 
are obtainable from a DER computer run. For the sample case input, "near the 
injector face" was taken as the first axial step printed by the DER program. 


The variable DT0XDM is the change in the oxidizer vaporization time delay with 
respect to mixture ratio. At the present time, a value of zero is recommended 
based on model verification cases (Ref. 1). 

The last oxidizer variable input is XIMP0X. This is the impingement point in 
inches for ^e oxidizer jets. 
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The next two cards contain data for the fuel. These data are obtained the same 
way as they were for the oxidizer and are input in the same order. 


The last card for this control section contains combustion gas data. The first 
variable input on this card, MWG, is the steady state molecular weight evaluated 
at the overall mixture ratio. This may be obtained from tables of molecular 
weight versus mixture ratio printed by the DER program. The next variable, CS, 
is the characteristic velocity evaluated at the overall mixture ratio. This is 
also obtainable from DER table output. The last three variables on this card 
are DRGDMR, DCSDMR, and DHDMR. These are the partial derivatives with respect 
to mixture ratio of the gas constant, the characteristic velocity, and the 
reference enthalpy respectively. DRGDMR and DCSDMR can be calculated from 
equilibrium calculations but a value of zero is recommended for DCSDMR based on 
model verification cases (Ref. 1). 


The variable DHDMR, and also the variables DELH0X and DELHFU, is calculated by 
curvfitting the steady-state energy equation with stagnation temperature/mixture 
ratio data calculated by an equilibrium program. The steady-state energy 
equation can be written as 








- (4 hox> 

r ‘ . I . 


- I 


MR (1+MR) - MR^ (1 


where 


+MR.) / 3h \ 

A 




(218) 


(219) 


and the subscript <t> indicates that the variable is to be evaluated based on 
the overall mixture ratio. 
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HYDRODYNAMIC INPUT 


This section describes data needed by the hydrodynamics subroutine, HYDRDY, 
to simulate the various feed system components. It is assumed that the 
feed system being modeled has been laid out on the generalized feed system 
schematic of Fig. 4 with an appropriate segment number assigned to each 
feed system component (or combination of components). 

Basic Feed-System Data 

To describe the basic feed system it is necessary to know the length, area, 
resistance and wall compliance of each of the numbered segments of Fig. 4 
which are being used. Also, the acoustic velocity and density of the 
fluid in each segment must be known. If there is dissolved an entrained 
gas in the system, then a preliminary calculation must be made for each 
feed system section to account for the effect of the gas on the fluid 
acoustic velocity. 

Specific parameters required for the numbered segments are: 

A - Segment cross-sectional area - in.^ 

CW - Segment wall compliance (AV/AP/V - in.^/lb 
L - Segment length - in. 

R - Segment linearized hydraulic resistance (aP/W) - sec/in. ^ 

V - Segment fluid acoustic velocity - in. /sec 
RH0L - Segment fluid density - Ib/in.^ 

Valves, Fittings, Orifices, Screens. Flowmeters, etc. 

These components can each be described in the model simply as lumped resist- 
ance at the end of a line segment. Rather than using all the attributes of 
one of the numbered segments (length, area, wall compliance, etc.) for one 
of these "resistance only" components, it is suggested that its resistance 
merely be added to that of the adjacent upstream pipe segment. The combin- 
ation canHhen be entered as one of the numbered segments -w4th the length, 
area and wall compliance values being primarily those of the pipe segment. 
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Accumulators 

A feed system accumulator can be represented as one of the side branch 
lines of the Fig. 4 schematic by specifying an appropriate length, area, 
acoustic velocity and fluid density for the fluid volume of the accumu- 
lator and also specifying an appropriate connecting resistance. The 
spring rate of the accumulator piston can be specified in terms of the 
segment wall compliance value. 

Propellant Tanks 

A large tank will have the effect of constant fluid pressure at its 
outlet and can be represented simply as the input to segment number 1, 

15, or 22. No descriptive parameters are required for these inputs. 

Small tanks can be represented as one of the side branch lines in a 
manner similar to an accumulator. Ullage volume in a small tank is 
represented by a reduced value for the segment acoustic velocity. 

Cavitatinq Venturi es 

The steady-state effect of a cavitating venturi is to have constant 
flow through the venturi as a function of variations in downstream 
pressure. For an oscillatory system, the vapor bubble downstream of 
the venturi throat makes the venturi look like a constant pressure 
boundary for small amplitude oscillations. To simulate this effect the 
effective acoustic velocity for the segment downstream of the cavitating 
venturi should be made very small (*10 inches/sec would be appropriate). 

The steady-state hydraulic resistance of the cavitating venturi can be 
lumped with that of the upstream pipe segment as described above for 
valves, fittings, etc. 

Reqeneratively-Cooled Thrust Chamber 

Regeneratively-cooled thrust chamber jackets can be represented as one or 
more of the numbered Fig. 4 segments. Because in most thrust chambers the 
coolant flow area changes continuously with length, as many segments as 
possibla- should be devoted to the jacket so as to improve- the simulation 
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accuracy. The fluid temperature also may change significantly along 
the chamber length thereby necessitating the use of several segments 
with different acoustic velocities to achieve accurate simulation. 

Lines, Ducts. Bends. Bellows, and Flex Lines 
These components are described in the model in terms of the basic 
numbered segment input parameters of length, area, fluid acoustic 
velocity, fluid density, wall compliance and linearized hydraulic 
resistance. For a duct or line of constant diameter, D, wall thickness, 
h, and wall material bulk modulus, E, the program input wall compliance 
value, CW, is simply D/LE. For a bellows or flux line of volume, V, the 
wall compliance value, CW, may be calculated from AV/Ap/V where AV/AP 
is the volume change per psi at the operating pressure. 

Injectors 

The hydrodynamics subprogram employs a separate set of equations to 
describe the hydrodynamic characteristics of the two injectors in the 
Fig. 4 generalized feed system schematic. The specific input parameters 
for the two injectors are the volume, linearized hydraulic resistance, 
orifice inertance (1/Ag), fluid acoustic velocity, and a structural parameter 
defining the change in injector volume per psi of injector aP. In terms 
of the program variable names the required injector parameters are: 

RF - Linearized hydraulic resistance for the "F" injector, 

(2 Ap/w) _ sec/in. ^ 

R0 - Linearized hydraulic resistance for the "0" injector, 

(2 AP/w) - sec/in. ^ 

VOLF - Volume of the "F" injector - in.^ 

VOLCI - Volume of the "0" injector - in.^ 

VF - Fluid acoustic velocity of the "F" injector - in. /sec 

V0 - Fluid acoustic velocity of the "0" injector - in. /sec 

ZF - Inertance of the "F" injector orifices (1/Ag) - sec^/in.^ 

Z0 — Inertance of the "0" injector orifices, (1/At) - sec^/in.^ 

KF - "F" injector deflection constant, (AV/Ap) - in.^/lb 

KQ - "(J" injector deflection constant, (AV/Ap) - in.^/lb 
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Tees, Splitters and Capped Lines 

No provision is made in the model for completely generalized input of tees 
and branched lines. However, a system of considerable complexity can be 
modeled by laying out an appropriate flow path on the generalized Fig. 4 
schematic. For example, a feed system with up to seven side branch lines 
can be simulated by choosing the main flow path through segments 1, 3, 4, 

14, 17, 21, 25, 26, 28, 29, 30 and the"F" injector in series. 

Input Variables 

Data input to the hydrodynamics subprogram is from three sources: (1) Via 

the argument list in the CALL HYDRDY statement, (2) Through labeled coniion 
block/C0MTAP/and (3) By use of the NAMELIST data read routine. 

The argument list variables, in order, are: 

IR - Data read flag - dimensionless 

INPHYD- Program function flag - dimensionless 

FRE - Single frequency for feed system frequency response 
calculation - Hz 

GIN0 - Output value of oscillatory oxidizer injector flowrate for 
input frequency FRE - dimensionless 

GINF - Output value of oscillatory fuel injector flowrate for input 
frequency FRE - dimensionless 

PCIN - Injector end thrust chamber pressure - psia 

W0IN - Steady-state oxidizer injector flowrate - Ib/sec 

WFIN - Steady-state fuel injector flowrate - Ib/sec 

Several HYDRDY input variables are transmitted via labeled common block/ 
C0NTAP/. ITAPH is the logical unit number of the output device on which sub- 
routine HYDRY tabulates output values of oscillatory injection flowrates 
for the specific frequencies (up to 100 separate values) given in the array 
FREQT in common block/CONTAP/ . The value of the FREQT in common block/C0NTAP/ 
is the total number of frequencies stored in the array FREQT. 

All other data required by HYDRDY, including all the feed system descriptive 
data, Ts read in by use of the NAMELIST routine. The local rules for using 
this routine should be checked to verify that the correct card or terminal 
format is being used. Table 3 shows a list of allowable FORTRAN names, the 
maximum values of subscripts, and a definition of the names. The name of 
the NAMELIST block is HYD. 
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TABLE 3. NAME LIST/HYD/DATA INPUT NAMES 


NAME 

A(30) 

CW(30) 


DEFINITION 

Array containing segment cross-sectional area values 
Array containing segment wall compliance values 


FREQ 

FREQT(IOO) 

ICRT 

ID 

IH(126) 


Frequency at which HYDRDY will compute feed system 
frequency response if INPHYD ;< 1 

Array containing frequencies at which HYDRDY will compute 
feed system freouency response of INPHYD > 1 . 

ICRT = 1; injector flowrate gain and phase will be plotted 
vs frequency. ICRT = 0; no plot (default) 

Dummy name to allow for data card sequence numbers 

Array containing control flags used by subroutine 
FREQD. Can be used to obtain printouts and plots of 
feed system frequency response for other variables in 
addition to the injector flowrates. 


IRFLAG 

ITERM 

ITYPE 


IWRITE 


KF 

K0 

L(30) 


IRFLAG = 0; read data from unit 5 (default) 

IRFLAG f 0; read data from unit ITERM 

ITERM = 0; no terminal data input (default) 

ITERM > 0; read data from terminal (unit ITERM) 

ITYPE = 1; both oxidizer and fuel feed systems are 
modeled simultaneously (default) 

ITYPE = 2; oxidizer feed system modeled on first pass 
of frequency response routine; fuel feed system modeled 
on second pass. 

IWRITE = 0; HYDRDY input printed on unit 6 (default) 
IWRITE = 2; extensive printout of all HYDRDY inout, 
intermediate output and final output on unit 6 
IWRITE = 1; printout of HYDRDY input and final output 
on unit 6 

Injector face flexibility constant for "F" injector 
Injector face flexibility constant for "0" injector. 
Array containing segment length values — 
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TABLE 3. (Concluded) 


NAME 

DEFINITION 

0MI 

Lowest frequency for injector flowrate gain/phase plot; 
default to FREQT(l) if not entered 

0MFL 

Highest frequency for injector flowrate oain/phase plot; 
default to FREQT(NFREQT) if not entered 

R(30) 

Array containing segment hydraulic resistance values 

RF 

Hydraulic resistance for "F" injector 

RH0L(3O) 

Array containing segment propellant density values 

R0 

Hydraulic resistance for "0" injector 

V(30) 

Array containing segment acoustic velocity values 

VF 

Acoustic velocity for "F" injector 

V0 

Acoustic velocity for "0" injector 

V0LF 

Volume of "F" injector 

V0L0 

Volume of "0" injector 

ZF 

Inertance of "F" injector 

Z0 

Inertance of "0" injector 
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Required input in the NAMELIST/HYD/data is a value of A, CW, L. R, RH^L, 
and V (see Table 3 for descriptions) for each numbered segment being in- 
cluded in the feed system. Values of KF and/or K0, RF and/or R0, VF and/ 
or V0, V0LF and/or V0L0 and ZF and/or Z0 are also required. It should be 
noted that, when possible, both oxidizer and fuel feed systems should 
simultaneously be laid out on the Fig. 4 schematic with the injector labeled 
"V being used for the oxidizer side (data values K0, R0, V0, V0L0 and Z0) 
and the injector labeled "F" being used for the fuel side (data values, KF, 

RF, VF, V0LF and ZF). If this can be done, a single call to subroutine 

FRESP will generate frequency response data for both oxidizer and fuel feed 
systems. If feed system complexity requires that the oxidizer and fuel feed 
systems be laid out separately on the Fig. 4 schematic, then two sets of 
input data must be read and subroutine HYDRDY must call FRESP twice - first 
for the oxidizer feed system calculations and second for the fuel feed sys- 
tem calculations. To specify this option, variable ITYPE must be set equal 

to 2. 


Variable INPHYD in the HYDRDY argument list controls the HYDRDY calculation 
process. If INPHYD < 1 the oscillation injection flowrates are calculated 
for a single frequency, specified by variable FRE in the HYDRDY argument list. 
If INPHYD > 1, HYDRDY calculates oscillatory injector flowrates for the num- 
ber of frequencies, NFREQT , which are contained in array FREQT. Both NFREQT 
and the array FREQT are stored in labeled common block/C0MTAP/ prior to call- 
ing HYDRDY. 

NAMELIST variable ICRT controls the option for generating CRT plots of the 
oscillatory injection flowrate gains and phase values as a function of fre- 
quency. If ICRT = 0 (the default value) no plots are made. If ICRT > 1 
plot output is written to the output file named SYSCRT. 

NAMELIST variable, ID, is a dummy name which can be used on each input card 
to provide an identification number field without violating the NAMELIST 
restriction that the entire card is read as data. For example, ID = 00000010 
could be^in columns 70-80 of a HYDRDY data card and ID = 00000020 in columns 
70-80 of the next card. The NAMELIST routine would then interpret each card's 
sequence number as a new value for the dummy variable, ID. The value of ID 
is not used in any way by subroutine HYDRDY. 
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NAMELIST variable, ITYPE, is used to indicate to HYDRDY the format of the 
feed system modeling. If ITYPE = 1 (the default value) it is assumed that 
both oxidizer and fuel feed systems are modeled simultaneously with only 
one set of HYDRDY input values (for the 30 segments and 2 injectors of the 
Fig, 4 schematic). If ITYPE = 2, HYDRDY will send two consecutive sets of 
input data; the first set will be assigned to the oxidizer feed system and 
the second set to the fuel feed system. For either value of ITYPE the pro- 
gram will assume that the injector labeled on Fig, 4 is the oxidizer 
flow outlet and the injector labeled "F" is the fuel flow outlet. There- 
fore, this convention must be followed when laying out the feed system model. 

NAMELIST variables IRFLAG and ITERM are optional HYDRDY inputs which in- 
dicate that data input will be provided from a timesharing terminal. If 
IRFLAG = 0 (the default value), input data will be read only from FORTRAN 
logic unit number, ITERM. It should be noted that the default values for 
IRFLAG and ITERM are set up so that the initial data input will always be 
card input on unit 5. 

After reading the initial NAMELIST data on unit 5, HRDRDY checks the value 
of ITERM and, if non-zero, proceeds to read additional first case data 
from the terminal on unit ITERM, Thus, the first case card NAMELIST input 
could consist of the single item ITERM = N, where N is the terminal logical 
unit number. If terminal input only is desired, block data program/F/ can 
be recompiled with the IRFLAG default value changed from 0 to 1 and the ITERM 
default value changed from 0 to the desired unit number. 

NAMELIST variable IWRITE (main program control variable IPRHYD) controls 
the printed output from HYDRDY. If IWRITE = 0 (the default value), only 
the NAMELIST input to HYDRDY is printed on logical unit 6. If IWRITE < 0 
the NAMELIST input is printed and the normal HYDRDY output is printed as 
well as being saved on an output device in binary form. If IWRITE = 1 both 
HYDRY input and normal output are printed. If IWRITE > 1 extensive print- 
outs of subroutine FRESP intermediate calculations are printed in addition 
to the normal HYDRDY input and output. __ 
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PROGRAM OUTPUT 

The output of the FSCSM computer program is provided as the usual 
tabular printout. A sample case is included in Appendix E which 
corresponds to the input dataset listed in Appendix D. As is also 
mentioned in the Program Operation section of this manual, the in- 
put case listed in Appendix D consists of two cases being run back 
to back. The output from the first case is given in Appendix E 
from pages E-2 through E-14. The first page of output consists of 
a title page identifying the current version of the FSCSM computer 
program. The input data are printed out as they are read in. This 
permits both a full documentation of the computer run conditions for 
later analysis as well as a convenient method to check for input 
errors if unusual results are calculated. Page E-3 of the listing 
in Appendix E gives the two alphanumeric cards identifying the case 
at the top of the page right under the program title. Subsequent 
to these two cards, the information on the cards read in by the main 
control section and the nozzle admittance control section are printed 
out. After reading and writing these cards, and since INPN0Z = 3, 

• the program proceeds to the nozzle admittance table calculations. 
Information pertinent to these calculations is printed on page E-4. 

The freouency table goes from 150 Hz to 400 Hz as specified by the 
input variables FREOMI and FREQMX. 

Since INPHYD = 2, the program proceeds to the hydrodynamic subroutines 
right after the nozzle admittance calculations. Input for this routine 
is in the form of NAMELIST data. The NAMELIST is output on pages E-5 
and E-6. A printout of the feed system response table computed by sub- 
routine HYDRDY and saved on file ITAPH is given on page E-7. 

The next set of input required is used in subroutine C0MBOY and STEADY. 
This is output on page E-8. The steady-state profiles are then computed 
and printed on pages E-9 and E-10. 
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The program then begins its search for solutions to the nozzle admit- 
tance boundary condition. The first one it finds is at 210.42 Hz. 

The program then outputs the combustion dynamic coefficients, the fre- 
quency and decrement, and the feed system response for this solution 
on page E-14. On page E-12, the oscillatory profiles correspond 
to this solution are given. The program then proceeds to the next 
case. 

Since the second case does not generate the data on files ITAPN and 
ITAPH (it only reads this information), these tables are not printed. 

The first page of output, page E- 15, in Appendix E, consists entirely 
of the data read by the Main Control Section and the Combustion Dynamics 
Control Section. Since the STEADY Control Section print code is zero 
(IPRSTE = 0), the program skips over the steady-state output (although 
of course, it still computes it) and proceeds directly to the section 
which solves the nozzle admittance boundary condition. The first root 
it finds above the input frequency of 265 Hz (given by the variable 
0ME6AR), is at 280.62 Hz. It prints out the frequency, decrement, 
nozzle admittances, and feed system response for this solution. Out- 
put of the combustion coefficients and oscillatory profiles is bypassed 
because the input flags IPRC0M and IPRCHM were set to zero. 

The final page of output is the title page. This indicates normal 
termination of the job. 
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PROGRAM OPERATION 

The FSCSM computer program is designed to read in an input case sequen- 
tially, perform the calculations for that case, and output the results. 

The program then transfers back to its beginning to read in the next 
case. In this manner, running jobs back-to-back is quite straightforward. 
The sample input case listed in Appendix D provides an example of two 
such cases run back to back. The first case, given by the first 24 cards, 
is run with no prior information residing on the hydrodynamic feed system, 
the combustion dynamics, or the nozzle admittance datasets. Since 
INPHYD = 2 and INPN0Z = 3, this case generates tables of the hydrodynamic 
feed system response and nozzle admittance versus frequency and saves them 
on files ITAPH and ITAPN, respectively. The subseouent case (the last 
13 cards in Appendix D ) will use the information stored on these datasets. 
Although these two cases were run back to back, this was by no means 
necessary. The second case is self-contained and could be submitted 
separately. Of course, if this were the situation, the user must be 
sure there are datasets on files ITAPH and ITAPN which are applicable to 
that second case. 

For the sample dataset run, the two input cases found solutions to the 
nozzle admittance boundary equation at 210.42 Hz and 280.62 Hz, respectively. 
If there are no other solutions between these two frequencies, the same 
effect could have been obtained by setting the input variable NR00T equal 
to -2 for the first case instead of -1. The program would have then looked 
for the first two roots above the input frequency 190 Hz (0MEGAR) and 
found both solutions automatically. The second case would not be input 
for this situation. 
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Program Size, Overlay Structure, and Timing 

Without overlay, the FSCSM computer program load module requires 
262.4 K Bytes of computer storage on the IBM 370 Model 165 computer. 

This storage does not count the buffers needed for input/output. 

If one allocates a 1 K Byte of buffer size for each of the three 
date sets used to store the feed system, combustion dynamic, and 
nozzle admittance data (which are all unformatted input/output) , 
uses two buffers for each data set, and adds in the buffer require- 
ments for his card input, printed output, and CRT output, then the 
total buffer soace should be well under 10 K bytes on a 370/165 
computer. With the overlay structure specified in Fig. 16 , the 
total program reouirement is 220 K bytes of storage on an IBM 
370/165 computer, including two buffers for each of the three 
unformatted datasets at 1 K bytes each. 

Comouter run time has only been checked for an IBM 370/165 computer 
where the subroutines were compiled using the IBM procedure AF0RTRAN 
with the optimizing parameter, 0PT, equal to one. For this situation, 
each iteration during the search algorithm portion of the program 
(when ISCNT equals one or four) averaged 3.7 CPU seconds. When ISCNT=5, 
each iteration is about twice as fast. For the cases run during model 
verification a five Hz step size for the search algorithim was used 
(DELFRQ=5). For these cases, each solution to the nozzle admittance 
boundary equations averaged 0.85 minutes of CPU. 

Program Input/Output Dataset File Information 

The case input dataset file number used by the FSCSM computer program 
is 5. The printed output dataset file number is 6. There are three 
auxiliary files used by the program. These are specified by the input 
parameters, ITAPC, ITAPH, and ITAPN, corresponding to the combustion 
dynamics, hydrodynamic feed system, and nozzle admittance datasets. 
Control of the reading from or writing on to these respective datasets 
is specified by the three input flags INPC0M, INPHYD, -^d INPN0Z. The 
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program uses unformatted input/output statements for transmitting 
information to and from these datasets. A convenient blocksize to 
use is 1 K bytes. 
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Figure 16. FSCSM Program Overlay Structure 
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Diagnostics 

The Feed System Coupled Stability Model computer program has been designed 
to operate as straightforward as possible with a minimum amount of user 
interaction for each case being run. There may be times however, when the 
program’s results appear questionable or the algorithm used to find solutions 
in the frequency space to the nozzle admittance boundary equation runs into 
difficulty or does not find solutions that were expected. Many diagnostic 
messages are coded into the program to warn the user of such problems. 

Also, there are certain dump codes which enable the user to obtain inter- 
mediate output in order to debug most problems that may arise. 

One of these dump codes is the variable IWSKP. When this is set to zero, 
no intermediate output is obtained. When it equals one, a certain amount 
of limited output will be generated. This output comes in two forms depend- 
ing upon whether or not the program is within its search algorithm portion 
or its two-dimensional secant portion. For the first case, the variable 
ISCNT has the value of one or four. In the second case it has the value 5. 

When ISCNT equals one or four and IWSKP equals one, subroutine SOLVW will 
print the following variables in the order given; the iteration counter 
(KNTR), the control flag (ISCNT), the counter (KSCNT4), the current values 
of omega (oj), the upstream and downstream nozzle admittances (CN0ZA and N0ZA), 
the absolute value of the error in the nozzle admittance equation (HN), the 
value of the test function (FTSTZ), and the determinant and condition number 
of the transposed Jacobian (DET2 and C0ND2). This printing will be performed 
every time the real part of w is incremented by 2tr*DELFRQ right after the 
imaginary part of uj has been chosen to minimize |F|, the absolute value of the 
error in the nozzle admittance boundary equation. The user can employ this 
output to determine if there is a region in the w-plane where a possible solu- 
tion may have existed (e.g., the error became small but the test function did 
not change sign). He can then rerun his case while taking smaller frequency 
steps through the narrowed range where he suspects a solution may exist. Also, 
the program may jump over a solution if there is a singularity within DELFRQ 
of that %Q.lution. If this is the case, the program will 5 £nse the singularity 
and not proceed any further in its search in that range. Rerunning the case 
with a smaller value of DELFRQ will solve this problem. 


R-9808/115 



When ISCNT equals five and IWSKP equals one, subroutine S0LVW prints after 
each two-dimensional secant method, these variables in the following order: 
ISCNT, KNTS, 0MEGA, CN0ZA, N0ZA, and HN. Although it did not happen for 
any of the cases performed during the checkout of the computer model, the 
two-dimensional secant method may diverge. The above computer output would 
be useful in determining the cause of the problem. 

When the variable IWSKP equals two, all the above output is printed plus the 
followi ng: 

1. When ISCNT equals one or four, intermediate output is obtained dur- 
ing the iterations to minimize |FN| with respect to Imag (w). 

For this case, one obtains the variables KNT, lER, XI, X2, FI, F2, 
0MEGA, FN, GN, and HN. These are all described in Appendix A. 

This output may be useful in seeing how the error is changing as 
a function of the decrement when the real part of the frequency is 
held fixed. Further, when ISCNT equals one or four, the variables 
FN, DFRDX, DFIDX, DFRDY, and DFIDY are printed along with the out- 
put obtained when IWSKP = 1. 

2. When ISCNT equals five and IWSKP equals two, one obtains the out- 

put for the case IWSKP equals one for the two-dimensional secant 
method plus the following variables in order: XRl , XII, FRl , FIl, 

XR2, XI2, FR2, FI2, XR3, XI3, FR3, FI3, XR4, XI4, FR4, FI4, and FN. 
These variables correspond to the current values of w and FN being 
used by the two-dimensional secant method. They can be used to 
trace which points the algorithm is replacing as the iteration 
proceeds as well as how the error is behaving. 

Another input variable which controls intermediate output is the FORTRAN 
variable IWRT. This variable is input as zero, no intermediate output is 
obtained. If this variable is input as a positive number, then intermediate 
output from subroutine CHAMDY is obtained. This output consists of the oscil- 
latory profiles for the variables P (pressure), RH0 (density), MR (mixture 
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ratio), and T (temperature) along with the current value of the complex 
frequency, oj. This output is printed every time subroutine CHAMDY is 
entered . 

The diagnostic messages that are coded within the FSCSM computer program 
may be printed for several reasons. 

Within subroutine S0LVW, there are three diagnostic messages coded which 
will appear when certain iteration counters are exceeded. The first is 

****^*******ir:k*******'k'kic**iriir**ir*itit*it**it**ic**i^i^**icic 

WARNING, POSSIBLE ROOT IN FREQUENCY RANGE: — 
************************************************* 


When this message appears, it means that a potential root was bracketed but 
the error did not decrease sufficiently within ten additional iterations 
to warrant the program proceeding further with its search in that range. 
Moreover, the determinant did not change sign and the condition number 
remained less than CTEST in that range. Rerunning the case over the speci- 
fied frequency range given in the message with IWSKP equal to one or two 
may prove beneficial if the user suspects there may be an actual solution 
in that range. 

The second diagnostic message printed by subroutine S0LVW is 

**** UNABLE T0 FIND R00T F0R IMAG PART 0F F **** 

Along with this message, the variables XI, FI, K2, F2, X3, F3, ANS, FANS, 

KNT, lER, and 0MEGA are printed in the order listed. When this message 
appears, it means the algorithm to minimize |FN| with respect to Imag (w) 
has failed. If this message appears, it usually means something is wrong 
with the input parameters. The only occurrence that the programmers are 
aware of when this is not the case is when the error attains a minimum as 
I Imag (o))_|-* ». Since this happens only in the most extraordinary situations, 
the procedure should be to rerun the case and not include the frequency 
range where that anomally is occurring. 
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The third diagnostic produced by subroutine S0LVW is 

**** EXCEEDED CONVERGENCE LIMIT **** 

Along with this message, the variables lER, KNTS, KNTR, ISCNT, XRI, XII, 
FRl, FIl, XR4, XI4, FR4, FI4, are printed. 

This message will appear if KNTR is greater than KNTRMX or KNTS is greater 
than KNTSMX. In the former case, the usual error is that the user input 
too small a DELFRO to cover the range between solutions to the nozzle ad- 
mittance boundary equation in KNTRMX stens or too small a KNTRMX to allow 
that range to be covered in steps of length DELFRO. 

In the case where KNTS is greater than KNTSMX, it would probably mean that 
the two-dimensional secant method is diverging. The job should be rerun 
with IWSKP equal to one or two to obtain more information concerning the 
problem. 

There is also a diagnostic message printed from subroutine C0MMAT. This 
is the subroutine that solves the four by four system of linear equations 
for subroutine CHAMDY. If any of the diagonal elements of the associated 
matrix are zero, then the message 

***** DIVIDE CHECK IN C0MMAT ***** 

will appear along with a printout of the row number of the zero diagonal 
as well as the complex matrix being solved. If this error messaae appears, 
then there must be something very wrong with the case being run, e.g, the 
input data is in error, or a dimension has been exceeded. One should 
recheck his input carefully and then, if necessary, rerun the case with 
IWRT eoual to one and IWSKP equal to one or two. 

Two similar messages as the one above are printed by subroutine C0GAEL. 

The fixst of these messages is 
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****ERROR IN C0GAEL SUBROUTINE, J AND JfWX EQUAL, RESPECTIVELY**** 


and the second is 

****MATRIX IS SINGULAR, EXIT FROM C0GAEL. THE PIVOT ELEMENTS ARE... 
The reasons for these errors are similar to the C0MMAT error message. 
Conversion to UNIVAC 

The following cards must be changed to execute the FSCSM program on a 
UNIVAC computer (see Appendix C for code listing): 

1. Change C0MPLEX*16. . .to C0MPLEX 

Routi ne Card Number 

CHAMDY 80 

C0MMAT 150 

2. In subroutine C0MMAT, change CDABS to CABS on card No. 24. 

3. In subroutine TDPL0T, replace card 11310 with 

31 CONTINUE 0011310 

4. In subroutine N0ZADM, replace card 1730 with 

8 N0ZA = CMPLX(SYR, SYI )... .00001 730 

5. In the main program, change ATAND( ) to 57.296*ATAN( ) 

on card numbers 3370, 3410, 3450, and 3490. 
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APPENDIX A 


FORTRAN SYMBOL TABLE 
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FORTRAN ENGINEERING 


VARIABLE 

NAME 

VARIABLE 

SYMBOL 

TYPE 

CONTROL 

SECTION 

DESCRIPTION 

A(IOO) 

A 

R 

C0MARE 

Axial distance array, 

X(I) = XO + (X-1)* DELX, m 

ADDFU 

®fu 

R 

C0MCBM 

Fuel liquid jet diameter exponent, 
unitless 

ADD0X 


R 

C0MCBM 

Oxidizer liquid jet diameter exponent, 
unitless 

ADVFU 

bfu 

R 

C0MCBM 

Velocity exponent for fuel atomization 
process, unitless 

ADV0X 

^0X 

R 

C0MCBM 

Velocity exponent for the oxidizer 
atomization 

AINJ 

n j 

R 

C0MARE 

Cross sectional area at injector face, 
m2 

AMA(4,5) 

^ — 

C 

CHAMDY 

Array used to store coefficients of 
finite difference equations approximat- 
ing the oscillatory differential 
equations 

ANGLEX 

— 

R 

C0MN0Z 

Nozzle convergence half angle, degrees 
(see Fig. 15) 

ANS 


R 

S0LVW 

Solution to minimization of error in 
nozzle admittance boundary equations. 
Returned to S0LVW from subroutine ZER0. 

CFUl 

CFU16 

C-i , . . . , 

'fu 

C 

C0MCBM 

Fuel combustion coefficients, unitless 

CMA(4) 


C 

CHAMDY 

Right hand side of finite difference 
equation approximating the oscillatory 
differential equations. In equivalence 
with CMA(1,5). 

CN0ZA 

''-'u 

C 

C0MCHM 

Upstream nozzle admittance 

C0 


R 

C0MCBM 

Sonic velocity evaluated at overall mix- 
ture ratio, m/sec (ft/sec) 

C0ND1 


R 

S0LVW 

Condition number of the transpose of the 
Jacobian of the difference between the 
upstream and downstream nozzle admit- 
tances with respect to w. 
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FORTRAN 

VARIABLE 

NAME 

ENGINEERING 

VARIABLE 

SYMBOL 

TYPE 

CONTROL 

SECTION 

DESCRIPTION 

C0ND2 

— 

R 

S0LVW 

Condition number of the transpose of 
the Jacobian of the difference between 
the upstream and downstream nozzle ad- 
mittances with respect to w. 

C0X1,.. 

C0X16 

• » c 1 

ox 

C 

C0MCBM 

Oxidizer combustion coefficients, 
unitless 

CRR 


R 

C0MN0Z 

Contraction ratio, cross-sectional area 
of chamber/throat area, unitless 
(see Fig. 15) 

CS 

c* 

R 

C0MCBM 

Characteristic velocity evaluated at 
the overall mixture ratio, m/sec (ft/sec) 

CTEST 

— 

R 

S0LVE 

Upper bound on the condition number of 


the transpose of the Jacobian of the 
difference between the upstream and 
downstream nozzle admittances with re- 
spect to the complex frequency. If the 
condition number of that matrix exceeds 
CTEST for a given frequency, then it is 
assumed that Jacobian is singular near 
that frequency and hence a solution will 
not be sought at that point, unitless. 


DA(IOO) 

dMdx 

R 

C0MARE 

Slope of area of chamber at XM(I), m 

DCSDMR 

dc*/dm 

R 

C0MCBM 

Partial derivative of the characteristic 
velocity with respect to mixture ratio 
holding, m/sec. 

DELFRQ 


R 

S0LVE 

Frequency increment used during the 
search procedure to solve the nozzle 
admittance boundary equation, Hz. 

DELHFU 


R 

C0MCBM 

Pseudo energy term for fuel , J/kg 
(Btu/lbm) 

DELH0X 

*''ox 

R 

C0MCBM 

Pseudo energy term for oxidizer, J/kg 
(Btu/lbm) 

DELMX 

— 

R 

S0LVE 

Maximum change in u allowed between 
iterations 

DELP 

AP 

R 

C0MCHH 

Oscillatory pressure at injector face, 
dimensionless ^ 
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FORTRAN 

VARIABLE 

NAME 

ENGINEERING 

VARIABLE 

SYMBOL 

TYPE 

CONTROL 

SECTION 

DELVFU 

— 

R 

C0MCBM 

DELV0X 

— 

R 

C0MCBM 

DELX 

Ax 

R 

C0MARE 

DETl 

— 

R 

S0LVW 

DET2 

— 

R 

S0LVW 

DFIDX 

-- 

R 

S0LVW 

DFIDY 

— 

R 

S0LVW 

DHDMR 

{dh/dm)^ 

R 

C0MCBM 

DM1 

— 

C 

ADARND 

DM2 


R 

ADARND 


DESCRIPTION 

Steady-state velocity difference be- 
tween fuel droplets and gas stream 
normalized to the sonic velocity at 
the overall mixture ratio, unitless 

Steady-state velocity difference be- 
tween oxidizer droplets and gas stream 
normalized to the sonic velocity at the 
overall mixture ratio, unitless 

Axial distance between successive X(I), 
m 

Determinant of the Jacobian of the dif- 
ference between the upstream and down- 
stream nozzle admittances with respect 
to 0 ). 

Determinant of the Jacobian of the dif- 
ference between the upstream and down- 
stream nozzle admittances with respect 
to 0 ). 

Derivative of Imag (FN) with respect 
to real (co). 

Derivative of real (FN) with respect 
to real (u). 

Partial derivative of gas reference 
enthalpy with respect to mixture ratio 
averaged over the mixture ratio range 
during steady-state operation, J/kg 
(Btu/lbm) . 

Used to store intermediate values needed 
to compute boundary conditions and co- 
efficients of finite difference equa- 
tions approximating the oscillatory 
differential equations. 

Used to store intermediate values needed 
to compute boundary conditions and co- 
efficients of finite difference equa- 
tion approximating the oscillatory 
differential equations. 
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FORTRAN 

VARIABLE 

NAME 

ENGINEERING 

VARIABLE 

SYMBOL 

TYPE 

CONTROL 

SECTION 

DM3 

— 

C 

ADARND 


DM4 

— 

c 

ADARND 

DM5 

— 

c 

ADARND 

DM6 

— 

c 

ADARND 

DM7FU 

— 

c 

ADARND 

DM70X 

— 

c 

ADARND 

DM8FU 

— 

c 

ADARND 

DM80X 

-- 

c 

ADARND 

DM9FU 

— 

c 

ADARND 

DM90X 

— 

c 

ADARND 

DM22 

— 

R 

ADARND 

DMRB(IOO) 

^MR/x 

R 

C0NSTS 

DRGDMR 

dR/dx 

R 

C0MCBM 

DRHQB(IOO) 

dQldx 

R 

C0NSTS 

DTFUDM 

St /StIR 
fu 

R 

C0MCBM 


DESCRIPTION 

Used to store intermediate values needed 
to compute boundary conditions and co- 
efficients of finite difference equa- 
tions approximating the oscillatory 
differential equations. 


Derivative of steady-state mixture ratio 
with respect to distance, m“' 

Partial derivative of gas constant with 
respect to mixture ratio evaluated at 
the overall mixture ratio, J/kmole/°K 
(ft-lb/lb/“R) 

Derivative of steady-state density with 
respect to distance, kg/m3/m 

Partial derivative of fuel vaporization 
time delay with respect to mixture ratio, 
holding the vaporization blockage term, 
droplet diameter, and Nusselt number 
constant, sec 
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FORTRAN 

VARIABLE 

NAME 


DT0XDM 


EPSF 


EPSFS 


EPSX 


ENGINEERING 

VARIABLE CONTROL 

— SYMBOL — TYPE SECTION DESCRIPTION 




C0MCBM Partial derivative of oxidizer vaporiz- 
ation time delay with respect to mixture 
ratio, holding the vaporization block- 
age term, drop diameter, and Nusselt 
number constant, sec 


S0LVE Relative error criterion used during 
the search algorithm for the portion 
of the program that minimizes the error, 
HN, in the nozzle admittance boundary 
equation with respect to the growth 
coefficient, coj. To obtain convergence, 
it is necessary that 


R 


g|HN|^ 

5 0)j 

Uni tl ess. 


<EPSFS 


Tightened relative error criterion used 
to determine if convergence has been 
obtained while iterating to solve the 
nozzle admittance boundary equation. 

To obtain convergence, it is necessary 
that ^ 


|Na - N. |/|N. I <EPSFS. 

Relative error criterion used during the 
search algorithm for the portion of the 
program that minimizes the error, HN, 
in the nozzle admittance boundary equa- 
tion with respect to the growth coeffic- 
ient, wj. To obtain convergence, it is 
necessary that 

. - Uj |/|o), 1 <EPSX 

n h h 

where the subscripts 1 and 2 refer to 
two successive iterations, unitless. 


R-9808/A-6 



FORTRAN 

ENGINEERING 





VARIABLE 

VARIABLE 


CONTROL 



NAME 

SYMBOL 

TYPE 

SECTION 

DESCRIPTION 

EPSXS 


R 

S0LVE 

Tightened relative error criterion used 





to determine if 

convergence has been 





obtained while iterating to solve the 
nozzle admittance boundary equation. 

To obtain convergence, it is necessary 





that 






|(i)l - u)2|/|w2l <EPSXS 





where the subscripts 1 and 2 refer to 
successive iterations, unitless. 

FI 

— 

R 

S0LVW 

Derivatives of absolute value sauared 
of the difference between the upstream 
and downstream nozzle admittances with 





respect to Imag 
XI, X2, and X3. 

(lo) corresponding to 

F2 

— 

R 

S0LVW 


U 

F3 


R 

S0LVW 


f 

FANS 

-- 

R 

S0LVW 

Derivative of squared error in nozzle 
admittance boundary equation at ANS. 

FI 



R 

MAIN 

Interpolating factor used in the main 





program. 


FIl 

— 

R 

S0LVW 

Used to store successive values of Imag 
(FN) during the 2-dimensional secant 





method . 


FI2 

-- 

R 

S0LVW 



FI3 

— 

R 

S0LVW 



FI4 

— 

R 

S0LVW 

1 

r 

FN 

— 

C 

FZER0 

Difference between the upstream and 
downstream nozzle admittances. 

FNF 

— 

R 

FZER0 

Imaginary (FN) 


FNR 

— 

R 

FZER0 

Real (FN) 


FRl 


R 

S0LVW 

Used to store successive values of Real 



(FN) during the 

2:jdimensional secant 


method. 
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FORTRAN ENGINEERING 


VARIABLE VARIABLE 

NAME SYMBOL TYPE 

FR2 - R 

FR3 - R 

FR4 - R 

FREQ 2Tru)p R 

FREQMI — R 

FREQMX — R 

FREQT(IOO) — R 

FRQMAX — R 

FTSTl — R 

FTST2 — R 

GIFU — C 

G10X — C 

G2FU — C 

G20X — C 


CONTROL 

SECTION DESCRIPTION 

S0LVW Used to store successive values of Real 
(FN) during the 2-dimensional secant 
method. 

S0LVW 

S0LVW t 

S0LVE Frequency, i.e., real (w), Hz 

C0MN0Z Minimum frequency, used for generation 
of frequency table FREQT, Hz 

C0MN0Z Maximum frequency used for generation 
of frequency table FREQT, Hz 

C0MTAP Table of terms used for computation of 
downstream nozzle admittance. 

S0LVE Maximum frequency above which no solu- 
tions to the nozzle admittance boundary 
equation will be sought, Hz 

S0LVW Test function used to determine if a 
solution to the nozzle admittance 
boundary equation has been bracketed. 

S0LVW Test function used to determine if a 
solution to the nozzle admittance 
boundary equation has been bracketed. 

ADARND Coefficient of oscillatory pressure in 
fuel oscillatory vaporization 
expression. 

ADARND Coefficient of oscillatory pressure in 
oxidizer oscillating vaporization 
expression. 

ADARND Coefficient of oscillatory density in 
fuel oscillatory vaporization 
expression. 

ADARND Coefficient of oscillatory density in 
oxidizer oscillating vaporization 
expression. 
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FORTRAN 

VARIABLE 

NAME 

ENGINEERING 

VARIABLE 

SYMBOL 

TYPE 

CONTROL 

SECTION 

DESCRIPTION 

G3FU 

— 

C 

ADARND 

Coefficient of oscillatory mixture ratio 
in fuel oscillatory vaporization 
expression. 

G30X 

“ “ 

C 

ADARND 

Coefficient of oscillatory mixture ratio 
in oxidizer oscillating vaporization 
expression. 

G4FU 

— 

C 

ADARND 

Coefficient of oscillatory velocity in 
fuel oscillatory vaporization 
expression. 

G40X 

- “ 

c 

ADARND 

Coefficient of oscillatory velocity in 
oxidizer oscillating vaporization 
expression. 

6AM0 

^0 

R 

C0MCBM 

Specific heat ratio evaluated at the 
overall mixture ratio, unitless. 

GINJFT(IOO) Ginjfu 

C 

C0MTAP 

Fuel feed system response table. Real 
(6INJFT) is the amplitude of the re- 
sponse and Imag (GINJFT) is the phase 
angle of the response. 

GINJOT(IOO) 

r 

C0MTAP 

Oxidizer feed system response table. 
Real (GINJ0T) is the amplitude of the 
response and Imag (GINJ0T) is the phase 
angle of the response. 

GN 

— 

C 

FZER0 

Variable used to store the value of the 
derivative of CN0ZA-N0ZA with respect 
to 0)1 or the value of CN0ZA-N0ZA itself. 

HN 

— 

R 

FZER0 

fn*Fn 

I 

— 

I 

- 

Used throughout the program as a do 
loop index. 

INPN0Z 


I 

C0MN0Z 

Code used for the downstream nozzle 
admittance term calculation. 

INRT 

— 

I 

DUMP 

Code used to determine whether or not 
intermediate output from CHAMDY is 
desired. 

I PASS 

— 

I 

S0LVE 

Internal code no longer in use. 
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FORTRAN 

VARIABLE 

ENGINEERING 

VARIABLE 


CONTROL 


NAME 

SYMBOL 

TYPE 

SECTION 

DESCRIPTION 

IPRN0Z 

— 

R 

C0MN0Z 

Code used for downstream nozzle admit- 
tance term printout 

IR 


I 

MAIN 

Flag set by MAIN program to indicate 
the first pass through it. After read- 
ing in a new case. 

ISCNT 

— 

I 

FZER0 

Code used to determine logical flow in 
subroutine S0LVW. 

ISLP 


I 

FZER0 

Code used to determine whether or not 
the derivative of FN with respect to 
Imag (u)) is needed. 

ISTRT 


I 

S0LVE 

Code used to indicate the first itera- 
tion after a solution to the nozzle 
boundary equation. 

II 


C 

ADARND 

The imaginary number i. 

ITAPC 

— 

I 

C0MTAP 

File number used to save combustion co- 
efficients table. 

ITAPH 

— 

I 

C0MTAP 

File number used to save hydrodynamic 
coupling term table. 

ITAPN 

— 

I 

C0MTAP 

File number used to save nozzle admit- 
tance term table. 

IWSKP 

-- 

I 

S0LVE 

Intermediate output dump code used in 
subroutine S0LVW. 

J 

— 

I 

- 

Used throughout the program as a do 
loop index. 

KNTMX 


I 

S0LVE 

Maximum number of iterations allowed to 
minimize the error in the nozzle admit- 
tance boundary equation with respect to 
the imaginary part of o). 

KNTR 


I 

S0LVE 

Counter used to control the number of 
iterations used during the search 
algorithm between solutions. 

KNTRMX 


I 

S0LVE 

Maximum number of times the frequency 
will be allowed to be changed by DELFRQ 
during the searching algorithm between 
each solution. 
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VARIABLE 

NAME 

ENGINEERING 

VARIABLE 

SYMBOL 

TYPE 

CONTROL 

SECTION 

DESCRIPTION 

KNTSMX 

— 

I 

S0LVE 

Maximum number of iterations allowed for 
the convergence of the two-dimensional 
secant method used in S0LVW. 

KSCNT4 

— 

I 

S0LVE 

Counter used to control the number of 
iterations used when ISCNT = 4. 

KWHERE 

— 

I 

MAIN 

Flag to control logical flow in the 
MAIN program after a call to subroutine 
S0LVW. 

MBFUI 


R 

C0MCBM 

Fuel injection mass flowrate, 
kg/sec (lb/ sec). 

MB0X1 

Kx'inj 

R 

C0MCBM 

Oxidizer injection mass flowrate, kg/sec 
Ob/sec). 

MG I 

n j 

R 

C0NSTS 

Steady-state gas flowrate at injector 
face, kg/sec. 

MR (100) 

MR" 

R 

C0MCHM 

Oscillatory mixture ratio, dimensionless 

MRB(IOO) 

HR 

R 

C0NSTS 

Steady-state mixture ratio, unitless. 

MRGI 


R 

C0NSTS 

Steady-state gas mixture ratio at 
injector face, unitless. 

MRNTFU 

-- 

C 

ADARND 

Mixture ratio integral in fuel 
oscillatory vaporization expression. 

MRNTPX 

-- 

C 

ADARND 

Mixture ratio integral in oxidizer 
oscillatory vaporization expression. 

MWG 

MW, 

0 

R 

C0MCBM 

Molecular weight of the gas evaluated 
at the overall mixture ratio, kg/kmole 
(Ibm/lb-mole) . 

NFREQT 

— 

I 

C0MTAP 

Number of points in frequency table. 

N0ZA 

ANj, 

C 

FZER0 

Downstream nozzle admittance 

N0ZAMR 

^f^MR=constant 

c 

FZER0 

Term used in computation of downstream 
nozzle admittance 

N0ZAT(1OO) — 

c 

C0MTAP 

Table of terms used for computation of 
downstream nozzle admittance. 
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VARIABLE 

NAME 

ENGINEERING 

VARIABLE 

SYMBOL 

TYPE 

CONTROL 

SECTION 

DESCRIPTION 

NR00T 

-- 

I 

MAIN 

Number of solutions to the downstream 
nozzle admittance boundary equation 
being sought. 

NRT 


I 

MAIN 

Used as do loop index for MAIN program. 
Counts the number of solutions to the 
nozzle admittance boundary equation. 

NUBFU 

^fu 

R 

C0MCBM 

Average steady-state fuel Nusselt number 
used in vaporization expression, 
unitless. 

NUB0X 


R 

C0MCBM 

Steady-state oxidizer Nusselt number 
used in vaporization expression, 
unitless. 

NXP 

— 

I 

C0MARE 

Number of points in axial distance array, 
inclusion between XO and the start of 
nozzle inlet. 

NXPMl 

-- 

I 

CHAMDY 

NXP-1 

0MEGA 

u> 

C 

C0MCHM 

Complex frequency. 

P(IOO) 

p" 

R 

C0MCHM 

Oscillatory pressure, dimensionless. 

PI 

TT 

R 

- 

pi (3.141593) 

PC 

P 

R 

C0MCBM 

2 

Steady-state chamber pressure, N/in. 
(psia). 

PINTFU 

— 

C 

ADARND 

Pressure integral in fuel oscillatory 
vaporization expression. 

PINT0X 

— 

C 

ADARND 

Pressure integral in oxidizer oscillatory 
vaporization expression. 

RBS0X 

— 

C 

ADARND 

Collection of terms used in oxidizer 
oscillatory vaporization expression. 

RCCX 

• ■“ 

R 

C0MN0Z 

Ratio of the radius of curvature at the 
nozzle inlet to the chamber radius at 
nozzle inlet, unitless (see Fig. 15). 

RCTX 


R 

C0MN0Z 

Ratio of the radius of curvature upstream 
of the throat to the chamber radius at 
nozzle inlet, unitless (see Fig. 15 ). 
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VARIABLE VARIABLE 

NAME SYMBOL 

TYPE 

CONTROL 

SECTION 

RG0 % 

R 

C0MCBM 

RHNTFU 

C 

ADARND 

RHNT0X 

C 

ADARND 

RH0(1OO) P' 

R 

C0MCHM 

RH0B(1OO) P 

R 

C0NSTS 

RH0GI 

R 

C0NSTS 

RH0INJ P'inj 

R 

ADARND 

RINJ 

R 

C0MN0Z 

RPSFU 

C 

ADARND 

ssvi(ioo) -- 

R 

C0NSTS 

SSV2(100) " 

R 

C0NSTS 

SSV3(100) — 

R 

C0NSTS 

SSV4(100) — 

R 

C0NSTS 

SSV5(100) " 

R 

C0NSTS 

SSV6(100) — 

R 

C0NSTS 

SSV7(100) " 

R 

C0NSTS 

SSV8(100) — 

R 

C0NSTS 

SSV9FU(100) — 

R 

C0NSTS 

SSV90X(1OO) — 

R 

C0NSTS 


DESCRIPTION 

Gas constant evaluated at the overall 
mixture ratio, J/kmole/'’K (Btu/lb 
mole/°R). 

Density integral in fuel oscillatory 
vaporization expression. 

Density integral in oxidizer oscillatory 
vaporization expression. 

Oscillatory density, dimensionless. 

Steady-state density, kg/m . 

Steady-state gas density at injector 
face, kg/m3 

Oscillatory density at injector face, 
unitless. 

Radius of the chamber at the injector, 
in. 

Collection of terms used in fuel oscil- 
latory vaporization expression. 

Steady-state parameters computed in 
subroutine STEADY for use by subroutine 
CHAMDY I 
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VARIABLE VARIABLE CONTROL 


NAME SYMBOL 

TYPE 

SECTION 

SSVIO(IOO) 

R 

C0NSTS 

SSVl 1(100) 

R 

C0NSTS 

SSV12(100) 

R 

C0NSTS 

SSV13(100) 

R 

C0NSTS 

SSV14(100) 

R 

C0NSTS 

SSV15(10) 

R 

C0NSTS 

T(IOO) T' 

R 

C0MCHM 

TAUBFU T, 

fu 

R 

C0MCBM 

TAUB0X T 

ox 

R 

C0MCBM 

TB(IOO) T 

R 

C0NSTS 


R 

C0MCBM 


R 

C0MCBM 

TITLE(18,2) __ 

R 

MAIN 

V(IOO) v' 

R 

C0MCHM 

VAPBFU(IOO) 

R 

C0NSTS 

VAPB0X(1OO) mvap^^ 

R 

C0NSTS 

VB(IOO) V 

VBFU V 

Jfu 

R 

R 

C0NSTS 

C0MCBM 

VB0X 

R 

C0MCBM 


R 

c?NSTS 


DESCRIPTION 

Steady-state parameters computed in 
subroutine STEADY for use by subroutine 
CHAMDY. I 


I 

Oscillatory temperature, dimensionless. 

Fuel vaporization time delay, sec. 

Oxidizer vaporization time delay, sec. 

Steady-state temperature, °K. 

Steady-state oxidizer drag time delay, 
sec. 

Steady-state oxidizer drag time delay, 
sec. 

Array containing 2 card records of 
alpha-numeric title information. 

Oscillatory velocity, dimensionless. 

Steady-state fuel vaporization rate, 
kg/sec/m*’. 

Steady-state oxidizer vaporization rate, 
kg/sec/m3 

Steady-state velocity, m/sec. 

Steady-state fuel liquid injection 
velocity, m/sec (ft/sec). 

Steady-state oxidizer liquid injection 
velocity, m/sec (ft/ sec). 

Steady-state gas velocity at injector 
face, m/sec. 
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VARIABLE 

NAME 

ENGINEERING 

VARIABLE 

SYMBOL 

TYPE 

CONTROL 

SECTION 

DESCRIPTION 

VINTFU 

-- 

C 

ADARNO 

Velocity integral in oxidizer oscil- 
latory vaporization expression. 

VINT0X 

-- 

C 

ADARND 

Velocity integral in oxidizer oscil- 
latory vaporization expression. 

VX0 

'-'x.o 

c 

C0MCHM 

Oscillatory velocity at X=0, 
dimensionless. 

X(IOO) 

X 

R 

C0MARE 

Axial distance array, 

X(I) = XO + (X-D* DELX, m. 

XO 

Xo 

R 

C0MARE 

Start of vaporization point, m 

XI 


R 

S0LVW 

Used to store successive values of Imag 
(w) during the iteration to minimize 
the error in the nozzle admittance 
boundary equation with respect to Imag 
(u) . 1 


X2 

— 

R 

S0LVW 

X3 

— 

R 

S0LVW 

XIMPFU 

ximpfu 

R 

C0MCBM 

XIMP0X 

xiupox 

R 

C0MCBM 

XII 


R 

S0LVW 


XI2 

-- 

R 

S0LVW 

XI3 

-- 

R 

S0LVW 

XI4 

— 

R 

S0LVW 

XKFU 


R 

C0MCBM 

XK0X 

“"ox 

R 

C0MCBM 


t 

Fuel jet injection impingement point, 
m (in.). 

Oxidizer jet injection impingement 
point, m(in.). 

Used to store successive values of Imag 
(cij) during the 2-dimensional secant 
method . i 


T 

Fuel Klystron distance, m(in.). 
Oxidizer Klystron distance, m(in.). 
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VARIABLE 

NAME 

ENGINEERING 

VARIABLE 

SYMBOL 

TYPE 

CONTROL 

SECTION 

XM(IOO) 

— 

R 

C0MARE 

XN0Z 

— 

R 

C0MARE 

XRl 

— 

R 

S0LVW 

XR2 

— 

R 

S0LVW 

XR3 

-- 

R 

S0LVW 

XR4 

— 

R 

S0LVW 


DESCRIPTION 

Axial distance midpoints, 

XM(I) = (X(I) + X(I-l))/2, m. 

Nozzle inlet point, m. 

Used to store successive values of Real 
(u)) during the 2-dimensiona1 secant 
method. I 


T 
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PROGRAM FLOW CHARTS 


FLOW CHART INDEX 

PAGES* 

Cross Reference Index 

1-10 

MAIN (Procedures) 

01-09 

AREA 

10-11 

CHAMDY 

12-14 

CHMC0N 

15-16 

COGAEL 

17-19 

C0MBDY 

20-22 

C0MHAT 

23-24 

FRESP 

25-38 

HEAD 

39-40 

HY (B10ck Data) 

41-42 

HYDRDY 

43-49 

L0CFAC 

50-52 

N0ZADM 

53-56 

RKTDIF 

57-58 

RKTZ 

59-60 

RKZDIF 

61-62 

S0LVW 

63-70 

STEADY 

71-72 

TADAMS 

73-76 

TDPL0T 

77-79 

XIMAGF 

80-81 

ZADAMS 

82-85 

ZER0 

86-89 


♦Flow chart page numbers in upper right-hand corner. 
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D6/a/T5 TABLE OT CONTENTS AFC fCTEfCNCES 

CABO JO PAGt/BO* NArt 


AUTOTLOH CHI«T SET - rsCS« 

WTCfipCES (SOURCE SEQUENCE NO. AFC PAGE/BOxi 


CHART TJTLC - subroutine PRESPt (CRT . |k«| TE . IX I 


00006010 
ooooeaio 
00006^20 
00006230 
00006230 
OOOC6260 
00006270 
00006270 
00006290 
00006300 
00006300 
0000631 0 
00006320 
00006330 
00006:^0 
00006380 
0000639C 
OOOC6NOO 
00007300 
000073NO 
000073UO 
00007160 
00007370 
00006H30 
00006H60 
00006'«50 
00006H60 
00006U90 
00006610 
00006520 
00006550 
00006560 
00006680 
00006690 
00006750 
00006760 
00006580 
00006610 
00006630 
00006650 
00006790 
00006820 
00006850 
00006860 
0ODO6S8O 
00006890 
00006920 
00006930 
000069*0 
00006990 
00007010 
00007030 
00007050 
00007070 
00007090 
00007110 
00007130 
00007150 
00007160 
00007190 
00007220 
00007230 
00007250 
— «:00007270 
00007280 
OOOQ7UOO 

ooooTmo 

00 007*. 20 
00 007*.*. 0 


25.01 TRESP 

25.01 10 

25.02 

25.03 2001 

25.03 

25.08 

25.09 2002 

25.09 

26.03 

26. ON 2003 
26. ON 

26.07 2006 
26.00 

26.09 

26.10 200N 

26.17 

26.18 1216 

26.19 12:8 

26.20 Nnw 
26.23 NN5 
26.23 

26.27 

26.28 n5 

27.01 1217 
27. ON 1215 
27. ON 

27.05 1 220 

27.06 1219 
27.09 

27. iC 1222 
27. n 1221 

27. 12 3 

27.13 1671 

27. IN IS75 
27, 15 N 
27. IE 5 
28.0! 1667 

28.02 1668 
28 . ON 1669 

28.05 1670 

29.01 8 

29.03 12 

29.05 23 

29.06 2N 

29.08 27 

29.09 25 

29.10 26 

29.12 II 

29.13 13 
29. IN 15 

29.15 16 

29.17 18 

29.18 17 

29.20 20 
29.22 30 

29. 2N 29 
29.25 35 

29.27 20 

29.28 33 

30.01 3N 

30.03 37 
30 . ON 31 
30.05 36 

30.07 38 

30.08 330 

30.10 32 
30. n N3l 
30.12 N32 

30.15 N3 


00005300 N0 . OB-X 

00006230 25 . 05 

00006230 25. ON 

00006270 26.01 

00006270 25.10 
C0006300 26.06 

00006300 26 . 05 

OOC062NO 25,06 
OOC063NO 26. 12 
OOCQ63NO 26.11 

00006360 26. IN 


00006350 26.13 

00007y*0 26. PN 
00007350 26.25 

00007300 26.20 00007100 29,23 0000729C 30.09 

0000639C 26.18 

OOOC638C 26.17 
00006N30 27.01 

00006N50 27. ON 
00006N90 27.06 

00006510 27.09 


00006550 27.1! 

00006780 27.16 

00006550 27. M 
OOOQ66NO 26. ON 

00006620 28.03 
00006750 27,15 
00006930 29.12 

000068NO 29. ON 

00006870 29.07 
00006870 29.07 
00006920 29.10 

00006930 29.12 

00006920 29.10 

00006670 29.05 

00007100 29.23 

00007060 29.19 
00007100 29.29 


OOOOG9BO 29.13 


00006880 29.06 


00007260 30.06 


00007020 29.16 


00007270 30 . 07 


000071NO 29.26 00007150 29.27 






000072NO 30. ON 

00007220 30.03 

00C07360 26.27 


00007080 29,2! 00007270 30.07 0000729C 30.09 


OOOOTNIO 50.11 00007510 30.17 


RACE 


3 


R-9808/B-4 



06/25/T5 T4BLX OT CONftNTS AhO WTERCNCES AUTOTlOM CHART SCT - rSCSM 

CARO 10 PACC/BOA NATC RCrCRCNCCS i SOLACE SCQUCMCE NO. AND PACE /BOX) 


0000‘THSO 

WOOTHflO 

OOOOTSSO 

00007500 

00007500 

000075NO 

00007560 

00007570 

00007580 

00007600 

00007630 

00007670 

00007680 

00007690 

OOOOTNO 

00007750 

00007770 

00007780 

00007780 

00007780 

00007800 

00007830 

00007830 

00007830 

00008010 

ooooeoBo 

00008110 

ooooeteo 

00007850 

00007860 

00007860 

00007870 

00007900 

00007^0 

00007950 
00007970 
0000»*60 
00008010 
00008030 
00007990 
00008050 
00008300 
00008360 
00009^30 
00009*90 
C0008500 
00008510 
00008500 
00008550 
00008560 
00009590 
00008610 
00008650 
OOOO068G 
00000700 
0O0087H0 
00008770 
00008780 
C0008790 
CO0C88OO 
00008800 
00008910 
00008870 
00008900 
0000891 0 
00008950 
00000960 
00a09l’*0 
00009190 
00009000 
00009060 
00009080 


30. IH 39 

30.15 **0 

30.16 NJ 

30.17 N0 

31.01 N** 
31.0^ 

31 .06 

31.07 301 

31.08 <*8 

31.10 300 

31.11 300 

3!. 13 M7 
31.1** 49 

3t.l5 60 
3t . 18 

31 .00 
31 .03 
31 04 111 

31 04 
31 07 

31 .30 

31 31 311 

31 .3! 

31 .34 

30.01 69 

30.03 55 
30.05 56 
30.00 57 

30.09 309 
30. 10 75 
30. 10 

30.13 

30. 76 

30. 16 60 

30.17 65 

30.10 67 

30.00 7 

30.01 64 

30.03 71 

30.04 66 

30.05 70 

30.07 770 

30.09 773 

30.30 63 

33.01 610 

33.00 61 

33-03 300 

33.04 660 

33.07 

33.08 

33.10 3l6 

33.11 310 

33.13 313 

33 17 
33.19 70 

33.01 73 

33.04 303 

33.06 
33 07 

33.08 304 

33.31 

33.30 3C5 

33.35 

33.37 

33.38 000 
33.40 

33.43 003 

34.04 

34.06 

34.08 006 
34.10 005 

34.13 80 


00007440 30.13 

00007470 30.14 

00007490 30.16 00007680 31.14 

00007500 31.01 

00007670 31 13 


00007560 31.06 

00007570 31 .07 

00007540 31.04 

00007710 31 . 16 

OOOOTNO 31.18 
00007750 31.00 

00007780 31 07 

00007780 31,04 

00007800 31 .08 


00007610 31.09 

00007610 31.09 


00007600 31.10 

00007660 31.10 


00007830 31.34 

00007830 31.31 

00007980 30. 19 

ooooeno 30.01 

Q0008ICO 30.04 

00007040 31.35 

00007850 30.09 

OOOO850C 33.04 
00008510 33.03 


00009690 35.00 


00008000 30.04 

00308100 30.34 00008170 30.07 00007930 30.15 00008350 30.06 00009*00 30.31 

00007940 30.16 00008040 30.03 


00008000 30.00 

00008400 30.31 

00008350 30 . 08 

03007930 30.15 

00008000 30.08 00008450 30.30 


0000865C 33.14 

00000650 33.13 

03008580 33.09 

00008500 33.09 

00008740 33.01 

00008690 33.18 

00008770 33.04 

ooooeeoo 33.31 
ooQoeeoo 33. ae 

00008760 33.03 

00008910 33.39 
00008910 33.38 

00008980 33.43 

00009060 34 , II 

00009000 34.08 


PAGE 


R-9808/B-5 



PAOC 


06/a5('75 TAB.C OT CONTENTS *fC RErERENCES AWTOTLOM CHART SET - FSCSH 

CARD 10 PAGE/BOX MATC REFERENCES ISCW»CE SCOUENCC NO anO PaGE/BOX' 


00009J50 y».l7 T70 

000095M >,I7 

0O0C»*?0 >*.t9 771 

00009>*e0 3^. 19 

oooo9«5o 

oooo»*6o j^.ee 79 

0000»«70 3V.B3 BC 

00009^60 35.01 91 

0000».90 35.03 

000095CO 35.05 

000Q95I0 35.06 7B 

00009530 35.06 83 

000095SC 35.09 

00009570 35.10 ^ 

00009560 35.11 86 

00009590 35. :a 90 

00009590 35. I a 

000096CO 35.13 96 

00009610 35. 87 

000C9630 35.16 960 

00009670 35.19 961 

00009670 35 19 

00009660 35.80 65 

00009700 35.81 95 

00009770 35.83 97 

00009780 35.8^ 90 

00009790 35.87 

00009600 35.88 

00009680 35.31 99 

00009880 35.31 

0000987C 36,01 I 00 

0000991C 36.03 

00009980 36 . Os 

00009970 36.08 


00009890 > . ts 

00009350 3S.17 

00009510 35.06 

00009S50 3S.81 

00009S50 3S.8! 

0000^80 35.0! 

00009S90 35.03 
0D009S7O 3s. 83 

00009870 3S.18 

00009570 35.10 


00009610 35. IS 


00009680 35. 15 

00009630 35.16 

00009670 35.19 

03009660 35 , 38 

0000978O 35 , 8S 

30009680 35.3! 

00009800 35,88 
00009770 35 83 
00010880 36.35 

00009990 36.10 
00009950 36.06 


00009990 

00009990 

00010010 

00010080 

OOOIOOSC 

000)0180 

00010190 

00010880 

00010830 

0OC108S0 


36.10 108 

36.10 
36. 16 
36.10 

36.80 103 

36.85 10S 

36.38 105 

36.35 soo 

36.36 999 

37.01 9 


000O99S0 36.05 
00009970 36 . 08 
OOOlOOOO 36. IS 
OOOlOOtO 36. 16 

000:0030 38. 19 
OOOlOilO 38.13 
00010080 38.16 

OOD06S80 86. 19 

ooooeseo 38 . 8o 


00007390 86.88 


OOOO6S0C 87.05 


000065SO 87. !0 


00 006570 87.: 8 


CHART TITLE - NON-PROCEOURAL STA^trCNTS 


CHART title - SUBROUTINE ^CAD 

00000030 39,0! (CAO 00000710 8.08-X 

CHART title ' mon-procedural STaTETCNTS 


CHART title - BLOCK DATA 


CHART title - NON-PROCCOURAL STATETCNTS 


CHART TITLE - SUBROUTlfC HVORDYt IR. lT*4frO.FRE ,GtNO,GM^.PClN.UO’N.HriNI 


OOOOOOtO S3. 01 HYOROV 

00001180 S3. 05 

00001170 S3, 06 5 

0000 It BC S3. 07 

00001190 S3.0B to 

00001190 S3. 06 

00001810 S3. 18 

00001880 S3. 13 80 

00001880 S3. 1 3 

OOOOI8SO S3. 17 30 

0000 1 8S0 S3. 17 


000O8S0O s . is-x 
0000! in S3. 03 
00005309 SB. IS 
00001190 S3. 10 

OOOOU90 S3.M 
00001880 S3. 15 

00001880 S3.1S 

0000 laso S3. IB 


R-9806/B-6 



06/2V75 T«L£ OT COMTENTS AND REFERENCES 

CJMD ID PAGereax nakc 


AUTOn.OW CHART SET - FSCSH 
REFERENCES (SOURCE SEOUDCE NO. AND PACC/BOXI 


OOOOtHtO 

OOOOIHNQ 

0000t<«90 

00001600 

00001610 

00001620 

OOOOI6SO 

0000I6B0 

00001660 

00001680 

ooooiees 

ooooiees 

00009430 

00001890 

oooonio 

00001710 

ODOOKSO 

0000S300 


00005306 

00005308 

00005303 

0000530>« 

00005310 

00005940 

00009400 


>«3>2e 35 
>«3.25 
N3.26 HO 
HH.OH H5 
HH.07 


HH. 1 1 

HH.IH 

HH.18 

HN.IS 100 

HH.18 

NH.19 999 

HH.21 

*44.22 101 

*44.22 
48.06 HOC 
48.08 500 
48.11 502 

48.13 503 
48.13 

46.15 

48.16 501 
48.18 510 
48.22 
48.26 700 


OOQOIH20 43.23 
OOOOtBOO 44.04 
00001610 44.07 

00001600 44.06 

00001620 44.09 

00001650 44.11 

00001660 44.14 

00001410 43.22 
00001680 44.16 

00001650 44.13 
00001685 44,18 

00001710 44.23 

00005280 46.06 
0000530: 46.09 

00005306 48.14 

C0005304 48.17 

00005306 48.lt 
00005320 46.19 
00005310 48.18 


00005305 48.17 


CHART TITLE - NCM-PROCEDURAL STATEftNTS 


CHART TITLE - SUBROUTirC LOCFAC UK . X . TX.NX , JX .FX) 


00000180 

00000610 

0000 02 70 

00000280 

00000290 

00000?«0 

00000250 

00000320 

00000940 

00009*00 

00000430 

00009*50 

00009*60 

00009*70 

00009*90 

00000500 

00000510 


50.01 LOCFAC 

50.06 

50.10 200 

50. 11 30 

50.12 

50.13 40 

50.15 

50.16 20 

50.17 50 

50. 18 70 

51.01 90 
51 .9* 

51 .06 100 

51.08 110 

51 .09 

51.11 120 

51.12 130 

51 . 14 
51.17 150 


00002150 

oaoociTo 

00000150 

00000220 

00000290 


3.23-X 
50.04 
50.02 
50.09 
50 . ; 3 


00002320 4.09-X 


00000440 51.05 00009*80 51.10 00000540 5l 16 


00000250 50.16 


00000280 50.12 00000240 50.15 


00000590 51.; 7 
00000190 50.07 
00000450 51 . 07 


00000300 60.13 00000260 50.16 


00000410 51.02 00000440 51.05 
0000049C 51. n 

00000480 51.10 
00000500 51.12 
00000520 51 . 15 


00000530 51.16 


CHART TITLE - NW-WOCEDURAL STATErCNTS 


CHART TITLE - SUBROUTINE N 02 AOMi|R.CA«t.C 0 ,FREO,NO 2 AI 


00000030 
00000290 
00009*10 
00009*70 
000009*0 
00000570 
00000620 
00000670 
00000700 
00000710 
00000710 
^ 00000720 
DOOQIIIO 
00001220 
00001440 
00001490 
000019*0 


53.01 NQ2A0H 

53.03 
53.08 2 
53.12 

53. 15 2000 

53.16 2010 
53.20 5 
53.23 7 

53.27 

53.28 20 

53.28 

53.29 25 
9*. 09 
9*. 12 30 

95.01 ^ 

95.03 40 
96.06 45 


00002230 3.25-X 

00000270 53.0! 
00000290 53.03 
00000460 53. 09 
00000490 53.13 
00000530 53.14 
OOOOOG60 53.22 
00000600 53.18 
00000680 53.24 

00001810 55.24 

00001220 54.13 

00001390 54.19 
00001430 54.20 


R-9806/8-7 



06/a/T5 table or contents AM> WETEW E NCES 

CARO ID PACC/80X NM« 


00001600 

00001690 

00001730 

00001730 

ooQoieio 

0000 I 610 
00001900 
00001910 


95.1V SO 

95.16 e 

99.16 
99. 2s 
95.2V 
95.27 
95.29 


10 


00001660 99.10 
00001530 55.05 
00001680 55.13 
00001710 55.15 
00001290 5V.IB 
OOOOIBOO 55.21 
00001090 55.25 
00001900 55.27 


AUTOFLOW CHART SET - FSCSrt 
REFCREICCS <SOURCE SEQUENCE NO, A^O paGE/BOX) 


OOOOI7VO 95.19 


CHART TITLE - NON^>ROCCCMUL STATETCNTS 


CHART TITLE - SUBRCXITIAC RKTDIFIP.C.CP) 


00002650 

00002620 

00002630 

00002920 

00002930 

00002950 


57.01 WTDIF 
57.0V 

57.05 25 

57.06 50 

57.07 95 

57.12 30 

57.13 35 
57.1V VO 
57.15 V5 


00003750 59.07-X 00003860 S9.20-X 00003950 59.31-] 


00002610 57.03 
0000291 D 57. 15 
00002620 57.0V 


00002650 57,12 
000029«0 57.05 


00002970 57.13 


CHART TITLE - NON-PROCEOLRAL STATEfCNTS 


CHART TITLE - SUeROUTlrt RKTZINU.M. T 1 .U.DUM. JOPT I 


00003560 
00003720 
00003730 
00003770 
00003790 
00003790 
00003790 
00003810 
00003830 
000038V 0 
00003800 
00003090 
0O0C390O 
00003900 
00003910 
00003930 
00003930 
00003950 
00003970 


59.01 RKTZ 
59.03 
99.0V 10 

99.09 15 

99.10 20 

99.lt 25 

59.11 
59. IV 
59. 16 

59. 17 35 

59.22 VO 

59.23 VS 

59.2V 50 

59.2V 
59.26 30 

59.28 55 
59.26 
59.31 60 

59.33 65 
59.3^ 70 

99,39 75 


oooomo 5 N. 09 -X 00001 no 75 07 -x ooooavic ev.oi-x 

00003730 59.05 

000037VD 59.06 

0000 376D 59.08 

C0003790 59. 12 
00003910 56.26 

000O38V0 59.18 

00003850 59. 19 

00003870 59.21 

00003900 59,25 


00003930 59.29 
000039^0 59.30 
OOD039VO 59.30 
00003960 59.32 
00003980 59.3V 


CHART TITLE - NON-PROCEOURAL STATEfCNTS 


CHART TITLE - SU8R0UTII« RKZOrFIP.G.GPt 


00003070 
00003210 
00003270 
00003320 
00003330 
000033«0 
00003350 
00003VVO 
0000 3V50 
0000 35VD 
00003370 
00003360 
0000 3V00 
00003V 1 0 


61.01 RKZ01F 
61.03 10 
61.0V 15 

61 .07 
61 .09 
61.10 22 

61.11 25 

61.12 50 

61.13 55 
61.16 20 

61.17 30 

61.18 35 

61.19 VO 
61 20 VS 


00003770 59.09-X 00003880 99.22-X 00003970 99.33-X 

00003200 61.02 
00003310 61 .05 
00003320 61 . 07 


00003330 61.09 

000C>»30 61.20 

00003260 61,03 
00003y«0 61 . 10 


00003370 61.17 

00003360 61.11 00003390 61.18 


CHART TITLE - non-procedural STATEHENTS 


CHART TITLE - 9U8R0UTI^C Sa.VW<W#CRE) 


R-9808/B-8 



06/»/ri table of contents and retepences 

CAAO 10 PACC/BOX NA^C 


AUTOFLOW CHART SET - FSCSN 
REFERCNCES (SOURCE SEOUCNCE NO. />iO PAGE/BOX> 


00000030 
000039*0 
000009*0 
D00009I0 
000007BO 
00000600 
OOOOOE70 
OOOOI3IO 
0O0CIO1O 
00001030 
0000 1050 
00001 IHO 
0000 laso 

00001530 

00001760 

00001780 

000019*0 

00001990 

Doooei'io 

ooooaieo 

00O0T780 

0000^360 

00009*70 

00005510 

00005560 

00005550 

000059*0 

00003000 

00003050 

000056CO 

0000579C 

00005930 

0000395C 

00003170 

00003550 

00003560 

0000357C 

0000365C 

00003670 

00003730 

00003880 

0000390C 

0000396C 

00003970 


63.01 SOLVW 

63.06 555 

9*. 01 117 

9*.o^ ne 

6H. 10 

6H.l^ no 

65.15 115 

6M.19 f*5 

65.01 150 

65.05 155 

65.03 135 

65.07 137 

65.10 ISO 

65.15 !S6 

66.01 

66. DM 

66. 13 IS7 

66.15 150 

66. 19 160 

66.51 170 

66.55 505 

66.57 1 BO 

66.59 1B: 

66.30 185 

66.33 1855 

67.01 175 

67.05 IT? 

67.03 !85 

67.9* 190 

67.10 183 

67 .19 104 

67 55 

67.56 550 

68.01 195 

68.05 
68.08 
68.17 500 

68.50 501 

60.51 505 
68.53 503 

69.01 10 

69 05 15 

69 03 5000 

69 . 04 6000 


00005980 5.36-X 

00000410 63.03 
00000780 64.10 
00001050 65.01 

0DO00750 64.07 

00000970 64 . 06 

00000750 64.17 

00000840 64.01 

00000670 64.15 

00001070 65.05 

C0000950 64.05 
00000550 64.15 
00001730 65.18 
00001760 66.01 

00005590 66 . 34 

00001800 66.06 
00001880 66.11 
00001810 66.07 

00001900 66.15 
00005140 66.19 

00005660 67 . 17 

00005740 67.09 

0DC05530 66.31 
00005180 66.53 
00005510 66.54 
00005640 67.15 
00003530 68.16 
00001790 66.05 

00005610 67.11 

00005890 67.55 

00005160 66.55 

00000460 63.04 

00003500 66.05 
00003550 68.05 
00003370 68. 14 

00003GOO 68 18 
00003610 66 19 
00005100 66.17 

00001470 64,55 

00001560 65,11 

DQ003950 63.06 


00000830 64.11 


00001050 €6.04 
00000850 64.05 


00001510 65.09 


00005950 67.55 


00003360 66.13 


00003660 68 SO 00003750 68.55 00003770 68.53 


00005650 67.16 


00005650 67.13 00003330 68.11 


00003670 66.51 
00005130 66.18 
00003160 67. Oe 
00003930 66.56 
00003930 67.56 


00005560 67.05 


00003890 69.01 0000391G 59.05 


CHART title - NON-PROCEDUfiAL STATEMENTS 


CHART title - SUBROUTINE 

00000030 71,01 STEADY 

00000440 71.03 

O00QO550 71. 07 15 

00000530 71 08 16 

00000590 71.10 

OOOOlllO 71.59 

00001540 71.34 30 


STEADY! IPRSTEl 

00005860 5.34-X 

OGOOG360 71.01 
00000470 71.04 

00000510 7; .06 

00CD1540 71.34 

OOODlOeO 71.57 
00001130 71.59 


CHART title - NON-PROCEOtRAL STaT£>CNTS 


CHART title - SUBROUTltC 


TAOAMSrN.H.X.Y.DY. IQZ.IO) 


00000030 
000001 10 
00000130 
00000150 
00000570 
00000580 
00001300 
00001350 
Q000030C 
00000310 
C0000330 


73.01 10 

73 03 
73 04 15 

73.08 17 

73.09 50 

73.10 100 

73.11 105 
73.15 55 
73.13 30 
73. 14 35 


00005500 84 13-X 
00000950 74 55 
00000150 73.04 


00000560 73.07 00000860 74.09 

00001590 75. 15 
00000570 73.08 

00000300 73.15 



R-9B08/B-9 



06/25/75 
CATO 10 

000003^0 

OOOOOHTO 

00000^90 

00000560 

00000570 

00000590 

00000630 

00000670 

00000700 

00000710 

DC000600 

00000730 

000007HO 

00000920 

00000930 

00001 no 

00001220 


table or CONTENTS AND RE F UC NCES 
PACE/BOX NArC 


73.15 HO 

73.21 

73.22 H5 

7H.02 52 
7H.03 
7H.0H 55 

7H.0B 60 
TH.IO 62 
7H . n 70 
7H, 12 75 
7S.13 65 
TH. I'* 90 
7H. 15 95 
7H.25 97 
75.01 90 

75.07 
75-10 95 


00000290 73.09 
OOOOQH90 73.23 


00000590 TV. 05 


00000670 TH.IO 


00000700 7H . 11 
00000720 7H.12 
00000660 7H.20 
00000620 7H . 07 
00001220 75.11 


AUTOTLOW CHART SET • fSCS« 

REFERENCES ISOUTOE SEQUENCE NO. AND PAGE/BOX* 

00000320 73.13 


OOOOQG90 7H . 1 3 


CHART TITLE - NON-PTOCCOJRAL STATEICNTS 


CHART title - SUeRO/TIIC TGPLOT^W.Y.^FP.TL.XR.LL^ 


00011002 
00011050 
00011070 
00011120 
0001 1150 
0001 1 190 
00011220 
00011260 
00011270 
00011280 
00011290 
00011310 
00011330 
00011350 
00011360 
00011380 
00011390 
0001 IHIO 
0001 1H20 

0001 mo 

0001 IH50 
O00I1H7O 
0001 1H80 


201 


77.01 
77.03 
77. OH 
77.07 
77.0! 

77.09 200 
77.11 10 

77.17 

77.19 

77.20 20 

77.21 30 
77.23 31 


77,25 

77.27 


76.01 

78.02 

78.04 

78.05 
78.07 

78.06 


50 


00010090 36.11-X 00010050 36.21-X 00010130 38,26-X D0010I70 36.30-X 

0001 1070 77. OH 


0001 mo 77.06 
DQOt I IHO 77,07 
0001 1270 77.19 
CO0U250 77.15 
00011260 77.17 

0001 12H0 77. IH 
00011230 77.12 


OOC113DO 77.22 
00011 3H0 77.26 
00011H30 78.03 

000113 
00011380 
0001 IHOO 77.32 
0001 1H30 70.03 


77.30 

77.30 


00011320 77. 2H 


0001 1H30 78.03 OCCIIHIQ 7803 


0001 mo 

coot 1H60 


78 . Oh 
78.06 


CHART TITLE - NON-PROCEOURAL STATE^tNTS 


CHART TITLE - FUNCTION XIHAGT (XIMOMGl 


OOOCCN60 80.05 10 

00000630 80.12 20 

O0OC0620 80. IH 

00000750 60.15 30 


00000370 80.03 
0000051 D 80.08 
00000610 80.10 
00DO062O 80. I H 


00000700 80.13 


CHART TITLE - NON-m>CEDURAL STATETCNTS 


CHART TITLE - SUBBOUTIhC 2A0 AMSin.H.X,Y.DY. 1QZ> 


00001360 

oooomo 

00001H60 

00001H60 

00001600 

00001610 

00001630 

oooomo 
00001660 
00001570 
00001790 
0000 1800 
00001870 


82.03 
82. OH 15 

82.08 17 

82.09 20 

82.10 25 

82.11 30 

83.01 35 

83.02 HO 
03.07 

03 08 H5 
83 1 1 52 


00001260 5H.15-X 
OCOO2210 83.33 
OC001H80 82. OH 


00001600 82.06 

0000 1630 92.10 
00001620 02.09 
OODOteOO 83.09 


00002590 8H , 1*4 


00001650 82.11 


R>980B/B-10 



m/am table or coktents mefdcnces 

CMC ID PMC/BOK HWC 


AUTOrUM OMIT VT - F9C91 
W UUtNah {SOURCE SEQUENCE NO. AfC PAOE/BOKI 


ooooiao 83.12 

OQOOISOO 83.13 SS 

00C01930 33.17 60 

00001970 83.19 6S 

00001980 83.20 85 

00002000 83.21 70 

00002Qtr 83.22 75 

00002030. 83.23 80 

00002C10 83. 2H 85 

00002210 83.33 87 

00O0T220 01*. 01 90 

ooor2mo v*.o7 

0QC02520 9t.l0 85 

00002610 0«.i5 tra 


oooeiaoo 83 . ih 


00001870 83.19 


83.21 
00001900 83.20 
00002150 83.28 
00001920 83.16 
00002520 84.11 


00002020 83.22 


00001980 82.07 00001860 83.18 


owrr Tin£ - NON-raoccuuRM. statetcnts 


PAoe ic 


CHART TITLE - SLBROLTrilC 2ER0(ri , Tl I ,T21 ,ri 1 .f2I .ANS.PANS .EPST.EPSX.AT. lOIT.MCN 


000009*0 

00000160 

00000200 


00000250 

00000270 

00000310 


88.01 

86.04 

86.05 
86.07 
87.01 
87. ON 

87.05 


00001190 e5.08>X 
000002HO 86.07 

00000130 86.02 
00000160 86. ON 
00000260 87.02 


00000710 87.3* 00000900 87.35 00000760 88.01 00000820 88.06 


00000960 87.12 500 


00000910 

00000930 

00000980 

00000550 

00000600 


00000720 

00000770 

00000830 


87.14 3DC 

87.15 NOO 
87.23 600 

87.25 230 

87.26 2ND 
87.28 2N5 
87.30 250 
87.35 270 

88. 01 255 

88.02 260 
80.07 265 


00000330 87.06 
00000390 07.09 
00000670 87.32 
00000500 07.19 
00000520 07.21 
000005NO 87 . 22 
OMOON80 87,17 
00000600 87.26 
00000620 87.28 
00000680 07 . 33 
00000610 07.27 
00000790 88.05 


00000780 88. ON 


CHART TITLE - NON-PROCEOLRAL STATE«>JTS 


R-9806/S-I1 



06/«/T5 


TA8LC or diagnostics 


autotlcm chart SC' 


rscSM 


RAoe t 


LOCATIOI OIAGNOSTJC 

CARD to PACE /BO* 

oooiotoo 36. le uACcn«:D - printv extomic rcftrcnce 

00010060 36. UrOCriTCO - PRINTV EXTTWML RErCREMCE 

00010070 36.23 UrCCflfCD - PRINTV' EXTEWAAL REFERENCE 

OOOIOIHO 36.27 UNDCFI(«D - PRINTV EXTCIMAL REFERENCE 

00010150 36.28 UNOCFIKCD - PRINTV' EXTCRMAt REFERENCE 

00010180 36.31 UAOEFlhCO - PRINTV EXTE»»IAL REFERENCE 

OOOlOlSO 36 32 UAOCrirvCO - PRlNTv EXTEf»«L REFERENCE 

DO0I02I0 36 3V UNCEFINCD - LA9LV EXTERNAL REFERENCE 

D0011310 77.23 UNOCrihCD - CA.“HAV EXTEWWL REFERENCE 

00011330 77 25 UNOEnhCO - "SETMtV EXTERNAL RE'TRENCE 

000U 350 77,27 UNOEF’tvtO - SETMlv ExTEWAL REFERENCE 

00011360 77,20 UNOEFISCD - Snxyv EXTERNAL REFERENCE 

COO11370 77, 2S U^OeFl^CD - GRlDlv EXTERNAL REFERENCE 

00011390 77 31 jrCCF I ^ED - ’APLOTv EXTERNAL REFERENCE 

OOOllNiO 79,01 UN0EF1»CD - APLOTV EXTCfMAL REFERENCE 

0001 1750 79 05 UNCCFIJCD • LASLV EXTEWAL REFERENCE 

00011^70 78.07 UNOEFirCD - 'LAeLV EXTERNAL REFERENCE 


R-9808/B-12 



06/25/7S 


«/rOFU3M 0«#rT SET - fSCSR OOWTER MWORVl 


PAGE Ot 


own title - iwnwxcTORT oomcnts 


miN CCKTICL pm»m for tic feo system copya stability moxl 
PROQRMI OCVELOPCD by nOCICTDYlC. A DIVISION €T NOOBCLL 
INTCmATIONAL. CAMOGA PAW. CALIF 9I3CN 
mOORATtCD BY n. D. 5CHLMAN. ROCKCTOYlC . HAY I9T?5 


ORIGINAL PAGE IS 
OF POOR QUALrrXi 


R-9S06/B-13 



06/»/T5 


AUTOFLOU CMWT SET - FSC»< CtmntR PMOOmn 


^4CC 00 


CH«T TITLE - PWOCECMCS 




R-9808/B>14 















0B/W/T5 


MmjFLOM CMKITT SET - rSCSM CCmjTER PWKRV1 


PAGE 03 


CHART TITU - moCEOURCS 



If 
















06 / 25/75 


M/TOTLOW OhART SET - FSCSw CO«»U’TEB PROGRAM 


PACE CN 


CHART title - PROCEDURES 


Z !2. 


T 

y 





NTERNAL rORKAT/ 
INTO ThC list / 


/ T 

ITAPH / 



LIS’’ • 
(GiNJOTr ! ) . 
CINJFTi [ ! .1 
1 .NFREO'i 



m 

Q 3 30 4 


n 

25 I 09 

I |3 

MYOROY [! 



tlR.l^PHVO. 

is j Ij 

‘6 

FRED. GINJOX, , 

lO ‘ LOCTAC 

P 

GINJFU.PCIN. I 

; (JK.FREO. ; 


MOIN.kFlNl I 

,0 I rREOT.tTREOT. i I 

i i 

-T— ^ 

; ; t ‘ n ,r; } ! 





CONTINUE ]| 


If" KGiN M Coop ' 

i | 56 : » 1 ■ NFREOT 


26 


GIT4J0TM1 - 
REAL I G : njOT 1 1 1 ) 
•CExP( lO.C. - 

.0) I 

•AIMAGCOlNJO’n ) ! I 

•Pi^iao.) ; 


I 

27 


; GtNJTTt : I . 

I REALfGlNjT'f M n 

I ’CEXPr ' 0 . 0 , - 

! . 3! 

I •aimad(g:nxtu n 


CtWJTt PEED SYSTEM 
RESPONSE PARAi^TERS, 
GINJOX AAO GIN^U. 
(T)«ouoH statement 
70 ) 


GINJOX - 
REALtClNjOXl 
•CEXPHO 0 . - 
t .0* 

•AIMAGlGtNJOX! 
*PI/1GO ) 


17 


CINJTJ • 

; realiginjtu' 

t •CEXPdO.O. - 

; t .0) 

•AIMACIOINJTLM 

•P!/ieo.i 


’'yes 


5.01 


5 . 02 ; 
^ 70 


R-9808/B'16 



06/25/75 


autqflou chart set - rsCSM 


COM^nCT PROORAM 


PACE 05 


CHART TITLE - PROCECUCS 















06 / 35/15 


AUTOFLOH OORT SET - F5CSM 


COWJTER program 


RAGE 06 


CHART TITLE - PROCEDURES 




R-9608/B-)8 













06/2S/19 


M/ran£M own srr - nc9i omncR pmoorw^ 


^M3C n 


owrr TITU - NQM-RWmjRAL statetcxts 


H9Q 

900 

500 

600 


5t0 

610 


620 

630 


DirOSION TITIXC18.2) 

coHux orcGA, p. no. v. m. t. cnoza. vxo, 

OOKI. C0QC2, 00K3, CWH. C0X5, COXB. CCM7, C0X8. C0X9. 

OQK10. COKII. C0K12. C0K13, COXIH, C0XI5, C0X16. 

crut. CRje. cruj. crw. crv«. cnje. cnr. crue. crus. 

CF\no. crun. cruiz. crviis. crui**. cfuib. crjis. 

NOZA. GtNJOK. GINJTU. 

NOZAT. OINJDT. GlurT. 

FN. NQ2AM). GN, 

COOCTtlS) . CFUTU6). 
coxTT(i6,2). crxrrrii6.a) 

REAL reOKI, reFXM. F«B. tMGI. MGI . NU80X. NUGRj. KT . KO. L. HC 
cortw /COMCBHy XKOX. MCFV). fOOXI , fCFUl, TAUBOX. TACBFU, V90X. 
VBTU. 0AI10. AGO. OEUCK. GEL>*\>, PC. CO. 

COKl. C0K2. COX3. OOMH. 000(5. C0X6. C0X7. C0X8, COX9. COXIO, 
C09UI, COX12. COK13. CO0(IH. 00X15. C(»16, CfUl . CFVg. 

(Tus. criA«. crcs. cr\j6. cnn, crue. crue. cfuio. cruii. 
cmi2. CFUI3, CTU 1 H. CFUI5. CTU16. X!fWU. XI^POX, 

CS. ocsom. ocm. ONGOm. aovox, aodox. toraco, oclvox, 

NU80X. OTOXCM. AOMTU. AOORJ. TORAGT. OELVTU, NUBFU. DTFUOM 
EGU1 VALENCE 1COXT (}> .000(1 > . (CTUTU ) .CTUll 
COOtO* /CONSTS/ tWBllOO). TBIIOO). «H08<100). VBHOO). 

□Henooi. □moeiiDO). oveciooi. vapsoxhooi. VAPeruiioo), 
ssvniooi. SSV2II0O1. SSV3U00). SSWUOO). 

SSV5MOO). SSV6MOOI. SSV7nOOI . SSVBHOOt. SSVSOXnOD). 

ssv9FU(iotn. ssvioiiooi. ssviniooi, ssviaiiooi. ssvisnoo). 
ssvmiioo). ssvtsdooi. ssvieuooi. 

FMX>I. VCl, m}\, NGl 

(XmON /COMOPV PMOO). I»OUOO). VCIOOI. wnoo). TU00>. 

VKO. OrCGA. CNOZA. OCU> 

Cm«N /COriARE/ WP. xn00>, XH(IOO). AClOOl, OAdOO). OELX, 

XO, XNOZ, AINJ 

COrtVM /CW*«2/ RCCX. RCTX. ANCLEX. CRR. RINJ, M^OZ, FRCO«. 
novii. iPfMoz 

cortw /FZERO/ MQZA.NOZAHO.QN.FN.F»«,FN1 .►♦I, ISCNT. ISLP 

coowj /app/ imT 

data JK/l/, JKX/1/, P1/3.1H1593/ 

COmjN /SOLVE/ FREQ. OaFTW. OELnX . EPSF, EPSX, EPSFS. EPSXS, 

FROrVU. CTEST. IPASS. KNTR. ISTRT. KSCVTV . 1W9CP, 
KNTTO. KHTSW. KNTWK 

COMO* /COMTAP/ IfREOT. FHEOTdOO), NOZATdOQ). GlKJOTllOO). 
CINJFTUOO), ITAPC. ITAPH 

COmM /W/ ICRT, IRFlAG.MEItl.lTYPE, IPPHtO. tX,AA{301 ,CM(301 .Xr.KO, 
UI30) .R1301 .»0.r30J.W(30J ,VF .VO.VOLF . vaO.ZF .ZO.OSAVEl 180) . 
FSAVCt|B8> 

FORNATl I8A>«) 

F(»UT(IH1.///.26X.‘FEED SYSTEM COUPLED STABILITY MODEL 
S0(.18AN./.9X. I8AHI 
FOWATt 1216) 

F0»1AT(/.5a(. 'IfTHYD •* , I2.5X, * lAPCOM - ' . I2.5X . ' IfPNOZ -'.12. 
5X.‘lTAf4l - .13.5X,*1TAPC . I3.5X. ‘ IT/PN -“.13./.5X, 

•IPfWYD -• .I2.5X.'IPIC0M ••.I2.5X. IPRNOZ .12.5X, 

MPBCHi --.la.SX.MPRSTE . 12.5X. NXP -'.IVI 
F0»*1AT(6E12.e) 

FORWTt/.SX.'XO .IPEll.H.TX. *XN02 .El 1 .H.5X. -RIKJ 

El l.H.5X.*CA«0 -•,CU.H./.27X,‘C0 •* .Ell.^^.TX, DELP •' . 

El 1 .>41 

format )/.5X,*0#»OT -'.IS.SX. • lt«T -■ . 12.7X . ' IMSKP - . 12. 6X. 

'ICKTra -•.IH.HX.'IWTTMC -■ .m.SX. 'IWTSMX -’.l>4) 

FCmAT(/.5X.'OCGA<R>-*.iPEIt.>*.2X.-ahEGAn)-‘ .El I -H.2X, 

‘FRQMAX •' .Ell,H.3X,‘0CLF1W .Ell .H./.27X. DCLriX 
EH .'♦.'4X. 'CTEST ••.Clt.H.//.5X.EPSF - .Elt-V.SX. 

‘EPSX •• .EH.'4.5*,‘ePSrS .El l.V.VX. 'EPSXS -'.Ell. *4) 



R-980e/a-20 



mmm 


Mjrvxxu cmrr wn - rscw omj.L, "•oown 


OMIT TITU - MM-PnOCOMUL STATOCNTS 


loao 

1081 

2000 


2010 


2020 


F<mat/.5X. PC - MPcn.^i.Tx.'MBOKi -•.cn.H,tit.'r«Fui 
EM.Ht 

ronrUT(M2.9CI2.8l 

ro»uTi/,5x.-i#TcoT -• .m.iox.'nean 

’FTCCrOC .CM .•*) 

r0»WlTt///,20X. XOeu5TI0N DtNAMIC CtWFICieMTS’ ./) 
FORmTISC.'COKf ■.12.') , IPCIl ‘ ‘ .EM .•» .5X . CFU( .12. 

*> •‘.EIl H.-.'.EM.m 

F0R1AT1///.10X. TBEOUENCY -'.P8.2.' KZ.'./.IOX. 

■DECJOCN7 -’.F8.9.//. 

lOX. NOZZLE AOniTTANCE -* .FB.S. ' . ' .2X.F9.5,// . lOX, 

'fEED system RESPONSE' ./ ,20X. '0X101 ZEB ,F9.S, ' . “ ,2X,F9-S, 

/.2HX, TUEU .F9.S. ■ . ' ,2X.rB.3.///) 

FORMAT nOX.^Cl IX. DSCIU-ATORT' I ./,SX, 'OISTAACE' ,6X. 

PRESSlflE RATIO ,7 X.' velocity RATIO* .6X. 

■temperature RATIO' .7X. ' hi X n«E RATIO' ,/.5X. • I INOCSr . 
m5X.'A#PLITUOE PMASE'I./I 
rOWIAT CSX .FB . >♦ . V 15X .FB .5.2X ,n .21 ) 


ORIGINAL 

OF POOR 


PAGE 'iS 
quality 









M/25/'n AUTOFLOM CHMT SET - FSCSM COmjTER PROOPW1 PACE 11 

CM«T TITUE - «N-PWCaXIUL SUTOCWTS 

CCrtfM /corwc/ MO». XdOOt. XMdOO), AdOO). OA(tOO), DCLX. 

xo, moz, AiNj 


ORIGINAL PAGE IS 
OF POOR QUALITY 



D6/25/7S 


<UTOn.W CMART SCT - F9CSM COtfVTCT PfWWH 


PACK I? 


CHART TITLE - SUBROUTItC OWCT 


/ / 



R>SM8/>-24 




























os/arn 


AUTOTLOW chart set - FSCSH comr/ter proorah 


PAGE t>» 


CHAm TITLE - NOM-PROCEDURAL STATETCMTS 


O0PPL£X*i6 AHA(H,5). CHAIHI 

COPFLCX CrdA, P, V, m, T. CNOZA. VXO. 

CCKl. Ctns, COX3, COKS, COX9, COX6. CO>a. coxe. COK9, 
COX]0. CCSUl. COXI2, COX13> coxm, COXIS, CQX]6, 

CfVJ, CFue. CTU3. CflH. CFL®. CFLC. CFl/T. CRJ0. Cni9, 
CFUIO, CFLIIl, CFLM2. CFUI3. CFUIH. CTU15. CFU16 
Cttf^LEX OHJ, 0«3. Dm. OHS. I>6. OtfTOX . OK7TO. DM80X. 

DMBnj. OlOx. oinj, caox, czru, G30 x. gstu. g^ox. gmhj, 
RBSOX. RBSFU. PINTOX. PJHTFU, RmTOX, RPNTFU. WNTOX, 
tWMTFU, VIHTOX. VINTTU, 1] 


REAL reOKI, rCTDI. mB. rnci, HGI. NUBOX, NJBnj. FMC 

comoN /ccrcsh/ xkox. xkfu, reoxi, hbfui. talbox, TAueru, vsox. 

VBTU. OAHO, RCO, DCLHOX. OCLPFU. PC. CO. 

COXl. C0K2. COX3. COXH . C(W5. COX6, C0X7. C0X8. C0X9. COXIO, 
COXIl, COXl?. C0XI3, COXl**, COXtS. C0X16. CTUl . CFU?. 

CFU3. CTlH. crus. CTU6. cruj. CFU0. CFU9, cruio, CFUIJ. 
CFUl?. CTU13. crum. cruis. cruis, pmg. xi*f>ru. xiw»ox. 

CS, OCSDPR. D»cm. DRQOMR, AOVOX. A£X»X . TQRAOO. DELVOX. 
►A«0X. OTOXOH. AOVru. AOOru. TORAGT. OELVTU, NUBFU. OTTlJOH 
CWON /C0H5TS/ r«8(lOO). TBHOO). RMOB<iOO). VBllOO>. 

OMBIIOO), DRHOB(IOO). OVBtiQO). VAPBOXnOOl. VAPBFUtlOO). 
S5V1IIOO). SSV?IIOO). 5SV3CIOO). SSV**UOOI. 
ssvstloo). SSV6HO0J. ssv7tioo>. ssvenooi. ssvooxuoo). 
ssvsFutiooi. ssvionoo). ssviuioo). ssvi?nooi. ssviKiooi. 
SSVI‘*M00J. SSViSdOOI. SSV16(tOO>. 
fMOGl . VGI , rfIGl . HGI 

com« /COMOM/ Pdooi, »«noo>. vd 00 >. miioo), tiiooi. 

VXO. OrCGA. CNOZA. OELP 

COHMON /COHAfC/ NXP, XdOO». XHdOO). AdOO). OAdOOl. DELX. 

XO, XNOZ, AINJ 

CO* 10 N /AOARND/ CJOX.G?OX.C30X.(?*OX.Gtru.C?FU.G3FU.GMFL.PINTOX, 
f**rrOX . mNTOX . V I NTOX .P I NTFX) . »»^rFU . mNTFU . V t NTFU . RSSOX . 
WSFU.DMl .DM3.Dm.OH5.OM6.OK7OX.OH0OX.Ori9OX.Dfr7ru. 

OHarU , DHSTU , I [ . OH? . OH?? . RHO I NJ 
comm /oit«>/ i*«T 
EQUIVALENCE (AHAU ,S) .CHAI 1 )) 

9000 FORMATl/// OTCGA - ' , 1PBE13.5//3X. • I ‘ . I IX . P , l9x . RMO , 
aCX. - W .aox. ■ V .?IX, 'TV/IX. 13. IPIOEll .*»J 
9100 FomATdx. i3.iPioen .**1 


R>980ai/i>2» 



AUTOTLOW CH«T STT - FSCSM COrW/TCR PROGRAM 


PACC 15 


CHART TinX - SUVCUT tic CMCmill 


/ ■= 7 °! / 








o'? ^ 



i 

C CXIT J 


R-9a0e/B-27 




















0K/S5/1S 


MiraruM OHWT 9CT - rsc» ctmiTCR pnxRm 


CMAirr Tmx - WOM wwcccuwc STATDCNTS 


C9f>lXX OrCOA. f>, FM), V. m, T. CNOZA, VXO. 

COX], COKE, C0X3, ocm. cots, COX6. C0K7, COXS, C0K9. 

COX10, COKM. C0KI2. COKIS, COXtH, COXI5. COXI6. 

crui . crue. cnc. crw». cns. cne. cfu 7. crue, cru9. 

CFUIO, crui 1 . CFU12. CFUI3. CFVIH. CTU15, CRJ16 
COt«>LCX oni, Dn3, om. ns, os. OrPOX. DKTai, otsox. 

OMSTU. GIOX. Giru. GEOX. GETU, G30X, G7U. G*«OX. OtFU. 

nasox, R8SFV, PiMiOtt. PiMiru, »p<tox. rwrox, 

rSNTTU. VIMTOX, VIHTFO. 11 
.OMBODC.OnSFU 

RtAL reoxi, rsFui. me, moi. mi, nubox. nubtu, r«c 

COTON /COMC91/ XKOX, XWV, tSOXI . fSTUl . TAU80X. TAUEFU. VBOX. 
VBFU, GAnO. ROO, OEliOX, OEl>Ct>. PC. CO. 

ccat. C0X2. C0X3. COKH. CODS. CODS, 00X7. COX0, C0X9. COXIO, 
COXU. C0X12, C0X13, COKIH. COXIS, C0X16. CFUl . CFUE. 

CFU3. CF\h. CflJ5. CTUB. CFl/7. CTUO, CTU9, CFUlO, CTUI 1 , 
CTUJ2. CFU13, CFUIH. CRIIS, CFUIG. F**. XlfPTU. Xirf»OX. 

CS. DCSCWt, DtCm. OROOH^, AOVOX, AOOOX. TORACO, CELVOX, 
tA«CW. OTODCM. AO^. AOOFU. TORAOF . OELVFU. NUBFTJ. DTFUOM 

Mf»04 /CONSTS/ meiioo). TBdOO). Fwoenoo). venoo), 

OHSitOO). ORHOBdOO), DVBUOO). VAPeoxitOO]. VAPeFUnOO^ 

5svt»ioo>, ssvEnooi, ssvaaoo). ssy^Kiooi, 

SSV5U00>. SSV6M00], SSV7II00>. SSveilOOJ, SSV90XtlOO). 
SSV9FU(I00>, SSVIOIIOOt. SSVIIUOO), SSV121JOO). SSVl3riOO>. 
SSVlHdOO). SSV15M00>, SSVlGnOO). 
fMx: . vGi . ^ncl . mi 

COPMON /coron/ pciooj. FMXioot, vnoo>. mriioo). tiioo). 

VX3. OrtCA. CNOZA, OELP 

comcN /COMAFE/ wcp. xuoo). xniiooi. a(ioo>. da<ioo), oclx, 

XO. XNOZ. AINJ 

COrPON /A0ARM3/ GlOX .CeOK.C30K.n*OX .CIFU.G2FU .G3FU .GM^U.PINTOX . 
ffr«TOX . FWTOX . V IHTOK .PI NTTU ,f»*(TFu . mNTFU . V J NTFU . Rasox , 
RBSTU.Om .0M3.0W».t*S.0rS.DK7OX. omox .OmOX.OMTFlJ. 
OrfiFU.DNgru, l ] .Om.DteE.RHOlNJ 



R-9806/S-28 




CALL EXIT 


CALL EXIT 














08/S5/75 


AUTonjM CMAirr set - rscsM cottuter ptooram 


PACE 19 


owr TtTlX - NOA-RROCEOPAL stateicmts 


OirCNSlOt A<B2.IS6) 

ja fWIAT tSIMlERflOR IN COOAO. 9U8WTI«. J AhO wTlAX EOUAL. RESPECT I 

«XY. . 2112 ) 

16 rORMAT reiHiriATRJX IS SIM9ULAR. EXIT FROM CXJGAEL. T>C PIVOT ELEICN 

TS ARE.2EI4.6I 


OMGINAI/ 

OF POOR euAure 


t- 





0B/25/'W 


MirorUM CHAPT SET - rSC91 Om/rCR PPIOORA^ 


PACE » 


CMAflT TIILE “ SUBROUTIhC OCr«Dr( IP.FRCO.GINJOK.OIHJFU, IPfCOM. IM>COn> 


/"ofcor / 


















oft/as/Ts 


<UTOn.OH CM*«T SET - rSCSM OOWTCR raOGfUn 


PACE at 


CHART TITLE - SUBROUTINE O»«0nr( IR.nCO.OlHjOX.OINJFV. IPRCOM, l»«OMI 



ON 

DUM ■ \ 

OLIi»Xlff>FU/XKru I 

RO « iaOvTU - 
0 . 5 "AOOFU»OUril 
•& 1 NJFU 

RTP - 0.5 j 

[ ot - : 

I a. «pi«rRE«TDRAor , 
05 

[ RFJ - 

[ 0.5/OELVFU*IDt» 
f 01 - 

(O.O.I.OI’On 
/IJ. + Ol'Dl) 

cnji • - Ru - 

a.«RD • tNUBTU - 
I a. )»C. 5 -« 0 /NUBFU 

—j. 














autotlom chart set - rscsM ctaMPurcR prograh 


OWRT TtTU - MM-RROCECUML STATETCNTS 


COrfUX GINJOK, OINJFU. DUM. NOZAT. OtNJOT. OINJTT, 

COXJ. CCXZ. C0X3. COXH. CO«. C0X6. COXT. COx9, CCW9. 

CQXtO, COMll, C0X12. COX13, COXm. COX15, coxie: 

cmi. cnje. crus. crvi. cno. crus. crvn. crue. crus, 
cruio. crun, cruiz. crui3. ctuih. crois, cruis, 

C09CT(I61. CFUTJ16J. RM. RU. RO. RFU 
REJV. tooxi. rflFul , ^wc. hubox. Mueru 

COmON /COHTAP/ FREQTClOO). NOZAT(IOO). OIHJOTUOO). 

GiMjrmoo). rTAPM. itapc, itaph 

CtttCN /C0HC8H/ XXOX, XXTU. reoxi. f«FUI, TAUeOX. TAUeru. VBOX. 
^FU. GMO. RGO. DELHOX, DEL»CU. PC. CO, 
coxi. coxa. COX3, coxh. co». C0X6, cor?, coxa. coxa, coxio, 
COXll, C0K12, COXI3. COXJH. C0X15. C0X16. CFUl . CTUe. 

CFU3. CRh. CTU5. Cru6. CTI/7. CnjB. CTU9. cruiO. CTUll, 
cnjia, CTU13, cfum. cruis. cruie. •«;. xihpfu, xitpox. 

CS. ocsorfi, DKW. onoow. aovox. aooox. tdraoo. oelvox. 

MU80X, DTOXDn. AOVFU . ADOTU, TORAOT. OQ-VTU, MJ8FU, OTfUDH 
COUtVALENCE tCOXTt n ,COX1 ) . {cruTui.cTun 
DATA PI/3. IHI99/ 

1 roniAT(flEi2.e> 

a FORMATt/SX. JOCOK , IPtn .H.5X. TAUBOX - . 

Ell. H.3X. ’VBOX -• .Elt.M.SX, 'DELHOX .Ell ,H .27X , 

•TORACW -• .Ell .H.3X. AOVOX - ‘ ,El I .V . // ,5X . AOOOX .Ell -H. 
MX. 'OELVOX O’ .Ell.H.lX.’NUBOX -' .Ell -V ,‘iX , OTOXDM 
Ell^./.aTX.'XlWW -'.Ell.V) 

1 roi*iATi/.5x.'XKru •■.iPEti.s.w.'TAuaru 

EJ1.H.3X.'VBFU -‘.EU.H.SX.’OELrt'U •' .Ell .V./ .27X. 

■TORAOT - .Ell.^.3X. AOVFU -' .Ell .V.//.5X. AOOTU <•' ,EI:.4, 
HX, ‘DCLVFU •' .EU .V.3X. 'MJBFU .Ell .‘♦.‘♦X. DTFUCM . 

Ell .H,/,27X, 'XlffTU -'.Ell. 41 
4 FOWWTl/.SX. 'WAJ -■ . IPtl I .4.6X. 'CS .Ell .H, 

TX.'ORGCPR .EM.H.3X. 'OCSCTR .Ell .4 ./ ,27x. OMDWl 
Ell. 4) 


R>960e/B-34 



Oe/25/75 


MfrOTLOW CHART SET - FSCSM COHRJTCR RKXWAM 


PACE 23 


C»MtT TITLE - SUBPWLTMhC COf1ATlA,t#TA,H» 


ROUTIIC SOLVES. AX - 
e. kiCNE A.X. AK) B 
AfC CCrfLCX. 

Ad.J) IS N BY N 

COrwxx matrix 
8d) IS N BY 1 
C9fL£X VECTOR. IT IS 
STORED IN A( I .N*!) 

FOR 

I - l.N 

AI«kCR IS RETUneO IN 
AII.N«-t> 1-1. N 

M%A IS T>c NireER or 

ROUS IN TIC OirCNSION 






R*98Q(/»-3S 














DB/asm 


MirarUM CHMIT XT - rSC9H CCrfVTERmR^ 


P»OL 


i 


owirr TiTte - tcN Fwo cr r u m. statdcnts 


CCmXX*IB AfMU.D.H 

9000 OCCK IM COmAT ••••••/ 

IR • ’.I10/‘ mTRIXAU.Jt • V 
(5X,ai6.IPCI9.6, ' .‘.lfCI3.61l 


R-9806/ft-3C 


I 



xm/n/’n 


JMJTOn.04 CHMT SET • rSC91 OmnCR mOORAH 


PMC » 


CHART TITIX - SUOKH/TItC nCSPnCRT.lMRITC, IX) 


nCQUCNCV RESPONSE 
nOUTirC HITH TI»C 
OCLAY 

PROOawtCD BY J. K. 


hunting, rockctottc, 


MAY 1975 





R>»806/»-37 





























































































AUTOFUW CM*HT SET - rSCSH CCr^UTUR moOMM 


PACE 36 


CMMT TITIX - «UBOfrilC nC9>MCRT. ItMIlt , IX) 










06/a/T5 


M/TOFLOH CHART SCT - FSCSR COrPUTER P»»GRAM 


PACE 17 


CHART TITLE - SUBROOTI^ FHESPf tCRT. IVRtTC, tX» 


ZZjZZ 



CeiiD 


CALL EXIT 


ZZZ7 

oe 

RETURN TO SYSTEM 


original 
OF POOR QUALITY 


R-9eoe/fi-49 



AUTQR.OM CHART SCT - FSCSH COffVTER PROORAM 


PACE 38 


CHART Tinx ’ tOt>PNOCOUML STATEICHTS 


DircNBiON QH<ioi>.a«CP5tioi)'.uiiou ,or6« lon.lOf laei .cRiea. ;26> 

OlrOeiGN LXl C3l . I M3) l3l .C7(3> ,C<62. lS6i .TD(K,83) .0( 106) 

0!rCN51QN )MN3(601.HOI60>.)aCAt<S0),XinAGI60).OAIN(lOII .rAZCUOl) 

oihi^GiaN VMR<G0i.vMm62).PHANi9.iou .oecscs.ion 

OPC^GltM ATDU26) . IOF4 10) 

CCmON/P/tM.NCTD.tA.MCQ.TRIO. IRPS.OH) .OMFL ,U . ID.C ,TD. ATO. 

FREQ. CIHXH .GtNJFU.PC .MO.tr 

CO^LXX NOZATUOO) .CINJOTUOO) .GIMJTT ( 1 001 .CIHJOK .CIKTU 
CtmM/COHTAP/tnCOT ,nCOT C 1 00 1 .NOZAT .0 ! MJOT . 0 1 MJFT . 1 TARN , 1 T ARC . 

ITAI^ 

STATEfCKT FUNCTION 0CF1 NIT ION, OELTA(OOOirL) -1.3 -ABSIQOOIFV I / AMAXl ( AaSIQOOlFL) . . lse93679>«C'38l 
STATEICNT FUNCTION OCTINITION. tNlT«OOOaFU -0.5 ♦ SIGNt O.S.OO-TAtQOOSFL) ♦ 

02Fl.)n 

STATDCNT FUCTION OCFINITICM. T>€TA(0003F1. .OOO+PL) 

CM)- SIONM .O.OOOWHBO^)* ATANI U ASSIOOO^TL I • <1 -O-OCLTA (0003 
FLI)> l/(AB5fOOO3FT.)«OE2.TA{Q0O3FL) n-nSO.0/3. tm5931I • (KO-OCLTA 
IQ003F1.))* 180.0*tXLTAIQ003FUl*U .OHJeLTAIQOtPtFLIMC .5 HJNITI-QOOS 
rui 

lao F0nf1AT(M2.13.19. l3.19.F0.O.nO.2rt0.O> 

100 FOMTmai-FUflER OF MJFFUI) IN FREOUCNCIES EXCEEDS 100) 

ISl FO»IATt6n0.O) 

103 FOFt1ATI30H-NUrfiCR OF OOLUPWS EXCEEDS 136) 

51 F0mAT(6lie> 

12H F0f*tAT(3lH-Nlf«R OF EOUAT10I6 EXCEEDS 6?) 

MS F0nnATI36M INITIAL FREOCCY IS ZERO) 

in FORMAT t iwntc NU«R OF FREOUENCIES EXCEEDS 1001 

SO FO»*WTt3U0,IS.I«.ri2.On 

M3 FOnuTlTeHl 

> 

« rOF*1AT(t6H 1.0. VECTOR /UHO.aHlHl) 

112 rORMAT(|>« EOUATION .13 /IHHO .1P1E1H.6)) 

310 rOR1AT(03HlHATmx OF TITC DELAY TERMS) 

116 FOmAT(?tH LS HAS A iCCATI^ VALUE) 

52 rOfVIATieriH) variable magnitudes UGINO A FRE 

QUENCr OF .IPEIV.6.94 RPS./ IHO .9«H VARIA8 

LE NUMBER MAGNITUDE PHASE /HHO . 

I27,IPeE30.6)l 

1001 FORMATITO Ml TVCSE ARE INTERTCDIATE RESULTS PRODUCED AFTER EACH MA 
TRtX INVERSION. /SMIirur FREQUENCY IN RADI ANS^SECQFO AlC CYCLES/S 
ECOO./SM) VARlAaX.OAIN(Oe),PHASE<DEOREES) VARIABLES 

TO. ) 

323 FOmATlTlMI INITIAL VALUES OF OOCFFICIENT MATRIX) 

321 F0f0lAT(//13.5(.lPeC20.6/ISM> .IPTE1H.6D 

1000 FOP01AT<//F20.*l.rW).5//llH.F8.2.F7.l ,m.r0.2.F7.l ,!V.FB.2,F7.l , m.F 

0.2.F7. I . m .FB.2.F7. n ) 

53 F0RMATt21HI FREOUEMCY .5d3H VARIABLE ,12,3H 1/2IH0 

RPS CPS .5(l»l DECIBELS PHASE )/> 

5N FOR>«TMPeE10.2.5taM /. OPaFS.2)) 


R-seoe/B-u 





m/totkm cmaht set - rscsri co^wter program 


PAGE 39 


own Tinr - suBwn^c «ao 


/ *« / 


BOUTirC TO PRINT 
KCADINO 

PROGRArtCD 6Y n. D. 
SOCMAN. Roa<£T0^r^C. 
MAY IST75 


/MJTE TO OCV ^ 
/ 6 / 
VIA FORMAT / 


C ) 






R-9eoe/6>5i 



06/SSm 


AUTOn.<M CM««T SCT - FSC» OOm/TW ntOORAH 


WlOE HO 


o«RT UTu - Nffi r»mrco.mL siatdcnts 


10 FQPMATnMI./////.i|»,*ilNW.VTICAL OCSOUPTION' .// . 

*mx,*rCEO STSTOi OatFLED‘»//.>*7X. ‘STieiLm rVOBL' ,////// ^ 

*»7x.*cc»»vrEn noccL‘.///.2Sx.'PnoaRMi nmc. fscw. rivvw. 

•SION. rwY 1975'.///.2St.'0GVEL(lFC0 Bt. n. D. SDtfVkN. . 

' J. K. HUNTING. MOK. H. FCRTIG' . / .Hex , ‘ AOVNCO) PROGRAHS. * . 

ROCKCTtWC - ./ .sax. -DIVISION or ROCXICLL INTE»**AT lONAL- . / . 

<*ax. -OWOCA PMC. CALir 9I3DH-,///. 

aSX.-SPWEERED 0Y. NASA/unOON B. JCMCON SPACE CENTER 
Hex.-mwsTON. texas ttow* ./.nax . -under contract nasb-ihjis- » 




QOau^ 


R-96(W»-52 



AUTCrUM 0*nt *T - COfWtCT ^WOOIUM 


CNMIT linx • MjOCK data 






R-9808/ft-U 



M/T0n.0M CMwn SET - rscsM ccwter pnooRAn 


PAGE «*2 


0¥m TITLE - MM-moCEZUUL STATOCNTS 


OOmW/KT/ ICRT , IBFLAO . I TE»»1. 1 TYPE , Ifc« I TE . IX . A .CW.KT . KO ,L .R .RMOL . V . 
vr . VO , VOLT , VQLO . zr , 70 .OSA\C f t as > . rSAVC ( 1 8B ) 
seal AOOI .CMC30).Xr.KO.LI30l .Rl30) .»CLt30l ,V(30I 
data lCRT/0/ . tSFLAG/0/ , ITERH/0/ . ITYRE/ 1/ . IJ^ITE/O/ 


R-9808/B-54 



06/25/75 


MTOTLOW CHART SET - FSCSH COP*^(TER PROGRAM 


PACE tS 


CHART TITLE - SAJ0«a/TII>C KVORDYC IR. lATHYC.nc .GINO.GIAF .PCIN.UOtN.triN) 


/ KTORDT 


7 


9UBR0UT1AC 

MVDMOVnR.IA^mi. 

nVQ.OINXX.OlHFU, 

PC.UO.AT) 


SUBPROGRAM TO 
C^OJLATE FEED SYSTEM 
RESPONSE 

PROGRATtCO BY J. K. 
hunting. ROCKETimC. 
MAY 1975 
AROUrCNT list 
VMIABLCS. 

IR - DATA FLAG, 
-O-NO DATA READ, 

• r-REAO NEki DATA 
CASE. 

•2-MOOtrr EXISTING 
DATA CASE 

lAPHYO - MOOe FLAG. 
•1-CALCULATE INJ, 

FLOW GAIN * PHASE AT 
I FREO. 

•2-CALC. ♦ SAVE INJ. 
FLOU GAIN « FHASE FOR 
FREQ. RANGE 
FREO - !M*LT 
FREQUENCY FOR CASE 
WITH II^HYO-1 
CINXK - OUTPUT 
cor^x NurcER with 
A rtVlTUDC AND PHASE 
angle of 

OXID. INXCTOR FLOH 
OSCIUATIONS AT 
FREQUENCY. FREO 
CINwFU - OUTPUT 
COrPLEX NUrOER WITH 

afplituoe and WASE 

ANGLE OF 

FUEL INJECTOR FLOW 
OSCILLATIONS AT 
FREQUENCY. FREQ 
PC - STEADY state 

0PEHAT!^C chamber 
PRESSURE <LB/IN«*2I 
UO - STEADY STATE 

OXIDIZER INJECTOR 
Flow (lb/scci 
if - STEADY STATE 

FUEL INJECTOR FLOW 
tLe/SCO 


LABELCO COmON 6L0CX 
/COHTAP/ VARIAa.CS, 
FREOT - array of 
FRCOUENCICS i IDO 
MAX. » FOR l»#>Krt>-2 
CASE 

rTREOT - NLMER OF 
FREQUENCIES IN ARRAY. 
FREOT 

GINJOT - TABLE OF 
COTPLEX NUrCERS WITH 
GAINS Al« PHASE 
ANGLES OF 

OXID. INJECTOR flow 
oscillations AT FREOT 
FREQUENCIES 
GINjrt - TAB.E OF 
COTPLCX NLMBERS WITH 
CAINS Arc PHASE 
ANGLES or 

FUEL INJECTOR FLOW 
OSCILLATIONS AT FRCQT 
FREQUENCIES 


NMCLIST /HYO/ 

variables, 

ICRT - »^0T FLAG. 
-O-NO R-OTS, -I -PLOT 
GAIN«fHASE VS 
FREQUENCY 

IMTITE - PRINT FLAG. 

2-ND PRINT. 

— I -PRINT OUTPUT 
Oi-Y. 

-0-PRirvr ir#»UT OM.Y, 
-I-PRINT INPUT ♦ 
OUTPUT. 

•2-PRIMT INPUT. 

OUTPUT ♦ INTERTCOIATE 
CALCULATIONS 
IRFLAG - READ FLAG. 
-O-CARO irPUT, 

• l-TCRMINAL VfMJ 
mon ITERM 
ITERM - TEWINAL 
ir«n;T FLAG, -O-NO 
TERMINAL, >0-TEfPllNAL 
IHIT NO. 

NOTE, ITEmiD CHATOES 
IRFLAG TO 1 
1TYPE - system TYPE 
FLAG. -1-BOTH QKIO. y 
FUEL SYSTEMS 
□ESCRIBED 
WITH I data read, 
-e-OXID. SYSTEM DATA 
READ AND calc- 
ulated FIRST. nCN 
FUEL SYSTEM DATA READ 

• CALCULATED 

ID - OUmr NATC 
TO ALL&4 USE OF 
SEQUENCE NUTBERS ON 
NATCLIST 




! 


I 

j 

I 

1 


I 

I 


i 

I 




I BEGIN DO LOOP 
I HO I - I. 30 


^.01 — »| ' 


I R<|i - 1 000000 

I V I n - Hoooo . 

AM 1 - I . 
CWM » • 0. 
RhOL< I J - -OH 


H O ; 

I L < U ■ - 00 1 



26 




R<9aOB/B-S5 












06/2V75 


AUTOFLOU CHART SCT - fSCSH COWTCR FROORAM 


PAOC 


title - SUBROUTItC WCROTUR . IK^^TD.FRE .CtNO.OlNT.PCIN.HOIN.vriNI 





R-980e/B>57 















06 / 25^79 


4yrOFLOW CHART SET - FSCSM COfVTER WCORAH 


PAGE HG 


€»«Rr TITLE > SUenOUTMC KVOVT^R,l^f>KTD.mE.OI^O,01^F.PClN.UOtN.vriN) 



^»*7 01 


ft>900S/»-98 






















os/»m 


mitofum 04 AHT scT - rsc» ctmrrtR piosram 


CWirr TITLE - SUWOWTirC KrtWOr<W.l»'PWO,fWE.OlNO.OI>r.PCIN.MOIN.HriH> 



l-SBOB/B-SS 
































06/»r» 


autotlou cnwt XT - rscsM cor^tnxR muorah 


0«WT TITLE - NON-»>ROCCXM»*L STATEJCMTS 


COftW/HT/ ICRT JXLAO. ITOW. ITWE. IW1TE . IX.A.CH.KF .ICO.L .R.WO. ,V. 
vr . VO .VOLT ,vaO ,ZT . 20. OSAVE n8B> .rSAVC* 18BI 

REAL x(iee) 

EQUIVALENCE {XU> ,AI I > » 

CO*«N/r/h«.NCTO.^«.KEO.TRIC.IRPS.OHI IM.C.TO.ATO, 

FXQ.G INJOX ,e IHJTU.PC .HO .»»■ 

COff>LEX NOZAT(IOO) . 01NJ0T( I00» .OINJTT* UK)> .OiNwOX .GIN^FU. 

CINO.CIfr 

COrtCN/Cmr AP/I^TCOT . rXOT ( 1 0O I .mozat .q i nxt .oinxt.itapn.itarc. 
itafm 

XAL L130I ,R«30) .A<50) .V«30» .CI62. 126) ,TD<62.63) .ATOt 126) .H( 101 > 

.CUl30I.MLt30>.KO,Kr 

INTEXR*'4 IMI126I 

NATCLIST/HVO/L .R.A.V.ONI ,0W\.,FREQ,^F1C0T , imiTE . icrt.zf .rf.volt, 

VF ,20 .flO . VOLO . VO.FREOT . IM. 1 TCTM. IHTLAC.XO .XT .CH ,RMQL . 1 0 . 1 type . IX 

90 FORMAK- IKPUT NAfCLIST •+ftOOATAV VARIABLES ARE . L ,R.A . V .CH.»C 

L.t#'«0T,rKQT.1H»lTE.lTCRM'/- ICRT .2F .RF , VO.XOLF .VF.KF .ZO.RO.VOLO 
.VO.IO.ITVPC.IHFLAG) 

66 FOWlATr £«3 OF *NYt) IWVr*l 

WO FOWIATIlMl.^//.iex.-FEED SYSTEM RESPONSE PARAfCTEHS ’ . / / . 

SX. TREOUENCr .‘♦X.’OXIOIZER IHJCCTION RATE .'iX. 

■FUEL IHJCCTKX RATE" . / .20X. ' A^f>UTUOC* .5X. PWAX' . 

6X. Am-ITUCC’ ,3X. -^<AX- ./.i3X,0PF9.3.5X, IPEII.V, 

OPT 9 , 2 , SX , I PE 11 . . OFF 9 . 2 » ) 




R-98QB/»>il 







06/25/75 


AUTOTLOH CHART SET - rSCSM COWER PROORAM 


PAGE 52 


CHART TITLE - NQN-PR0CEI3LRAL STATEICNTS 


OirETGlON TK(M 

nWUTdHl mx 27« RROR IN TA8LE> 

FORHATUMO HIX 2TMCFER TO SteROUTPC LOCFAC // 

5X >« - IPEI5.R / **X SMTX - 6EI5.H / (8X 6E15.H) I 
FORIATIIMI 22X 8^ RROR - EXTRAPOLATION OF TABLE IS BEYOW R 

easonable limits > 


R-980e/»^ 



o^/zsm 


AUTOTLOW chart set - rSCSH computer PROC3RAH 


PAGE M 


CHART TITLE - SU8R0UTl»C NOeACTK IR.OAW.CO.fWQ.NOZA) 



R>9a06/i^ 


I 

























oB/esm 


m/totlcm o*rt sn • rscsn computer prooium 


R*OE » 


CMwrr Tinjc • fiUOTOuTirc moz«m( w.ow.co.nco.NOW* 



R-MOa/B-67 














06/2V7^ 


AUTOTLW CHART SET - FSCSM COUVTtR PROGRAM 


PAGE » 


CHART TITU - NOH-PROdDURM. STATOCMTS 


ir^lCIT fCM.*eiA-H.O-Zl 

COf*>l£X HOZA. NOZAT. OIMJOT. OIM^T 

REAL nCO. OttK, «3Qt. RCTX. ANOLEX. CHR, RINJ, CO. EREttW. 

nconi, FREQT, FlEQl. ^CQ^ 

OtrOCION OTCS.M). G15>. 0PT51. Yt5) 

CWON /Xt /CAM. SVN, ANGLE .RCT.RCC /X2/T .RT .Q.R I .R2,MC . IP 
/X3/ZIR.Z1I /XH/Oi 

COrttJN /COMTAP/ WTTEQT. FREQT(IQO). NOZATCIOOI. GlNJOTdOOI. 

GIMJFTCIOOI. ITAPN. ITAPC. STAPH 
CtmjN /COTfOZ/ RCCX, RCTX. ANGLEX. CRR, RINJ, IUSNQZ. FRE01X , 
EREOHl. IPRWZ 
9000 FonmTiGCia.S) 

9010 FORMAT I /.5X. -RCCX .ITCH .H.5X. RCTX • ■ ,E 1 1 .*♦ ,5 X . -ANClEX • , 

Ell.H.JX.CRR -‘.eil.Ht 

toco FC»1AT(1HI.^///.30X.SOKT>CORET1CAL WZZLE AOt1lTTANCES.//.Z3x. 

tMWWCM *Lt«ER • .F3.Z.-JM 5VN - .F6.H.»I GAWA - ,F5,3.// , 
7X.15W0ZZLE AN3LE « .FV. t .ex.SlKflAOI I OF CURVATLRE. 
.9KTMTOAT • .F6.4. 1 2H ENTRANCE “ .F6.N.// .9X.2^^C, 
7X.2HTR.e5t.2HY|.ex.l»r.ex.>4SYR.8X.3HSYI . 

6X . 9HALPMA . 5X . »«t«ETA . / I 
1005 F0RMAT16X.F6. t .2FI0.5.ri0.2.«*FI0.5) 


ft-9808/B-6e 











aS/2V75 


AUTOrtOM OiART SCT - FSC» COMPUTCR PROGRAM 


«*ACC 58 


CH«T T|ll£ - WCT) PWQCiaURAL STATDCKTS 


lt«.ICIT RCAL^BCA-H.O-Z) 

COrtOt /XI/OAM.SVN.ANQLC.nCT.RCC /JG/T.RT.O.Rl .R2,WC. IP 

oircreiON 0(5) . OPtsi 

9 raWAT(3X. PRINTING FROM CARD '♦570' ,/ ,3X, R- .E15.8. 

3X. “Rl-‘ .EJ5.0.3X.-RT--.E15.9) 




R-9808/B-70 











06/25/75 


AUTOrVOW OWRT SCT • TSCSM com^CT program 


PAOC 60 


0«RT TITLE - NCII PWOCCOJRM- STATTTCNTS 


IWLICIT fCM.«eiA-H,0~Zl 
COmON /ie/T.ir,Q.Rl,FC.HC,IP 
ortoeiCN u(5i. a(S). tizisi, rzin.S), oum(5> 
1000 rareiATm6x,F6.H,tx.rio.s.3x,FtQ.5) 


lt-980e/i.72 













oe/a/w 


AUTOTLOW CHART SET - TSCSH CCmHER PROGRAM 


PAGE 68 


CHART TITLE - MON-RROCECUUC STATETCNTS 


tT*>LIC[T REAL^SIA-H.O-Z) 

COrtCN /XUGAM.5 VN, angle. RCT.RCC /X2^T .RT ,Q.Rl . R8.HC . IP 
/X3/Z1R.2U 
OirCNSION 0(9) . GP<9) 

16 F0RHAT(3X, ’PRINT INC FROM CARO 5000 “./ .3X , R- E 19.8. 

3X. Rl»- .EIS.0.3X. RT-’ .C15.8> 



R-M08/»-74 



06/25/^ 


M/TOn.OM CM«T SET - FSCW COm/TBR PNOORAfl 


PAGE 63 


0«Mt TITLE > SUOna/n»C SOLlMflCMCJC > 




sl bf wo ota m solves for 

TVC CXPFIXX FICauCNCY 
BMED ON T>C 
UPSTNCAH MO 
OOM«TVCAM NOZZLE 
AONtTTMCE 
PROORMfCD BT K. H. 
FDtTIC. ROCttTDYIC. 
NLY jg75 


MN I9CNT > l.H. GR 
5 THIS PROQRAN 
COtfVTES CNOZA-fOZA 
FOR A 

CI\CN OrCOA. MN 
ISCNT • t OR I. THIS 
R0UTI»C ALSO Flics 
inAOIOFCOAl THAT 
mNiniZES CNOZA-NOZA 
ICEPING RCALIOrCOAJ 
COST ANT 


FUmCR. lACH ISCNT * 
I GR H, TESTING IS 

PErroRrco to 
DETEW1IIC 

IF A ROOT IS fCARBT. 


XI - AirWG(OKXM) 


XIHAOr com/TES 
F>#»»I0, 1) 
•FNI-OO2A-N02A, 
HV<AeS (CNOZA-NOZA) . 
htCN ISCNT * I OR *4 
AlO ISLP - I , 
XlrUOr«OtHV««2) 
/OUMAOlOFCGAl 1 




NAS 

JUST SEEN COrRUTED IN 
XIHAGF. GO TO t>46 TO 
FINISH CCFPUTATION 
or JACtBlAN AM) T>CN 
CtmiTE FTEST2. 








06 / 25/75 


MAVLOM 0«MT «T - 


OUrfVTUt mOQMH 


CHART TITUE 


SU8R0UTIIC SOLVHflMCRCI 



DOfC Kf> TO STATETCNT 
l**5. 



R-9606/l-7§ 










06/2Sm 


AWTOnjOM CM«T 8CT - FSC9I CW«ITO» MttORAn 


owr TiTi£ - a»wuTi« sxvwr»*ciC) 


CSTIIMTC »OT OUCSS 
USIfO rM.SC POSITION. 

6H.OI* — N 

t» i 0> 

X3 • IF2*X1 - 

ri*)»i/«r2 - rn 


IS 



•.^aiGINAL page is 

OF POOR QUALTTSJ 



WOT IS BRACKETS SUT 
CONItCRGDCC NOT YET 
RCAChCO. CALL ZERO 
TO nWJ ROOT WITHIN 



ctmffc ocuivATi^s. 

rrcsrecrcsT 

ruNcnoNJ. 

OETZIOETERNIMAMT 
OP .MCCBIAN). NO 
COCeiCOMMTUM 
Nlf«ER OF JACOBIAN) . 



I !? 

|AA1 - DFR0K»«2 ♦ I 
I DFRDY»*2 I 


AA3 • 

fOPKWoriox ■* 
orPOT*DrioYi«*z 

^ 


16 


coroa • <AAl ♦ 
AAZ • SORTUAAI - 
AA^>••^ ♦ 
H*AA3))/<AA1 ♦ 
AA2 - SORTMAA) - 
AA2)»*2 ♦ H*AA3>) 


SCRTlOOKei 


17 


KNTR - IMTR * I 




list • i9on. 

»JTR. «CNTH. 

, O^EGA. 0«ZA. 
HOZA. KN. FTSTe. 
DET2, COK» 




R-980a/B'77 














MJTCFLQN CHMtT SET - FSC91 


ammj* proormi 


Rage b£ 


OWRT Tinx - MROUriPC 90LW<MKIC) 





(SCNT-1 Al« 
ISTRT"l(NOT FIRST 
Tire nWOUOH OR FIRST 
Tire T>«OUCH 
AFTER A SIMSULARITY 
IN T»C JACOBIAN HAS 
BEEN CR05SEDI . TEST 
FOR 

EXtSTANCE OF ROOT IN 
CURRENT FREQUENCT 
RANGE 

<FRE0-0ELFRa,F«0> . 


FT5T2*FTSTl ,L£. 
0. .AND. 

GET? •OCT I .GT. 0. 
.AW. 

AHAXl tC(N>! . 
COMKI .LE. CTEST 



S.9i H 

l>*7 ' 


FIl - FI2 
Xll • xie 
FRi - ne 
XRl - XR8 


FRS “ FR3 
XR2 - XRS 
ri2 - FI3 
xie - XI3 
FTSTl - FTST2 


cowi - coce 

XTl - 0ET2 


IF ISCNT - 1. UPDATE 
FRED BY OEIFRQ. IF 
ISCNT - FREQ HAS 

already been changed 

BY rETHOO OF FALSE 
POSITION. 


A8S(FR3I 

/CABStNOZA) .LE. 

EPSFS .AfC. 
AHINI (A8SI <XR3 > 
XRl >/XR3l . 
ABSU)«3 - 
XR2}/KR3n LE. 
EPSXS 



FRI - FR3 
XRl - XR3 



nCO - FREO ♦ 
OELFRO 


LIST 

lER 

KNTS. 1 

KNTR. 

ISCNT. XRl . 

XII, 

FRI . 

FI 1 . 

XR2. 

XIE. 

FRE, 

FIE. 

XR3. 

X13. 

FR3. 

ri3. 

XfN, 

XIR 

r«R 

fir 


tentative root mas 

BEEN BRACKETED. 
continue with ISCNT - 
•♦n -E. 

FREO MILL BE OHANCCD 
By false POSITION) 

until either it is 
OCTEWTITCO That t^c 
JACO eiAN IS SINGULAR 
OR THAT A ROOT IS 
actually there . The 

LATTER IS ASSLPEO BY 

default ir 
derations have been 
PEWOHMEO HIJhcuT ’"HC 
JACOBIAN Chang INC 
SIGN 

AND ThC condition 
NU!«R REfTAlNS LESS 
Than ctest a.*® thc 
ERROR 

!CA8SICnoa2a-noZa) i 
IS less Than 2. 5* 

EPSF 


IP • 0 

OrCGSV • OTCGA 
DET25 - DET2 
FT5T2S - FTST2 


I cowes - cowe 

FR3S - FR3 
j XR3S - XR3 

j FI35 - FI3 

X13S - XI3 
. KSCNTH • 0 

67.17- - 

181 29 

I FTSTO • FTSTl 

' DCTO - OETI 

CONDO - COWl 

IP • MAXOI tP.Ol * 

I 


67.09 

!82 

*~itSCNTN - KSCNTN < 


tF7ST2«XR2 - 
; FTSTfXR3i 
I /(FTST2 - FTSTl) 



Of€CA * XfN ♦ 

to. .1.) 

•AlHACIOrCDA) 

FRED • 
XRH/(2.*PU 


0 ) 

^IN7 


R-9806/B-78 








OB/ 29/75 


AUTOTLOH chart set - rSCSM COHVTEJ* PfWOlUH 


PACe 67 


0 «irr TITLE - SUeKX/T|^C SOL>Min*CRC) 




S'JAH PAGE IS 
:-OE QUALITY 




























DB/29/79 


AUTOTLOM CMAflT SET - FSCW COWVTW WWJRAM 


PAGE 69 


CHAMT TITU < SUBROUTlhC SaVMOMdC) 


MMCfC CCTtRHIMIS 
TWCTEH LOACATION IN 
CALLING PROGRAM 
UPON RETURN. 



/ 5000 / 


».M* >i 


03 


KH-CRE - 5 


EXIT j 




06 / 25/79 


AUTOTLOH CHWT SET - rSC» COr**UTCT PAOGIUH 


PAGE 70 


CHART riTlE - N O N- PR OCE DURAL STATE^ITS 


CCr*»tEX OrCOA, P. »WD. V. W. T. CNOZA. VXD. 

NOZA. crcoo. 

oori, 21. 22. 23. FH, NOZATR. ON. 

CNOZO, NOZO. FHO. ONO. 0FCO5V 

ccmjN /coHOti/ piioo). HHOdooi. v(joo). wnoo). Tnoo>. 

VXD, orcoA. CNOZA. oeup 

COmON /rZERO/ H0ZA.N02>W.0N.FN.F»«.rN( .HN, JSCNT.ISLP 
EXT17MAL XIMAOF 
DATA Pl/3-l»<1593/ 

COMMON /SOLVE/ FREQ, OELFRQ. OELTK, EPSF. CPSX, EPSTS. EPSXS. 

mOMAX. CTEST. IPASS, KNTR, rSTRT. KSCNT4. lUSKP. 
KNTreC. W«TO«. KNTWtt 

V>t90 F0RmTl/21S.lPNC13,5/2X.3nPEl3.5. ■ . ■ . IPE12-5) . IPEIS.S* 

860) F0R1ATI//' UNABLE TO ri»0 ROOT TOR IMAO PART Of f •••• 

//‘ XI ,Fl ,X2.F2.X3.F3./»6.FANS.KNT.IER.0«:CA - V 
3X,1P8E13.5/3X.2I10/3JC,IPE)3.5. ‘ . ’.1PE53.5) 

8500 FORMAT t/JIS.SUPEm. 5.' . ‘ . )PE12.S)/2X. )P*fEl3.5» 

8S10 F0RMAT(2X. IPEm.5, ■ . * JPC)2.5.2X. tP^EI3.5) 

»*5C FOWtATl//- 

MARNINO. POSSIBLE ROOT IN FREQUENCY RANGE. . 

IP2E15.6//' 

8515 format (/ 2I5. 3 MPE 14.5. . ‘ . 1PE12.5) , JPE14 .5) 

8520 rORMAT(2X.4(lPEl4.5. ‘ . •.lPC12.5)/2X.miPE14.5,- 

IPE12.5)/2X.)P£IH.5. ■ .-.IPEtB.S) 

8802 FORMAT)//- •••• EWEEO CONVEROENCE LIMIT »•••■/ 

lER.KNTS.KNTR.ISCNT • '.HUO/- X.F FROM i -H - ■/ 
t3X.lPE15.5. ■ . -.IPE15.5.)PE18,5. . '.1PE15.9)I 


R-90Oe/B-<2 



PACE 71 


CMMIT TinX • 9MRaJ1\*€. STOOnrUPRSTC) 





OF POGB QUAKT? 


£-«lfla/B-83 






















06/«5/T5 


MiTQFLQM OWTT StT - fSCSM COm/TtR moORm 


PAGE IZ 


0*tRT TITLE - NOH-PROCEtXlUL STaTETCNTS 


Ctmxx COKl. CW2. COM3. COKH. C0K5. C0M8. COTT, COW. COX9. 
COKIO. COKIJ. C0K12. C0K13. COKIH. COXI5. C0K16. 

CTui, cne, cnn. cfv». cfvs. crus. cn/r. crxje. cfu 9. 
cruio. cruM, CFVI2. cpuis. cfuih. ctuis. ctliis 
peal rWBtNJ. mOK\, «R>| . mB. mcl. MOl. NUBOX. 
ctm» /cac&i/ nem. fsoxi. tmjbox. jmov. vbox. 

veru. OAMo. rco. delmok. OEL>rg. pc. co. 

COXl . COJC. cow. COM**. COX5. C0X6. C0X7. C0X8. COX9. COXIO, 
coxn. C0W12. C0X13. COXl>*. COX15. COXJ6. CTUI. CTue. 

CFU3. CELT*. CEU5. CFU8. CTU7. CTU8. CTU9. CFUIO. CFUl I , 

cruia. CFUJ3, ctuih. cruis. cruie. mho. ximpfu. xl^Pox, 

CS. ocsom. OKHT. ORGOMR. AOVQX. ADOOX. TORAOO, C3EIVOX, 
M^OX, OTOXCM, AOVfU. AOOFU. TORAGF. OCLVFu. NUBFU. OTFU>i 

com>j /^coNSTs/ rfeitoo), tbciooi. rhqbuooi. vefioo). 

DffOMOO). DRHOBUOO). OVSflOOl. VAPSOXHOO), VAPBFU<100>. 

ssvKiooi. ssvsmoo). ssvsnoo). sswdOO), 

SSV5(I001. SSVS«>00». SSV7nOOl. SSV9<IOO!. SSV90XUOD). 

ssvgFunoo). ssviouooi. ssvnuoo). ssvipdoo). ssvijtiooi. 
ssvimioo). SSV15HOO). ssvietioo). 

fMXji . vGj. rtwi. nci 

CtjmON /COMAHE/ NXP, X<t00>. XMUOO). AdOO). OAdOO). DCLX . 

XO. XNOZ, AIMJ 

F0RMAT( JHl ,///,3 ix.-STEAOV STATE SOLUTION'.//. 

7X. -distance- . 

‘♦X. • TE»f>CRATU«; ■ .a. - VELOCITY* ,3X. -MIXTURE’ . 

6X, - DENSITY - .7X, 'PERCENT VAPORIZED' ,/.7X. ’ ( INC>C5I ' ,5X, 
tRAMCthC) - .**X.- IFT/Sl • .SX. RATIO' ,5X, 

• (l9m/FT«* 3) • ,6X. 'FUEL - .a. ' OXIDIZER' . / ) 

FORMAT C6X,F9.*t. a. 2Fl I ,Z.Fd •■*. 1PE15.5.0PF1 1 .5.FS.2) 


R-9808/B-64 






















06/»/75 


autotlow chart set - rscsM ccmrrER program 


RAGE *79 


CHART title - SLBROUTjrC TAOA«<N.H.X ,Y .OV . t02. 105 



t-980B/B>87 


ORIGINAE PAGE IS 
OF POOR QUALIIY 










06 / 29/79 


AUTOFLOW CHAHT SET - rSCSM Ctm^TCH WWOlUM 


FACE 76 


CHART TITLE - MON-RROCEDURAL STATDCNTS 


im>ICIT ICAL«e(A-H,0-Z> 

COmON /Xl/0A/1.Sm.ANQL£.nCT.RCC /X2/T.RT,Q,RI .RS.HC.IP 

ctmx /XH/01 

0!^e>6ION C0R<91. DP<5», DY(5,m . PREDtSi . T(51 . 0(9), GP<5) 
16 FX3RmTI3X. ’PRINTIMD nW! CARD 2180 ‘ , / ,3X . '»N • .E19.0. 

3X, 'Rl- ’ ,E19.8.3X, 'RT*' .E19.8> 

61 rO»8UTl3X, ’PRINTINO FROM CARO 2970', /.3X, R«’ ,E19.0, 

»,'Rl-' .El9,B.»t.*RT--,El9.e) 
lOOO F0niAT4H€X,F6.HJX,P10.Q,3X,ri0.5> 


























M/tcyuM cMwrr set - rscw co»wict wooram 


PAGE T9 


04WT Ttns 


- NCH-mXtDURAL STATDtNTS 


OlfCNSiON MUOl I .VdOU 


ORIGINAL PAGE IS 
OF. POOR QUALITY 


R>9B08/B-91 


06/Z5/T9 


AOTOFUJH own 9CT - rSCW COfVTOI PflOOUM 


f*MC m 


own TiTi£ - ntcTicN xiMorixnof}) 


WJTlfC TO COWTE 
MiZZLC AOn. DIFT AS 
function of 
imciNARY PART or 
OrCOA 

procrattco by K. H. 

FCRTJG. ROCKETCnnC. 
HAY 1975 


fcamnc ocga using 

ICM IHACIWRY PART. 



CCmiTE 
DICN0ZA-IC2A) 
/0(|HAC(0rCGAI>, - 
CALLED GN. AhO 
DfCAaS(CNOZA-NOZA» 
••2l/0<IHAC(0rCCAn. 
called xinAcr. 



GN - <CN - FNI/DH [ 

xihagf - I 
- 


R-9808/B-92 










oe/2sm 


M/rOFLOH 0«MT KT - F9C91 UfW/TCT FHOOWT 


CHtftT TITLE - KN-WOCOMWL STATDCMTS 


omxx 0»COA.HOZA.<»C2A.P.»«.V.»«.T,VltO,rN.««r«.ON, 

COOCl. COG. CCOCl. COKH. COO. COO. COa, COOO, CCD«9. 
coatio. coni, cotia. coxi3. con**, cons, coxis. 

CfUI. OrLB, cru3. CRAt. CFVO. CfXJB. cf\n. crvo. CRS. 
cruio. cron. crui2. cruij. cfuih. cruis. cruie 
.orcoo 

SEAL fOOXI. TOrui. NUBOX. W^TU. f*0 

COmON /C0rc»1/ WO(. »CFU. »«0X1, mtoi. JMOK, TMOTU, V9«. 
VBTU. oa«, ago. OELHOX. OCL^FO. pc, CO, 

COXl. COG. cots. COP*, COO. coo. coo. COO. coo. cotio, 
awn. cotis. coci 3 , coxih. ookis. coti6. crvn, cn«. 
cr03. CFV>*. CFU5. CTU6. CF07. CTUB. CTU9. CTUIO. CFUll . 
CTU12, CTU13. crun, cruis. crois, mjg. xiwtu, xiwoot, 

CS. DCSCf«. OMOR. ORCOH?. AOVOX. AOOOK. TWAOO. OCLVOX. 
NLBOX, DT0W3M. ADVTU. XCOTO. TORACT . DELVTU. M^U, OTFUJM 

OmiN /FZ0»/ NOZA.NOZAm.ON.rN.F>«.rNl.»tJ.JSCNT.lSLP 
comoN /coMcm/ pciooj ,rho< lOO) .vnooi .»onooi .t< ioo» . 

VXO. OCCA .CNOZA .OCtP 

corrON /coM«c/ Nxp.xnooi .xnc ioo» . auoO) . o*noo> .ocut. 

xo.moz.AiNj 


Pact? to 


R'9eQB/I-93 







































oK/s^’n 


AUTOFUJM 0««T SCT - FSCSW Ctm^W FKXMM 


PAOC 05 


CHWT TITLE - MON-mXEOlVUL STATTJCMTS 


Jh«.ICJT REAL'S fA-H.O-Z I 

comw /X1/0AM.SVW.ANQLC.RCT,«CC /XS/T.RT.O.RI .R2.HC.IP 
/XH/CH 

OltCNSION CORtS). 0Pf5). OVfS.HI. PHEO«51. Y(5> , 0(5). 0P<5) 
16 rOWATlSX. PPlNTlNO FROM CARO >30' .E15.0. 

3X.<RI-’ .Ei5.0.3X.'RT-' .E15.0) 

61 FORMATO, ‘PRINTING FROM CARD 3790 ‘ 3X . *R»* ,E 15.8. 

3X.'Rl-‘ .E15.8.3X. ‘RT*- .E15.0) 

1000 FORMATIHex.Fe.N.lX.FlO.S.SX.FIO.S) 



R-9a08/B-97 





























0S/a5/79 


CHARt TITLE - m i W M C Ca WW. stateicmts 


STATOCKT RfCTION OCTIMITJCW. TCK) • SM^UXI 


1 ' 

_<gsS'" 


R.90oe/i-ioi 


DO or MITOrUM CHMT KT J.JTJi IPWT STaTDCMTS FfBCXSStD 


E5CCUT10N TlfC - 


t mw at SEC 


OFPM^ B 
^ quality 


R-9806/8>102 



APPENDIX C 


COMPUTER CODE LISTINGS 


INDEX 


PAGE 


MAIN 

AREA 

CHAMDY 

CHMC0N 

CO0AEL 

C0MBDY 

C0MMAT 

FRESP 

HEAD 

HY (B10ck Data) 

HYDRDY 

L0CFAC 

N0ZADM 

RKTDIF 

RKTZ 

RKZDIF 

S0LVW 

STEADY 

TADAMS 

TDPL0T 

XIMAGF 

ZADAMS 

ZER0 


C-1 

C-12 

C-I3 

C-20 

C-24 

C-27 

C-32 

C-34 

C-48 

C-49 

C-50 

C-67 

C-68 

C-75 

C-77 

C-79 

C-81 

C-94 

C-98 

C-103 

C-1 05 

C-1 08 

C-112 


R-9808/C-1 
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